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Dynamic Problems Associated With 
Satellite Orbit Control 


Many future satellite applications will require provision for changing the orbital 
parameters of the satellite once it is in an initial orbit. Such changes of orbital parame- 
ters can range from simple ‘‘orbit trimming,” to correct inttial guidance or propulsive 
3 inaccuracies, to relatively large-scale changes in the orbit for maneuvering purposes. 

The magnitudes of impulsive applications of force required to perturb the orbit by 
desired amounts is calculated for some typical cases, under the assumption of a central 
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force field. From these results, the case of continuously applied forces may be derived. 
The treatment is perfectly general, and permits analysis of the effects of arbitrary 
variations of an arbitrary initial conical motion. 


For a few typical cases, the corresponding requirements of propellant consumption 


| Introduction 


HEN A BODY moving in a conical orbit (that is, a 
circular, elliptic, parabolic, or hyperbolic orbit) in a central force 
field is subjected to the action of impulsive forces, the effect is to 
produce instantaneous changes in the elements of the orbit. 
Thus, in the case of elliptic motion about a focus, the elements 
can be taken as the length of the major axis, the eccentricity, the 
position of the apse line in the plane of the orbit, and the inclina- 
tion of the plane of the orbit to a fixed plane in space. The plane 
of the orbit and the fixed plane will intersect in a line through the 
focus—the line of nodes of classical astronomy—and the position 
of the apse line can be defined in relation to the line of nodes. In 
classical astronomy we further define the direction of the node 
line in relation to some fixed direction, and specify the time of 
passage of the body through apogee or perigee, thus furnishing 
the six elements of classical astronomy. 

Tn turn, these orbital characteristics are fixed, for a given cen- 
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Nomenclature and Definitions- 


of chemical rockets is shown. 


tral force field, by specifying, at some time ¢, the body’s initial 
conditions of distance from the focus of the central force field, the 
magnitude of the velocity, and direction of velocity, the plane of 
the motion being determined by the plane through the radius 
vector and velocity vector. The characteristic of the application 
of impulsive forces is that instantaneously the radius vector does 
not change, while the velocity, in either magnitude or direction, 
or both, does change. The changes in velocity vector then change 
the elements of the orbit. 

Our first problem is to find the relationships between these 
changes. The second problem is, given the changes in the magni- 
tude of the velocity, to find the corresponding propulsive require- 
ments; for example, the propellant consumption. 

With the idealization of impulsive thrust applications, this 
second problem is fairly trivial compared with the first problem. 
It is possible to dispense with the idealization of impulsive 
changes in velocity, but few simple, general results can then be 
described; and, in any event, the gross features of the motion are 
adequately exhibited by this concept of impulsive changes, al- 
though specific practical cases can require that this idealization 
be dispensed with. 

The accelerations during the application of thrust are always 
finite, of course; in cases where these accelerations are low, the 
body may move an appreciable angular distance during the thrust 


gravitational acceleration on 


r apse = the line passing through an 
surface of earth = 32.199 lial? dimensionless radius line apse and the earth’s center 
ft/see for central field as- vector € = eccentricity of a conical orbit 
sumed a = semimajor axis of ellipse, or (e = 0 for a circular orbit, 

Ry = radius of earth = 20,925,500 real principal semiaxis of 0 < € < 1 for an elliptical 
ft hyperbola orbit, € = 1 for a parabolic 

V; = velocity of earth satellite at orbit, and 1 fora 
surface of earth = (goRo)'/? a= ." dimensionless semima- hyperbolic orbit) 
= 25,960 ft/sec jor axis of ellipse, or dimen- perigee = minimum distance from 

V = velocity of body in arbitrary sionless real principal semi- center of force to orbit on 
point in conical orbit axis of hyperbola apse line; = a(1 — e) for 

~~ eins vedi & — apse = point on an orbit at a mini- elliptical orbit. For a hy- 

7 A mum or maximum distance perbolic orbit the corre- 
bh = goRo? from the center of force sponding distance = a(e — 
r = radius vector from center of (the center of force is the 1) 


earth to arbitrary point on 
conical orbit 


case) 
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earth’s center in the present 


apogee = maximum distance from cen- 
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application, and the idealization of impulsive thrust departs 
more and more from the actual case. In such circumstances the 
total period of thrust application may be divided up into several 
smaller periods of undisturbed motion separated by impulsive 
thrusts. The body then moves on a succession of instantaneous 
conical orbits, each of which may be found by the method of 
Section II. This procedure is known to the astronomers as the 
“variation of the elements’; in the limit this procedure leads to 
the case of differential equations for the changes in the elements 
when continual disturbing forces are present. Practical at- 
tempts to mechanize the guidance system required to effect orbi- 
tal changes will almost certainly require computer solution of the 
differential equations, for example. 

Several techniques are possible for treating the first problem. 
The explicit perturbation analysis of the effects of impulsive forces 
on the elements of, for example, an elliptic orbit goes back at 
least as far as Lagrange. In this treatment, which is the standard 
one, a basic orbit is given, and the small effects on the orbital ele- 
ments of small changes of the velocity vector are found by a per- 
turbation treatment. A very clear exposition of this case is given 
in Ramsey’s text [1]! as well as the limiting case of continuous 
force applications (and the associated differential equations); this 
latter can also be found in Moulton’s book [2]. Such an analysis 
considers only first-order effects; for example, the changes in 


A 
orbital elements are found in terms of the small parameter Vv? 


where AV is the change in the original velocity vector WV. For 
many cases such an analysis may be inadequate; for example, in a 
highly elliptical orbit the change in velocity possible near apogee 
may be a multiple of the original velocity and yet permit elliptical 
motion; also, it may be necessary to consider the change from 
elliptical to hyperbolic motion, ete. 

Rather than use the perturbation method of approach to 
variations in conical motion, it is possible to use the exact solu- 
tions for conical motion with perfectly general initial conditions. 
In this case it is clear that the exact effects of arbitrary variations 
in velocity vector can be considered by using the new velocity 
vector as a new initial condition. To attempt this purely ana- 
lytically is extremely cumbersome, but it can be done simply and 
conveniently by graphical methods. The approach then is 
simply to plot generalized charts of the properties of conical mo- 
tion, and consider arbitrary variations of orbital parameters in 
these charts. The construction of such charts is given in Section 
II. 

It is interesting to note that this second technique goes back to 
Newton [3]. In Book I, Prop. XVII, Prob. TX of the ‘Principia”’ 


' Numbers in brackets designate References at end of paper. 


Nomenclature 


he gives a geometrical construction for the general solution of the 
problem, which we may call Newton’s planar problem, of finding 
the orbital characteristics of a body started in a given position in 
a given central force field with given velocity and direction. In 
Cor. III he describes the application of this solution to the case of 
impulsive changes of the initial orbit, while in Cor. IV he de- 
scribes the case of continual changes of the initial orbit. Al- 
though Newton’s solution as given involves geometrical construc- 
tions, it is quite certain that he had the analytical counterpart of 
his geometrical solution also. In the meantime, analytical par- 
tial solutions of this problem have appeared from time to time in 
the recent literature [4, 5], and some attention has also been given 
to the problem of transfer orbit determination in simple cases [6, 
7, 8,9, 10, 11]. A very simple treatment of the general problem 
is described in Section IT. 

The treatment of Section II is quite general and permits arbi- 
trary conical orbits to be treated. Among the classes of problems 
to which this treatment is applicable are the following: 


A With given initial conditions, describe the motion of the 
body. This is simply the solution of Newton’s planar problem. 
B Orbit trimming. These cases may be of two kinds: 


1 Correction of errors in desired orbital elements brought 
about by guidance or propulsive inaccuracies. 


2 Correctionof earth oblateness effects or other perturbation 
effects. For example, earth oblateness generally rotates the in- 


stantaneous plane of the orbit and rotates the instantaneous line 
of apses [12], except at very specific orbit inclination angles. It 
may be required to find the impulsive thrust vectors periodically 
required to approximately maintain the initial orbital plane and 
direction of line of apses. 


C Maneuvering between arbitrary conical orbits, with no 
time constraints. 

D_ Interception of bodies moving in arbitrary conical orbits by 
other bodies which must be impulsively deflected into transfer 
orbits from other arbitrary conical orbits. It may be required to 
find the magnitude of the impulsive deflections required to effect 
interception. Another problem of this nature is the impulsive 
deflection of a body to evade a pursuer attempting interception. 

E Rendezvous of two bodies moving initially in two arbitrary 
conical orbits. The term “rendezvous” implies that the two 
bodies are finally in contact with zero relative velocity. 


The distinction between D and E occurs because in EF the planes 
of the orbits and the final velocities at the end of the transfer orbit 
must be matched, whereas one or both of these conditions may be 
absent in D. 


ter of force to orbit on apse 


line; = a(1 + €) for ellip- = » defined by equa- 
tions (4) 


tical orbit 
6 = polar angle between apse line 


and velocity vector; 8 is 
measured clockwise from 
the radius vector to the 
velocity vector; Bapse = 90 


and radius vector to arbi- V. = local circular velocity = deg 

trary point on orbit. @ is gokey\ y = 6B —-—7/2 

measured counterclock- = angle between two orbital 

wise, and 6 = O deg at 25,960 Vv, planes 

perigee ft/sec I.) = specific impulse of propellant, 
t,t, = time in orbit between 6 = 0 , p in seconds 

deg and = @ deg. The Vv m, = mass of body 

time from 0 deg to 6 deg ”% = V. = the dimensionless ratio of my, = propellant mass required to 


equals the time from 360 — 
6 to 0 deg 


velocity 
point on the orbit to local 
circular velocity at that 


vitrary 
at an arbitrary give a body of mass m, an 


impulsive velocity incre- 


ment of AV ft/sec 
, defined by equa- point 
o = —", defined by equation (9) 
tions (4) 8 = angle between radius vector Mm; 
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Il Analysis 


We begin with three basic equations: 
The energy equation: 


++) 
r a 


with the + sign for hyperbolic motion, and the — sign for elliptic 
motion. 


The parametric equation for a conic section: 


(1) 


a 1 + €cos 0 


with same sign convention. 
The equation expressing the constancy of angular momentum: 


r = [+ya(e? — (3) 


Equation (3) can be immediately integrated to find ¢ = (8, ¢, 


bt, a), and we get: 
3 1 1 
(4) on tan tan 9 6 


Pat 1 + €cos 0 


(4) 


ee? — 1)'/? sin 


1 + €cos @ 


+ + (e — 1)" tan 6 


—In 


(e + 1) — (e — 1)” tan 


These equations are plotted on Figs. 1, 2. From these, we can 


10 


Fig. 1 Dimensionless time parameter for elliptical orbits 
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obtain the time required to travel between two values of @ on 
elliptic or hyperbolic orbits. 


2a 
From equation (1), we have V? = u (==), consider now 


the elliptic case and define a velocity: 


which occurs for r = a (that is, at the ends of the minor axis in 
elliptic motion). This is the local circular velocity at the value of 
radius vector r = a. We then have: 


(la) 


r 
v 


From this we get: 


Elliptic motion: = 2 — p?, (v2 < 2) 


(5) 
= = : (v? > 2) 


Hyperbolic motion: 
r vy? — 2 


| 


= yp? — 2, 


These results are plotted in Figs. 3, 4. If we are given initial 
values of v and r, these equations determine a, for a given undis- 
turbed conical motion. Once the value of a is thus determined, 
the equations also give the relationship between v and r for the 
remainder of the undisturbed conic motion. 


Now, from the relation tan 8 = —r a we have immediately 
r 


on using (2): 


+e cos 


(6) 


This equation holds for 0 < €<1,1<e< (Bis the angle 
between the velocity vector and radius vector). By using the 


simple trigonometrical identities for sin 8, cos 8, in terms of tan 8, 
and, on using (2) and (5), we also have: 


UM 


SS. 


os 
aah ae 
AA 80 120 160 


Fig. 2 Dimensionless time parameter for hyperbolic orbits 
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Fig. 3 Curves for determination of a and r for elliptical orbits 


e? = cos? B + (v? — 1)? sin? B, (7) 


which again holds for both elliptic and hyperbolic motion. 

These forms of the general relations between 6, €, v, and @ are 
selected because of their adaptability to numerical solution. 
There are, of course, a great many other forms of such connecting 
relations. 

Equations (6) and (7) are plotted simultaneously to give Fig. 5 
for € < 1 and Fig. 6 fore > 1. 

Equations corresponding to (1)-(7) can also be derived for para- 
bolic motion, if desired. No use will be made of parabolic motion 
in the present paper, however. 


180 


Fig. 4 Curves for determination of a and ; for hyperbolic orbits 


The equations (5), (6), and (7), and the corresponding figures, 
give the general exact solution of Newton’s planar problem. 
Thus, if we launch a body at a given value of r with given V and 
given 8, the subsequent motion is entirely determined by the 
following procedure: 


Vv 
a From (la) find V,; then we know vy = —~. 


b From (5) find a (which must be positive). 

c From the value of v and the given value of 8, there is a 
unique value of € and @ which is found from (6) and (7); the 
value of @ uniquely determines the position of the apse line. 

d Starting from these initial conditions, the undisturbed mo- 


160 


(A 
7 


| 


Fig.5 Relation between 8 and v, with « and 6 as parameters, for elliptical orbits. 


€ values for the nine 


inner closed curves are, reading inward, « = 0.30, « = 0.25, « = 0.20, « = 0.15, « = 0.10, « = 0.08, 


€ = 0.06, « = 0.04, « = 0.02. 
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éfound inc. At any point of the € = constant curve, we can then 
find the value of 8 and y, from (6) and (7), and the value of r from 
(5). In particular, we can in this way find the conditions at the 
end of the transfer orbit. 


e From b,c, and (4), the times taken between any two values 
of 6 during the description of the orbit can be found. 


Conversely, suppose that at some point given by 7, 9 of the 
undisturbed motion, the body receives an impulse which changes 
its local velocity coplanarly from Vo to W’ and its direction from 
Bo to B’; the same procedure as before, using now as initial values 


ro, VW’, B’ 
then defines ‘he new orbit, as characterized by new values: 
the difference 6 — 4 gives the rotation of the apse line due to the 


impulse, etc. The magnitude of the required velocity increment 
is 


AV = {Vo2 + — cos (B’ — Bo)} 


The components of this, in the plane of the orbit, tangent to and 
normal to the original orbit, are VW’ cos (B’ — Bo) — Wo and V’ sin 
(B’ — Bo), respectively. One can also measure the components of 
AV (again in the plane of the orbit) in directions perpendicular 
and parallel to the radius vector; these components are: 


AV, = V’ cos (B’ — 1/2) — Vo cos (Bo — 1/2) 
and 
AV, = V’ sin (8’ — 7/2) — Vosin (Bo — 
respectively, or simply: 
AV, = V' cos y’ — Vo cos yo, AV, = V’ sin y’ — Vo sin Yo} 
(y=B 2/2). 


Next we must treat the case of deflections of the orbital plane 


tion is characterized by a constant value of €, namely, the value of 


itself. Any two orbital planes will intersect in a node line passing 
through the focus of the central field, and this is the line on which 
we will consider effecting transfer between two initially different 
orbital planes when it is desired to make the final orbital planes 
coincide, as in the rendezvous problem (if rendezvous is desired 
at an arbitrary time, two rotations of the orbital plane can be per- 
formed, at the beginning and end of the transfer orbit; this elimi- 
nates the waiting time required to reach the initial node line). In 
cases where only interception is required, a single rotation of the 
orbital plane can take place at any time about the instantaneous 
radius vector as an axis. In either case, call the initial angle be- 
tween the two planes ¢. Now, in general, we must satisfy 
several conditions; not only must the orbital plane be rotated by 
¢, but also in the new orbital plane the orbital elements may 
have to be varied to satisfy supplementary conditions, such as the 
body’s passing through a given point in the new orbital plane at 
a given time. Thus, suppose that in the original plane the ve- 
locity is Wo and the direction of this velocity is specified by Yo = 

Bo — 2/2, while in the new plane, rotated by @ degrees, the 
velocity and direction are specified by V’, y’ in the new plane. 

We can then characterize the required velocity increment by AV, 

and its three components by AV,, AV,, AV,, where AV, = 

component perpendicular to the radius vector, AV, = component 
parallel to radius vector (both in the original orbit plane), and 

AV,, = component normal to the original orbit plane; AV,, 

AV,, AV,, are therefore orthogonal. The defining relations are: 


AV = {V.? + V2 — 2V.V’ [sin y’ sin Yo 
+ cos cos Yo cos (8a) 


AV, = V’ cos y’ cos @ — Vo cos Yo (8b) 
AV, = V’ sin y’ — Vo sin Yo (8c) 
AV,, = V’ cos y’ sin (8d) 


(8b) to (8d) then specify the direction in which the impulse must 


be applied. 
One other simple result may be of interest. Suppose it is de- 


sired to simply rotate the orbital plane without changing the 
orbital elements in the plane. Consider the plane normal to the 


120° 
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Fig.6 Relation between 6 and v, with « and 0 as parameters, for hyperbolic orbits 
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radius vector; then the projection of the original velocity vector 
on this normal plane is 


Vo sin Bo = Wo cos Yo 


The magnitude AV, of the incremental velocity in this normal 
plane which, when added to the original velocity in the 
normal plane Vo cos Yo, produces a resultant of magnitude Vp cos 
Yo, but rotated by an angle ¢ is then: 


AVr = (Vo cos Yo) (2 sin $) (8e) 
Note that AV, is not normal to the original orbital plane. A 
number of interesting consequences may be drawn from (8e); 


thus, for example, for a satellite in a circular orbit characterized 
by 


=0 

near the earth, and provided @ is large enough, it may actually be 
more economical to first enlarge the orbit to one of large € = €’ 
and a = a’, rotate the orbit plane at apogee, and shrink the en- 
larged orbit back to the original value of a = a, € = 0 at perigee, 
rather than try to rotate the original orbit plane with unchanged 
orbit elements. 

By use of (8a) to (8d) we have a general technique for treating 
arbitrary three-dimensional changes involving all orbital elements 
simultaneously—for, given the vector change in the original 
velocity vector, we can determine a set of new orbital elements, or, 
given the initial and desired new orbital elements, we can deter- 


sion, AV is the sum of all the individual changes in velocity. The 
mass of the vehicle before the maneuver is then m,(1 + @). 

In many cases of practical application of these results to bodies 
moving in the earth’s gravitational field, it can prove convenient 
to introduce as fundamental constants Ro, the radius of the earth, 
go, Surface acceleration of gravity of the earth, and V,, the circular 
velocity at the surface of the earth. These constants are con- 
nected by the relation: 


Vv, = = 25,960 ft/sec 


We can then write: 


r 
= 


Ry 


=u; 


Therefore: 


t, = see = min, 
t, = see = min. 


Some Examples 


The results of Section II, and the figures associated with those 
results, are intended as basic working tools for treating a large 
class of problems. By the use of these results any problem of 
three-dimensional transfer orbits between arbitrary conical or- 


mine a required vector change in the original velocity vector to bits, which may be noncoplanar, can be treated, under the as- 


make the transition from the initial to the desired new orbital 
elements. An example of this is given as Example 3 in Section 


sumptions that the impulsive approximation for the deflecting 
forces can be used and that the motion takes place in a central 
field. The charts, if drawn to a larger scale, can give entirely ac- 


Finally, we can assess the propellant consumption required to ceptable accuracies for defining the major features of the original 


effect the changes in orbital parameters. For, given the required 
change in the magnitude AV of the velocity V, we can write for 
the ratio o of propellant mass consumed, m,, to the final mass of 
the body after completion of the required maneuver, m_ the 
equation: 


(9) 


where /,,, is the specific impulse of the propellant. Equation (9) is 
plotted as Fig. 7. If several maneuvers are performed in succes- 


and disturbed conical orbits. It is also possible to machine pro- 
gram the procedures governed by the equations of Section II, if 
high accuracy is desired. 

Because the variety of problems which can be treated is so ex- 
tensive, we will content ourselves here with some simple illustra- 
tive cases which can serve to outline the routine application of the 
equations and charts of Section II; the procedures involved can 
then be readily applied to more detailed and specific problems. 
Example 1 A Case of Newton's Problem 
A body is launched with the following initial conditions: 


7 


AV (FT/SEC) 


0.005 0 


0.05 O. 


05 
c= 


Fig. 7 Relation between velocity increment and dimensionless propellant consumption parameter 
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vy=1.2, p=3, B = 110 deg; 
therefore, 


a = 5.35, 6 = 60 deg, ¢€ = 0.53. 


Thus the apse line is 60 deg behind the initial radius vector. 
Suppose now that the body continues on this undisturbed orbit, 
characterized by € = 0.53, until @ has attained the value of 160 
deg, at which time it is desired to enlarge the orbit to a value of 
a = 6.5 without rotating the apse line. What must be the im- 
pulsive velocity increment, in the orbit plane, which must be 
given the body to accomplish this? 

The condition of no rotation of the apse line is a constraint 
which shows that only such combinations of 8 and v are possible 
as lie on the 9 = 160 deg line in Fig. 5. Just before the impulse 
we have vy = 0.75, 8 = 110 deg from Fig. 5; from Fig. 4, we 
then get p = 7.6. From Fig. 4, we then get the value of v re- 
quired after the impulse, to make a = 6.5, asy = 0.91. 

The value of 6 associated with @ = 160 deg, vy = 0.91 is B = 95 
deg, and the new € = 0.19. The magnitude of the required ve- 
locity increment Av = 0.27 and has components Ay, = 0.13 and 
Av,, = 0.23 tangent and normal to the original orbit at 9 = 160 
deg. The absolute magnitude of AV is 


AV = 0.27V, = (0.27)(25,960)(7.6)-"/2 = 2500 ft/sec 


From this we can use Fig. 7 to find the propellant consumption; 
for example, if J,, = 400, we get o = 0.22, 

Tf, on the other hand, we had required the velocity increment to 
be tangential to the orbit at 9 = 160 deg, then to get the same 
orbit enlargement to a = 6.5 at 8 = 110 deg = constant, we 
would have after the impulse 6 = 136 deg, thus showing that 
now the apse line is rotated by 24 deg, while the new € = 0.38. 


Example 2 A Case of Coplanar Rendezvous 


Two satellites are in the same circular orbit in the same orbital 
plane, but one is leading the other by A@ = 30 deg. For defini- 
tiveness, suppose that the initial location of the pursuer is at 9 = 0 
deg, while the initial location of the pursued is at 8 = 30 deg, and 
that @ = p = 6, and the initial value of 


v=1.0, 
B = 90°, 


corresponding to the circular orbit. ; 

It is desired to rendezvous; what velocity increments are neces- 
sary to achieve rendezvous locations of @ = 45, 180, 360 deg, while 
the pursued body moves through angles of 6, = 15, 150, 330 deg, 
respectively? (Since this is a rendezvous problem, two velocity 
increments are necessary, one to get the pursuer out of the circular 
orbit into a transfer orbit in the same orbital plane, another to 
get out of the transfer orbit into the circular orbit once again.) 
The transfer orbits are symmetrical in the present instance, so 
that the apse line of the transfer orbit bisects the rendezvous angle 
in each case; in this event, the two velocity increments men- 
tioned previously are also equal. The procedure then is as fol- 
lows. Consider the case where the rendezvous location is at 0 = 
45 deg; then, since the transfer orbit is symmetrical, the pursuer 
in the transfer orbit must initially be at a @, value of #, = 360 — 
45/2 = 3371/. deg, move through @, = 360 or 0 deg as he crosses 
the apse line of the transfer orbit, and intersect the circular orbit 
at 0, = 22!/.deg. The subscript identifies the values as applying 
to the transfer orbit. Therefore, select a series of €, values on the 
6, = 3371/2 deg line, and corresponding values of v,, and by using 
Fig. 4 obtain the corresponding values of a, the initial value of 
p = 6 being known. Use of Figs. 1, 2, then gives a series of values 
of ¢ for the transfer orbit, since a series of values of a, is known, 
and 6, goes from 3371/2 to 22!/. deg. The time for the pursued 
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body to move from 6 = 30 deg to the rendezvous location of @ = 
45 deg is 


(806.1) (4) (6'*) see = 3100 sec = 51.7 min. 


The value of €, for which the transfer orbit takes 3100 seconds is 
found to be €, = 6.9; and so v, = 2.89, 8B, = 70 deg, at this value 
of €, on the 0, = 3371/. deg line, thus defining the initial condi- 
tions of the transfer. (Note that a more general set of solutions 
can be generated by following this procedure, but then for 
each value of ¢ found computing what the initial position of the 
pursued body must be for rendezvous—that is, determining @,. 
A systematic exploration of the rendezvous problem can be ac- 
complished in this fashion.) During the transfer orbit a constant 
value of €, = 6.9 is maintained, so that values of p, v, 8, ete., dur- 
ing the transfer can be found. For this particular case, the trans- 
fer orbit is hyperbolic, and the magnitude of the original velocity 
increment is Av, = 1.98; this same (in absolute magnitude) 
velocity increment occurs at the end of the transfer orbit to match 
the velocity of the pursued body in the circular orbit, so that the 
total rendezvous requirement is 2Av, = 3.96. 

Proceeding in this way, we can construct Table 1: 


Table 1 

6 = rendezvous angle........ 45° 180° 360° 
6) = pursued’s angular travel... 15° 150° 330° 

t = time to rendezvous, sec.. 3100 31,000 68 , 000 
2Asse = 

= 26.0 0.50 0.048 


o = — for I, = 400....:... 
my 


This example shows how very rapidly the propellant require- 
ments increase as the time to rendezvous is decreased. The 
propellant requirements can be made arbitrarily small by per- 
mitting the pursued body to make multiple revolutions before 
rendezvous, but at the expense of greater and greater times 
to rendezvous. 


Example 3 A Case of Noncoplanar Interception 

A first body (the pursuer) travels in an elliptic orbit in the 
equatorial plane. To fix the angular position of the body we will 
use the geographical co-ordinates of the intersection of the radius 
vector with the surface of the earth, which is assumed to be non- 
rotating. The major axis of the ellipse is characterized by the 
value of a, = 2, while €¢, = 0.3, and the apse line at perigee is at 
latitude 0 deg, west longitude 90 deg. At the time ¢ = 0 the 
pursuer is at perigee, moving east, and we find »; = 1.14, p: = 1.4, 
8B, = 90 deg. In the meantime, a second body moves northward 
on a polar circular orbit with a2 = p: = 2 and longitude 0 deg; 
the value of vp = 1.0 and 8, = 90 deg = constant. At ¢ = 0 the 
second body is at latitude 0 deg moving north. It is desired to 
intercept this body at 30 deg north latitude, 0 deg longitude, at 
a time ¢ = (806.1)(2%/2)(7/6) = 1200 seconds, by giving the first 
body the proper transfer impulse at ¢ = 0. What are the charac- 
teristics of this transfer impulse? 

In general, the angular travel A@ of the pursuer required to 
effect interception must be found from the formulas of spherical 
trigonometry, but in the present case this angular travel is clearly 
90 deg. The problem, therefore, is to find some conical transfer 
orbit on which an angular travel of 90 deg requires 1200 seconds, 
and for which p; = pinitiat = 1.4, Py = Prinat = 2.0. The position 
of the apse line of the transfer orbit is, for the moment, undeter- 
mined, as is its eccentricity, so that we must, for example, re- 
peatedly select a value of €, and attempt to select the position of 
the apse line so as to effect interception, or vice versa. However, 
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since @,, €, are constant during the transfer orbit; the requirement 
that the transfer orbit pass through the values of p; = 1.4, py = 
2.0, used in (2), gives a relationship between the initial and final 
values of 0, (8;, 8,); in the present case this relationship is: 


Selection of a series of values of 0; then gives 6, = 6; + 90 deg; 
for each value of 6; we compute ¢, from the preceding equation. 
Each set of values of 6,, €, then gives a value of v, and 8, from Figs. 
5 and 6, and a value of a, from Figs. 3 and 4. We can then com- 
pute a@,'/? and, knowing 6,, 6,, a,"*, the times for the traversal of 
the 90 deg conical arc, from Figs. 1 and 2. 

We find, in this way, that interception is effected by choosing 
the approximate values 


6; = 330°, 6, = 60°, e, = 2.35, 
vy, = 1.88, a, = 0.93, B, = 69°. 
The perigee distance of this hyperbolic transfer orbit is then 
Pperigee = 1.26, 


or a perigee altitude of 0.26 earth radii above the surface of the 
earth. This set of values completely specifies the orbit elements 
in the plane of the transfer orbit, which had to be rotated by ¢ = 
30 deg out of the equatorial plane. 

Introducing the notation of Section II, we have: 


V. = 1.14 
vy’ = = 1.88 
Yo = 0° 
y’ = -21° 
= 30° 
Use of (8a) to (8d) then gives: 
Av = 1.17 
Av, = 0.38 ; 
Av, = —0.67 
Av, = 0.88 


These values completely specify the magnitude and direction 
of the impulsive force required for the pursuer to effect intercep- 
tion. 

Finally, we have: 


AV = (Av)(V,) = (1.17\(25,960)(1.4)~'/? = 25,600 ft/see 


Use of Fig. 7 for J,, = 400 then gives for o the required value of 
o = 6.2. 
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In a similar fashion we can treat the general classes of problems 
relating to three-dimensional conical transfer orbits between two 
arbitrary initial conical orbits with specified boundary conditions. 
As experimentation with the procedures will show, it is not al- 
ways possible to find real solutions to interception or rendezvous 
problems with fixed initial and terminal conditions; for example, 
the required transfer orbit may pass through the earth, ete. Such 
cases may require the terminal conditions to be relaxed, or other 
modifications to the problem to be made. 

The same techniques may also be applied to cases where the 
earth is not the focus of the central force field. For example, 
transfers between planetary orbits may be treated, taking full 
account of orbital eccentricity, inclination to the plane of the 
ecliptic, varied positions of the apse line, etc., and such specialized 
problems as the calculation of the velocity increment induced in a 
vehicle by the attraction of a massive moving body [13]. For such 
problems, and for reasonably accurate solutions of problems such 
as the preceding ones, the figures must, of course, be drawn to a 
larger scale than the scale in which they appear in the present 
paper, and the value of 4 must be appropriate for the central 
force field under consideration. 


References 


1 A.S. Ramsey, “Dynamics,” Part II, second edition, Cambridge 
University Press, London, England, 1945, pp. 14-19. 

2 F.R. Moulton, ‘‘Celestial Mechanics,’’ second revised edition, 
The Macmillan Company, New York, N. Y., 1914, pp. 366-406. 

3 F. Cajori, ‘Newton's Principia,” University of California Press, 
Berkeley, Calif., 1946, pp. 65-67. 

4 K. Shutte, “Die Bahnbestimmung aus dem Vektor der Bahn- 
geschwindigkeit und der Einfluss einer Aenderung desselben auf die 
Bahnelemente,” in ‘‘Space-Flight Problems,’’ Laubscher & Cie, Biel- 
Bienne, Switzerland, 1954, pp. 89-102. 

5 M. Vertregt, ‘‘Die Bahnbestimmung aus dem Vektor der 
Bahngeschwindigheit,’’ Astronautica Acta, vol. IV, fase. 2, 1958, pp. 
135-137. 

6 B. H. Paiewonsky, ‘Transfer Between Vehicles in Circular 
Orbits,” Jet Propulsion, vol. 28, February, 1958, pp. 121-123. 

7 OM. Vertregt, ‘Interplanetary Orbits,” Journal of the British 
Interplanetary Society, vol. 16, no. 6, March-April, 1958, pp. 326-354. 

D. F. Lawden, “The Determination of Minimal Orbits,’’ 
Journal of the British Interplanetary Society, vol. 11, no. 5, September, 
1952, pp. 216-224. 

9 H. Preston-Thomas, ‘‘Generalized Interplanetary Orbits,’’ 
Journal of the British Interplanetary Society, vol. 11, no. 2, March, 
1952, pp. 76-85. 

10 D. F. Lawden, ‘Fundamentals of Space Navigation,’’ Journal 
of the British Interplanetary Society, vol. 13, no. 2, March, 1954, pp. 
87-102. 

11 D. F. Lawden, “Corrections of Interplanetary Orbits,”’ Jour- 
nal of the British Interplanetary Society, vol. 13, no. 4, July, 1954, pp. 
215-223. 

12 D. G. King-Hele and D. M. C. Gilmore, ‘‘The Effect of the 
Earth’s Oblateness on the Orbit of a Near Satellite,”” Tech. Note No. 
G. W. 475, Royal Aircraft Establishment, Farnborough, October, 
1957. 

13. D. F. Lawden, ‘Perturbation Manoeuvres,” Journal of the 
British Interplanetary Society, vol. 13, no. 6, November, 1954, pp. 
329-334, 


Transactions of the ASME 


f 
Py 1+¢€,cos); 
Pi 1 + €, cos 0, 
vid 
. 
: 


H. TAKEYAMA 
E. USUI 


Research Engi s, Gov t 
Mechanical Laboratory, 
Tokyo, Japan 


Study of Machinability of Metals 


If the definition of the so-called machinability is simplified, a rather scientific treat- 
ment of machinability is possible. 
bility of a metal in terms of chip formation, which is one of the most important aspects, 
is closely related to the ratio of the shearing strain to the shearing stress near the fracture 
point of the shearing-test curve of the metal. 


The experiment to be described reveals that machina- 


The larger the ratio, the poorer the 


machinability from the viewpoint of chip formation. 


Wir 1s machinability different among various metals 
even though the machining conditions are entirely equal? 
For example, when copper is machined, the chips are very thick 
with a very small shear angle, whereas when an alloy steel is 
machined, the behavior is opposite. 

The object of this experiment is to investigate the reason 
for this. However, the term ‘“machinability’” has various 
aspects, that is, chip formation, cutting force, tool life, surface 
finish, and so on, which are associated with each other to some 
extent. In this paper the definition of machinability is con- 
fined to that in terms of chip formation, which is one of the 
fundamental aspects, in order to find some sensible relationship 
between machinability and the physical properties of a metal. 

In an earlier paper! by the authors, machinability in terms of 
chip formation was demonstrated to be controlled by the tool- 
chip contact area, and the equation cot @ = cot dy + KA/(A’ 
cos @) was introduced, where @ is the shear angle, @p the shear 
angle when the tool-chip contact area is assumed to be zero, A 
the tool-chip contact area, A’ the area of cut, a the rake angle, 
and K a constant of proportionality. Therefore the analysis 
of machinability thus defined must be possible by means of the 
theory of tool-chip contact area previously derived by the 

authors. 

In this experiment, the machinability data, when various 
types of metals were machined under definite conditions, were 
analyzed and compared with the material testing data for these 
metals. As a result, it was found that the slope of the shearing 
test curve near the fracture point is closely associated with the 
behavior of the metallic fracture ahead of the tool face during 
machining, and, furthermore, with the tool-chip contact area or 
the machinability in terms of chip formation. 


Test Procedure 


In order to investigate the relationship between machining 
characteristics and physical properties of metals, various types of 
metals such as aluminum, copper, brass, mild steel, alloy steel, 
Meehanite cast iron, and various grades of brass were machined 
under definite conditions, and then the foregoing metals were 
subjected to several types of material testing, such as compression 
tests, hardness tests, shearing tests, and model machining tests. 

Machining Test. Tubes of aluminum, copper, brass(Cu:60 
per cent, Zn:40 per cent), mild steel (C:0.2 per cent), alloy 
steel (C:0.44, Ni:1.76, Cr:0.77, Mo:0.19, Mn:0.73 per cent) and 
Meehanite cast iron, whose external diameter and wall thickness 
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were 2.8 in. and 0.12 in., respectively, were machined two- 
dimensionally without any cutting fluid. For each material, 
shear angle, two components of the cutting force, and tool-chip 
contact area were measured. The shear angle was obtained 
either from the chip thickness or the weight and length of the 
chips depending upon the shape of the chip section, and the 
tool-chip contact area was measured by a photomicrographic | 
observation of the tool face. Whether an observed tool-chip 
contact area is a real one or not was discussed in the previous 
paper.! 

For such a basic study of metal machining, the influence of the 
built-up edge should be minimized as far as possible, and as to 
chip formation the so-called flow type is desirable. These are 
the reasons why the machining tests were conducted at a high 
cutting speed (330 fpm). Carbide tools of a steel-cutting grade, 
which were prepared under definite conditions, were used for the 
turning tests, and the tool shape was 0, 0, 6, 0, 0, 0, 0. 

Plates (5 in. X 3 in. X 0.12 in.) of brass, whose zine content 
was varied such as 20.7, 30.5, 34.5, 36.6, 38.4, 40.2, and 42.2 
per cent, were machined at a very slow speed by using a milling 
machine, the type of cut being two-dimensional without any 
cutting fluid. Tools used in that case were of high-speed steel 
and the tool shape was 0, 0, 6, 0, 0, 0, 0. 

Material Test. Vickers or Brinell hardness, which was as- 
sumed to be the flow stress of the metal, was measured for each 
metal. 

A shearing test also was conducted for each material mentioned 
by using an Amsler testing machine for obtaining the relation- 
ship between the shearing stress and the shearing deformation as 
well as the shearing strength. A test piece, 0.08 in. diameter and 
2.0 in. length, was fitted closely into the hole of a jig as shown 
in Fig. 1, which was held between the jaws of the testing machine. 


1H. Takeyama and E. Usui, “The Effect of Tool-Chip Contact 
Area in Metal Machining,’”’ Trans. ASME, vol. 80, 1958, pp. 1089- 
1096. 


Fig. 1 Jig for shearing test 
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Fig. 2 Setup of model machining test 


From a physical point of view, shearing strain is more desirable 
than shearing deformation, but by such a simple shearing test it 
was impossible to determine the shearing strain. 

Model Machining Test. In a conventional machining test the 
shearing manner in the shear zone is greatly affected by the 
shear angle as well as by the physical properties of the metal. 
The purpose of the model machining test is to isolate the in- 
herent shearing pattern of a metal resulting from physical 
properties themselves, by removing the influence of the shear 
angle. Fig. 2 shows the setup of the model machining test, in 
which the test piece, having a projection, is fastened to the face 
plate of a lathe and turned so as to be sheared off at the shoulder 
of the projection, by a tool of zero rake angle, in the direction 
of Y-Y. Although, strictly speaking, the shearing in such a 
model test is not quite identical with that in a conventional 
machining test, the shearing fracture of the projection can be 
qualitatively compared to that in the shear plane during actual 
metal machining. 


Test Results 


As reported in the previous paper,! theoretically and experi- 
mentally the following equation was derived. 


COs a A’ cos @ 

where @ is the shear angle, A the tool-chip contact area, A’ the 
area of cut, @ the rake angle, and K a constant related to the 
ratio of the normal stress on the rake face to the shearing stress 
in the shear plane, or, the ratio of the frictional stress on the rake 
face to the normal stress on the shear plane, which is considered 
to depend upon the material cut. Therefore, as long as both the 
rake angle and the material cut are constant, the shear angle 
should be a function only of the tool-chip contact area per unit 
area of cut as understood in Equation (1). As a matter of fact, 
when cot ¢ is plotted with regard to A/(A’ cos @) for various 
metals including aluminum, copper, low-carbon steel, brass, 
alloy steel, and Meehanite cast iron, the relationship is linear 
within a certain amount of error as shown in Fig. 3. This 
means that AK in Equation (1) does not vary much with regard 
to the material. This fact enables us to discuss machinability 
in terms of tool-chip contact area per unit area of cut only to a 
first approximation. 

It is probable that the tool-chip contact area per unit area of 
cut increases when the “flow stress” of a material cut is small 
compared with the shearing stress of the material, if it is as- 
sumed that the normal load on the rake face, which is a function 
of the shearing stress on the shear plane, is equal to the product of 
the flow stress and the tool-chip contact area. 

Plotting the relationship between the tool-chip contact area 
per unit area of cut and the ratio of the shearing stress to the 
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COTANGENT OF SHEAR ANGLE, Cot @ 
a 


4 6 8 o 6 8B 
A/(ACo8 a) 


Fig. 3. Relationship between cotangent of shear angle and A/(A’ cos a). 
A,A’ and a denote tool-chip contact area, area of cut, and rake angle, 
respectively. Material Meehanite cast iron, Ni-Cr-Mo steel, brass, 0.2 
per cent carbon steel, copper, aluminum; tool tungsten carbide of steel- 
cutting grade; tool shape 0,C,6,0,0,0,0; cutting speed 330 fpm; feed 
0.010, 0.0046 ipr; type of cut two-dimensional; cutting fluid dry. 
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Fig. 4 Relationship between tool-chip contact area per unit area of cut 
and ratio of shearing stress to flow stress of material. Cutting conditions 
same as in Fig. 3. 
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Fig. 5 Shearing test curves for various metals 


flow stress for which the Vickers hardness is substituted, 
there results Fig. 4. Apparently it shows no consistent cor- 
relation between both variables. 

Since the metal removal ahead of a tool face in metal machining 
is caused mainly by shearing, the characteristic of the shearing 
test curve of a metal should be associated closely with the 
machinability characteristic of the metal. Fig. 5 shows the data 
of the shearing tests carried out with the setup as shown in 
Fig. 1 for aluminum, copper, low-carbon steel, brass, Ni-Cr-Mo 
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alloy steel, and Meehanite cast iron. Although shearing strain 
is more desirable as abscissa than shearing deformation for a 
physical analysis, it was impossible to obtain because of the type 
of shearing test. 

As seen in Fig. 5, each metal has its own shearing character- 
istic, especially near the fracture point. In the case of such 
metals as aluminum or copper the shearing fracture takes 
place accompanying a considerable increase of shearing de- 
formation without an appreciable increase of shearing stress, 
so that no one can tell exactly when the shearing fracture occurs. 
On the contrary, in the case of such metals as Meehanite cast 
iron or alloy steel, the shearing stress increases rather steeply 
compared with the shearing deformation until the point of 
shearing fracture is reached where the stress suddenly falls as 
soon as the fracture occurs, so that it is easily recognized. 

It is not unreasonable to consider that the shearing fracture 
along a shear plane in metal machining corresponds to the 
fracture point in the shearing test which is identical with the 
uppermost end of each curve in Fig. 5, and that the stress or 
strain pattern in front of a tool face, which will control the 
machinability characteristic of a metal, is closely associated with 
the behavior of the shearing test curve near the fracture point. 
Plotting the tool-chip contact area per unit area of cut, which 
is denoted by A/A’, for each metal against the slope of the 
corresponding shearing test curve near the fracture point, which 
can be found by the tangent of the curve at the fracture point, 
the relationship is seen to be parabolic as shown in Fig. 6. This 
fact demonstrates that the ratio of the shearing deformation to 
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Fig. 6 Relationship between tool-chip contact area per unit area of cut 
and shearing deformation per unit shearing stress near fracture point for 
various types of metal. Cutting conditions same as in Fig 3. 


Ni-Cr-Mo STEEL 


MILD STEEL ALUMINUM 

Fig. 7 Photographic observation of shear zone in model machining test. 
Upper part of each photograph shows the shear zone in model machining 
test. Material Meehanite cast iron, Ni-Cr-Mo steel, 0.2 per cent carbon 
steel, brass, copper, aluminum; tool high-speed steel; tool shape 0,0,- 
6,0,0,0,0; speed quasi-static. 
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Fig.8 Composite data of shear angle, tool-chip contact length and width 
of shear zone in model machining test. Cutting conditions same as in 
Fig. 3. 


WIDTH OF SHEAR os IN MODEL TEST, IN 
TOOL-CHIP CONTACT LENGTH. IN 


COTANGENT OF SHEAR ANGLE 


the shearing stress near the fracture point has a definite relation 
to the machinability characteristic in terms of chip formation. 
The larger the shearing deformation per unit shearing stress 
near the fracture point, the larger the tool-chip contact area 
per unit area of cut or the poorer machinability in terms of 
chip formation. This fact seems to mean that a metal 
whose shearing behavior is diffuse, produces a large tool-chip 
contact area in metal machining or poor machinability in that 
sense. From this point of view, it will be interesting to observe 
the proper shearing pattern of each metal. 

The upper part of each photograph in Fig. 7 shows the pattern 
of the shear zone for each metal which was obtained by the model 
machining test explained before. The maximum width of the 
deformed zone in the model machining test varies with regard 
to the material. In other words, diffusibility or sharpness of 
shearing slip, which is inherent for each metal, is associated with 
the slope of the shearing test curve near the fracture point. 
As seen in Fig. 8, which is the composite data showing the 
machinability characteristics for the afore-mentioned metals, 
the machinability in terms of chip formation which can be 
expressed by the cotangent of shear angle (cot @) has a good 
parallelism with the tool-chip contact area per unit area of cut, 
and the width of the shear zone in the model machining test. 

Now let us show another example of the afore-mentioned rule 
of machinability. The machinability characteristic of brass 
conspicuously varies with the content of zinc therein. Generally, 
brass having more zine is less sticky and results in a smaller 
tool-chip contact area or a larger shear angle in machining as 
seen in Fig. 9. Since the relationship between the cotangent of 
shear angle (cot @) and A/(A’ cos q@) is nearly linear as seen in 
Fig. 9, K in Equation (1) should be nearly constant for varying 
grades of brass. Therefore the machinability in terms of chip 
formation can be expressed, to a first approximation, only by the 
tool-chip contact area per unit area of cut as in the case men- 
tioned before. 

Fig. 10 shows the relationship between the tool-chip contact 
area per unit area of cut and the shearing deformation per unit 
shearing stress near the fracture point on the shearing test 
curve for various grades of brass, containing 20.7, 30.5, 34.5, 
36.5, 38.4, 40.2, and 42.2 per cent zine. The machining data 
were obtained by two-dimensional machining tests at a very 
slow speed and a rake angle of 20 deg, and the shearing test data 
were obtained by the previously mentioned method. The 
relationship is apparently parabolic as in the case of Fig. 6. 
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Fig. 9 Relationship between cotangent of shear angle and A/(A’ 
cos «) for various grades of brass. A, A’, and a denote tool-chip contact 
area, area of cut, and rake angle, respectively. Material various grades 
‘of brass (Zn: 20.7, 30.5, 34.5, 36.5, 38.4, 40.2, 42.2 per cent, Cu: rest); 
tool high-speed steel; tool shape 20,0,6,0,0,0,0; type of cut two-dimen- 
sional; cutting fivid dry. 


The numerical difference between Fig. 10 and Fig. 6 is due to the 
fact that the machining conditions, especially the rake angles, 
are different between the two tests, and, furthermore, some 
lubricating components such as zine and lead which have a ten- 
dency to decrease the tool-chip contact area will shift the curve 
downwardly as seen in Fig. 10. 

Generally speaking, tool-chip contact area or machinability in 
terms of chip formation is determined by two factors; that is, the 
characteristic of the strain or stress field within a material cut 
in front of the tool face and the affinity between chip and tool 
face. Therefore such metals as free-cutting steel or brass, 
which contain lubricating components, cannot be generalized 
by the foregoing rule only, since the affinity between chip and 
tool face must be considered as well, although the rule holds well 
within the same group of metal as seen in Fig. 10. 

Now let us consider the physical meaning of the relationship 
‘of Fig. 6. The assumptions for this purpose are as follows: 


(a) The shearing stress in metal machining reaches the maxi- 
mum attainable static value stress in the shearing test on the 
shear plane. 

(b) The stress situation in the shear plane in metal machining 
corresponds to the stress near the fracture point of a metal 
which is denoted by the uppermost end of each shearing test 
curve in Fig. 5. 

(c) The plastic zone within a material cut in front of the tool 
face determines the size of the tool-chip contact area. 


In the case of such a material as aluminum or copper, where 
the slope near the fracture point of the shearing test curve is 
shallow, the shearing slip or fracture of the part to become a 
chip occurs over a wide range across the so-called shear plane 
because of the large probability of so-called after-slip or fracture 
due to the diffuse shearing characteristic as shown in Figs. 5 
and 7. In other words, the shearing slip or fracture can take 
place over a wide plastic zone of the material ahead of the tool 
face because the stress near the fracture point can be regarded 
to be nearly constant over the wide range of shearing deformation. 
This apparently results in a larger tool-chip contact area per 
unit area of cut. On the contrary, in the case of such a metal as 
Meehanite cast iron where the slope near the fracture point of 
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Fig. 10 Relationship between tool-chip contact area per unit area of 
cut and shearing deformation per unit shearing stress near fracture point 
for various grades of brass. Cutting conditions same as in Fig. 9. 


the shearing test curve is steep as seen in Fig. 5, and, once 
fractured, the metal is separated infinitely, the probability of 
further shearing slips after the part to become a chip has passed 
through the shear plane along the tool face must be small. 
In other words, the shearing fracture or slip must be completed 
to a large extent at the shear plane, because otherwise, considering 
the steeper slope near the fracture point of the shearing test 
curve, the shearing stress of a larger part which has passed 
through the shear plane must be appreciably higher than that in 
the shear plane, and this apparently contradicts the assumption 
(a). At any rate, such a shearing characteristic makes the tool- 
chip contact area small. 

As understood by the foregoing explanation, the slope near 
the fracture point of the shearing test curve determines the 
probability for a differential element, which is to become a 
chip, to undergo further slips or fractures after it has passed 
through the shear plane, or the range of the stress field. The 
smaller the probability, the smaller the tool-chip contact area 
per unit area of cut, because if the probability is very small, 
the shearing slip or fracture must be completed mostly at the 
shear plane and, at the same time, the force of the deformed chip 
exerted on the tool face above the shear plane will be almost 
released. 


Every metal has its own machinability characteristics, and the 
machinability in terms of chip formation which can be expressed 
by shear angle or chip-thickness ratio has a good parallelism 
with the tool-chip contact area per unit area of cut. The tool- 
chip contact area per unit area of cut in metal machining is 
closely related with the slope near the fracture point of the shear- 
ing test curve. The larger the shearing deformation per unit 
shearing stress near the fracture point, the larger the tool-chip 
contact area per unit area of cut, because the foregoing slope 
can be considered to determine the probability for the shearing 
slip or fracture of the element to become a chip to take place 
after it has passed through the so-called shear plane, or the range 
of the stress field in front of the tool face. Thus shearing test 
data can be used for evaluating the machinability in terms of 
chip formation except when some lubricating components are 
contained within the material cut. 
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treatment. 


rR THEIR published discussions of Freudenstein’s 
excellent paper,? de Jonge, Kenngott, and Nickson indicated the 
desirability of developing a simple mathematical method for syn- 
thesizing a four-bar linkage, a method in which the physical rela- 
tionships are easily recognized in the mathematics. The following 
mathematical development for four precision points is based on 
Hein’s graphical solution.® 


Nomenclature 
The following nomenclature is shown in Fig. 1. 


i = 1,2,3,4 = subscripts which refer to precision points (or 
positions) 

lengths of links 

y = f(z) function to be generated 

a = f(8) = function in terms of angular positions 

1 ZB-A-C, measured counterclockwise from A-B 

to A-C 

y ZB-A-D, measured counterclockwise from A-B 
to A-D 

6; = the smaller angle between links a and d 


a, b, c, d 


The remaining notations are shown in Fig. 1. 


Graphical Solution 


The graphical solution will be considered briefly. The com- 
plete graphical solution, which is based on inversion, is shown in 
Fig. 1, where, for clarity, the angles between the precision points 
for the required function are shown separately at the upper left of 
the figure. The procedure for making a graphical solution is as 
follows: To obtain point-position-reduction (Punktlagenreduk- 
tion),* the @ range is oriented so that the angle between two arbi- 
trary precision points (in this case 3 and 4 were chosen) is bi- 
sected by the extension of link a. The 8-range need not be prepo- 
sitioned because only the magnitudes and relationship of the 
angles are required when inversion is made with link b stationary. 
For convenience, inversion is made about position 1 which is 


1 Based on a graduate project completed at Michigan State Uni- 
versity, College of Engineering. 

2 F. Freudenstein, ‘‘Approximate Synthesis of Four-Bar Linkages,” 
Trans. ASME, vol. 77, 1955, pp. 853-861. 

3 Kurt Hain, ‘‘Angewandte Getriebelehre,’”’ Hermann Schroedel 
Verlag, Hannover, Germany, 1952, pp. 362-363. 

Contributed by the Machine Design Division and presented at 
the Annual Meeting, New York, N. Y., November 30-December 5, 
1958, of THe AMERICAN Society OF MECHANICAL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, June 17, 
1958. Paper No. 58—A-130. 
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Four-Bar Linkages—Approximate Synthesis 


An analytical method for synthesizing a four-bar linkage as a function generator is 
presented. The method, which permits the arbitrary selection of four precision points 
and finite angular ranges, is based on a graphical solution. This permits a preliminary 
graphical investigation of the six possible linkages before selecting one for analytical 


shown horizontal. Links a and d for the inverted position cor- 
responding to precision point 3, are shown as A-H and H-E where 
the length of link d is not yet determined. Links a and d for the 
inverted position of precision point 4 are shewn as A-G and G-E. 
The intersection of these two inverted positions of link d is at E. 
From symmetry, G-E = H-E (this resulted from bisecting a; with 
the extension of link a). To utilize point reduction, we choose H-E 
for the length of d. With the length of d determined, the in- 
verted position for precision points 2 and 1 can be drawn as 
A-J-F and A-B-C, The perreadicular bisectors of C-F and F-E 
intersect at D which completcs the required linkage A-B-C-D. 


Derivation of Basic Equations 


The derivation of the equations will be based on the graphical 
solution, Fig. 1. After the precision points are assumed for one of 
the links and the corresponding points calculated for the other 
link, the a-symbol is assigned to the angular positions of the link 
covering the larger range and placed at B as shown. (In the work 
that follows, six possible solutions for each set of assumed pre- 
cision points will be indicated. If the @ and §-ranges are equal, 
one of the solutions will be indeterminate.) The 6-symbols are 
assigned to link 6 which has a smaller displacement range. S;, the 
diagonal A-C for the linkage in position 1, is at angle [’, with link a. 
Link 6 forms an angle W from the horizontal, and is measured 
positive in a counterclockwise direction. Since link b was fixed 
in the inversion process, Y is a constant. 6, is the angle A-B-C in 
position 1. @ is always the smaller angle between links a and d. 
The values of 6 are 


6, 


6: 


6, = 0 = 180 


Applying the law of cosines to G-A-H and G-E-H, and equating 
sides G-H, gives 


d? + d? — 2d? cosw = a? + a? — 2a? cos B; ( 


where w is obtained from the expression for the interior angles of 
quadrilateral A-G-E-H 


Bs + (180 = + (130 = +w = 360 


@ = a; — By 
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H = 180 + +) (1) 
or 


Graphical solution used for the derivation of equations 


Replacing » in Equation (2) by the foregoing value and letting 
a = 1, dis found to be 


1 — cos 
— cos (as — 


With links a and d and the angles 6,,.2.3,4 known, there are three 
remaining unknowns; namely, links b and c, and angle y. By 
noting the trigonometric relationships of these unknowns in posi- 
tions 1, 2 and 3, three equations can be written and solved simul- 
taneously. These equations are: 


Position 1, triangle C-A-D 
c? = b? + S,? — 2bS, cos (yy — T)) 
where (YW — T',) = angle C-A-D. 

Position 2, triangle F-A-D 
c? = b? + S,? — 2bS, cos [py — (8; + T:)] (5) 
where [WY — (8, + T,)] is angle F-A-D, and T, is the angle J-A-F. 

Position 3, triangle E-A-D 
c? = b? + Sj? — 2bS8; cos [W — (8: + + (6) 


where [Y — (8; + 8. + T)] is angle E-A-D, and T; is the angle 
E-A-H, 

Eliminating c? and 6? by subtracting Equation (5) from Equa- 
tion (4) gives 
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0 =0+ 8S? — — 2b{S, cos (W — T)) 
— cos [W — (8: + T2))} (7) 
Letting (8; + TI.) = @ and solving for b gives 
S;? — 
2[S; cos — Ts) — cos — @)] 


Subtracting Equation (6) from Equation (4) gives 
0 = 0+ S,? — S;* — 2b{S, cos — T)) 
— S; cos (6+ (9) 
Letting [8; + 8. + 3] = B and solving for b gives 
S,2 — S;? 
2[S, cos — S; cos (y — 


(10) 


Equations (8) and (10) are equated, giving 

S? — S,? 

2[S; cos — — S: cos — @)] 
2[S; cos — Ts) — cos — B)] 


(11) 


This can be written in the form 


reos(W +s cos(W — )=tecos(y — @) 


where 


(12) 
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Table 1 Numerical values for the six solutions of the example 


1-3 | | 2-3 | 2-4 | 3-4 


@.| 148° 45° |130° 130%00") 111915" | 98°45" 
—— |167°30'| 167°30' | 
| ——  |167°30' | —— _, |173°45' 
e,| | | 13000°} —— | —— 


d |0.3824 |0.5786 |i.0000 |0.9843 | 1.7/66 |2.6505 


$1 |1.3676 | 1.5245 |1.8126 | 1.7984] 2.2784 | 2.9718. 
—— 1.5764 |(.988! | }.9725 | 2.6233 | 3.4430 
$3} 1.3033 — | — 2.1129 | 3.6118 
S4|1-1236 | 1.3244 | |1.7984| —— | —— 


Q5965 |-1.1142 | 0.0000 |-1.1806 | |.0864 |3.5379 
|1-326! 1.1806 | 4.5870 }0-9177 
t |0.7923 | |!.326! | 0.0000 | 5.6886 |14.5064 


r 11°42.39 2 | | 6 19°$0.63' 32° 983, 

Ty |19°39,36' | © 15.0" 290735! 3°57.00' | 13°38.64 


0.36" 8° 1.27), |16°33.68 | 16 *30.72' |30° 9.31" 4228.51 
B }30°20.64' "50.37 |16°27.75' | 2512.69 |26°39,75" |36°2/.39 


-89°20' |16°13.22'| 205°0.17't 51°49,56 


b 10.6766 |0.7973 |1.6955 |-2.1524 | 1.3611 |1.5629 
[1.5726 |1.4739 |0-2900)0.353g |o 7818 |3.8196 


r= S,(S3? So?) 


t 


— S,?) 


When expanded, Equation (12) is 


r[cos cos T; + sin sin T;] + s[cos cos ® + sin W sin B] 
= t[cos cos @ + sin sin Q] 
. Collecting terms and solving for gives 


tsin @ — rsinT, — ssin® 


From Fig. 1, triangle C-A-B for position 1 is used for reference to 
write the genera! equation 


I, = ken (15) 


S,;2 = a? + d? — 2ad cos 6; where 1 =1,2,3 (16) 


and 


It is now possible to solve for b in Equation (10). With b known, 
c can be determined from Equation (4). 


Numerical Example 


Design a four-bar linkage to generate the function y = 2*}, 
where z is to range from 1 to 5 with precision points at 1, 2.5, 3.5, 
and 5. From the equation, the y-range corresponding to the 2- 
range is from 1 to 29.366. The a-range corresponding to the 
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S3(S2? — S,?) (13) 


Fig. 2. The six solutions of the example shown in position 1. 


z-range, is to be from 0 to 100 deg, and the B-range is to be 


from 0 to 50 deg Then 


Substitution of these values of x and y in the original equation 


gives 


For the precision points, the corresponding angles are 


The numerical values obtained in the solution of this problem are 
recorded in the right-hand column of Table 1, and the linkage, 
drawn in position 1, is shown at the lower right in Fig. 2. 


B = 1.7627 


37° 30.00’ 
25° 0.00’ 
37° 30.00’ 


Tmax — 5 - 1 
Qmax ) ad ( 100 ) sed 


max min 29.366 — 1 


25 


50 


Bi 
Bs 


B; = 27° 17.26’ 


10° 18.68’ 


12° 24.07’ 
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Table 2 Equations for the five alternate solutions 

1-2 2-3 a- 4 
(54-52) 5,(S,2-52) |$,(57-52) 6,(52-52) |5,(55-5) 
Su (53-5?) $4 5;(S2-S;) 5, (52-5) 5;(5£-5,7) 
[80-(=2=*) 190-4) 180-425) 


| 190-(—253) 


izt,3,¥ (= 42,4 1,2,3 t=1,2,4 i=1,2,3 


| 


2 [S. cos ¢y-a) 


4 
2 


-cos(4, +8; +A, ) 


SA 


Coan 

~ 


(B, +B, ) 
- cos(B,+f3) 


Q 

if 


si 

$,*- 


)-52 ] 
d: 


| 

cos(W-I,")~ Sy cos(¥-B ] 


d 
b 
d 
b 
d 


Alternate Solutions Table 2. The equations for Y and c? are the same for all com- 
binations of point reduction. 

Using the original precision points in the example, solutions 
sary to make a graphical solution before solving for the linkage | were made for the remaining five combinations of point reduc- 
mathematically. However, a graphical solution can be made ina tions. The results of these are shown in Table 1, and the linkages, 
few minutes and examined for its suitability. If the linkage has each in position 1, are shown in Fig. 2. 


From the foregoing example it can be seen that it is not neces- 


undesirable characteristics, a second solution can be made by 
. . . H 
considering the point reduction of another combination of any Conclusion 
two precision points. If desired, all six solutions of the linkage for 
the assumed precision points can be made graphically before 
selecting one to be solved analytically. The derivations of the 
five other combinations are similar to the one included in this 


By choosing different sets of precision points, there can be an 
infinite number of solutions for each function. In the method 
presented here, preliminary graphical investigations can be made 
rapidly to determine a suitable linkage before extensive calcula- 
paper. The final equations needed to solve problems are listed in tions are made. 
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DISCUSSION 


Ferdinand Freudenstein? and G. N. Sandor® 


The authors are to be commended for supplementing existing 
mathematical methods for function generation with a readily 
understood graphical method providing insight into mechanism 
geometry, proportions, and ranges. The discussers would look 


forward to the extension of the authors’ work to path generation 
sometime in the future. 


T. P. Goodman® 


The authors are to be commended for this contribution to 
kinematic synthesis which combines the graphical and analytical 
approaches. 

It should be noted that the six solutions shown in Fig. 2 are 
actually six of the infinite family of solutions which can be found 
by constructing Burmester’s circle-point and center-point curves.’ 
By choosing one solution from this infinite family, an additional 
condition such as a fifth precision point or a most favorable trans- 
mission angle® can be fulfilled. The authors’ six solutions give a 
survey of possible solutions to the problem from which the 
general appearance of the other solutions can be inferred by inter- 
polation and extrapolation. 


B. L. Harding’ 


The authors have performed a service to kinematics by con- 
verting Hain’s geometrical solution for four-phase, function- 
generator design into an analytical solution. Only one class of 
four-bar function generator is discussed, however, and it would be 
interesting and of value to know what would happen to the 
given equations where (a) the two cranks rotate in opposite 
directions, and (b) where the angle 6, as illustrated in Fig. 1 of 
the paper, is decreasing rather than increasing. The writer has 
no further criticism, but would like to elaborate on the subject. 

Hain” also describes a similar geometrical solution for five pre- 
cision points. The authors have charted the course for conversion 
of these geometric techniques to analytical solutions, but their 
method cannot be applied to Hain’s five-phase construction 
until someone can convert the following expression 


sin sin — 


2 2 


into an equation of the form 


Bos 
as = (6, Bis’ 


Lacking this expression, the technique under discussion is pres- 
ently limited to solution for four precision points, except in the 
special case discussed in the following: 

When the function to be generated is of the form f(z) = ax? + 
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Electric Company, Schenectady, N. Y. Assoc. Mem. ASME. 

7 R. A. Kraus and K. S. Rangasami, ‘Graphical Linkage Analysis,” 
Machine Design, vol. 25, 1953, pp. 177-182. 

8K. Hain, ‘“Drag-Link Mechanisms,’’ Machine Design, vol. 30, 
1958, pp. 104-113. 

® The Heald Machine Company, Worcester, Mass. 

1 Kurt Hain, ‘‘Angewandte Getriebelehre,”” Hermann Schroedel 
Verlag, Hannover, Germany, 1952, pp. 362-366. 
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bx + c, the authors’ method can be applied, and even simplified, 
to yield a solution for six precision points because of the following 
mathematical property of such functions: 


— In — 


— Lm = — > Lm > Ly > Ly 


INVERSION RELATIVE TO LINK Ley, 
IN PHASE * 1 


6 
2 
a= P + q Pr 
Fig. 3 
Design Procedure. Using notations similar to those of the 
paper, and referring to Fig. 3, herewith 


y = f(x) 
a = f(p) 


3 Chovse 21, 22, 23, 26 < a3 < 2 


4 Be = Bs = Bis 28x 
Bis = Bis Bus = Bis 2Bis 
is = = 


Hae) — fixe) — f(s) 


6 S, = 2a cos Bre Bu 


; = 2a cos S; = 2a cos — 
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+ 4 


Si(S;? — cos + S,(S,? — S;?) cos (cs + + — cos (cs Bs) 
7 y =tan- — 
Bas Bu 
S,(S3? sin + S.(S,? S3?) sin | Qe + + S;3(S2? S,?) Sin Q43 + 4 
rien S.2 — S8,? four-bar mechanism and which approximates the prescribed 


Bw\ 
+ d? — cos (v - 


Bernard W. Shaffer!! 


The authors have performed a valuable service by mathe- 
matically expressing Hain’s graphical solution for the synthesis of 
four-bar linkages. 

Use of the authors’ equations enables the designer to find a 
four-bar mechanism which satisfies the requirements of four 
precision points. Nevertheless, the designer has no assurance 
that the mechanism can generate a function containing these 
precision points. Even though the function actually generated 
contains the four precision points, very frequently it deviates 
considerably from the required function. The deviations are 
small, however, whenever the precision points also belong to a 
nonprescribed function which actually can be generated by a 


11 Professor™of Mechanical Engineering, New York University, 
New York, N. Y. Mem. ASME. 


function. 

The writer feels that before attempting to synthesize, one 
should investigate the existence of a suitable mechanism. An 
existence test is available for the four-bar mechanism, in the 
form of the compatibility equation developed by Shaffer and 
Cochin."? If after proper substitution the compatibility equa- 
tion is found to be satisfied identically over a finite range of varia- 
bles, one is assured of the existence of a four-bar mechanism 
which can generate the prescribed function exactly in that range; 
otherwise, only discrete points of the function can be reproduced. 

Obviously, the problem of generating a function is not the 
same as the problem of generating precision points. 


J. C. Wolford!* 


The authors have developed a simple analytical procedure for 
designing four-bar function generators for which the generated 
function fits the desired function at four precision points; that is, 


12 Bernard W. Shaffer and Ira Cochin, ‘‘Synthesis of the Four- 
Bar Mechanism When the Position of Two Members is Prescribed,” 
Trans. ASME, vol. 76, 1954, pp. 1137-1140. 

13 Associate Professor of Engineering Mechanics, University of 
Nebraska, Lincoln, Neb. Assoc. Mem. ASME. 
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four crank positions of the four-bar mechanism are co-ordinated. 
The development is based on a graphical method given by Hain 
which makes use of point-position-reduction, a technique appar- 
ently first used by Kraus.!4 Because of the simplicity of the 
authors’ analytical method and the fact that it can be used in 
conjunction with the graphical method as suggested by the 
authors, it should prove very useful to designers. 

The method does have some limitations, perhaps the principal 
one being the fact that only six solutions can be obtained out of 
the infinite number of possible solutions. Some solutions given 
by this method may result in mechanisms for which the actual 
function generated differs from the ideal function between pre- 
cision points by an amount which is too large. Other solutions 
may be unsatisfactory because of a poor transmission angle or 
inability of the linkage to operate over the entire required range 
of motion. However, if the function is one of the type which is 
suitable for approximate generation by a four-bar mechanism, 
usually one or more satisfactory solutions can be expected by the 
use of this method. 

The authors’ numerical example was checked using a graphical 
method by means of which all possible solutions can be obtained.!® 
The steps used in this graphical procedure for determining four- 
bar mechanisms with four crank positions co-ordinated are: 


1 Locate the relative poles Riz, Ris, Ris, Res, Res, and Rs, for 
the motion of link d relative to link b, Fig. 4, herewith. Three 
contrapole quadrilaterals can be drawn with these six relative 
poles as vertexes. Contrapoles are poles with no like ciphers 
in their indexes. For example, Ri; and Ra are contrapoles. The 
contrapoles are at the opposite vertexes of the quadrilaterals. 
In Fig. 4, RisR:sResF2s is one of the three contrapole quadrilaterals. 

2 Using one of the three quadrilaterals, the center-point curve 
for the motion of link d relative to link b is next constructed. 
The center-point curve is a focal curve and can be constructed as 

4 R. Kraus, ‘‘Lagenreduktion als Hilfsmittel zur Synthese ebener 
Kurbeltriebe,’’ Reuleaux-Mitteilung, vol. 3, 1935, p. 637. 


1s Rudolph Beyer, ‘‘Kinematische Getriebesynthese,’’ Springer- 
Verlag, Berlin, Germany, 1953, p. 100. 
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such, or it can be constructed making use of the property that 
points on the curve “‘see’’ opposite sides of the contrapole quadri- 
lateral under angles which are equal or supplementary. 

3 Next, pin joint D is chosen arbitrarily on the center-point 
curve. 

4 Pin joint C is determined making use of the fact that poles 
Ry; and Ry (or Riz) ‘see’? DA and CB under equal angles. 


Fig. 4 of this discussion shows this method applied in the deter- 
mination of the authors’ linkage 2-4. 

It is necessary that pin joint D lie on the center-point curve 
for the motion of link d relative to link b. All joints D of the 
authors’ linkages do lie on this curve with the exception of joint 
D of linkage 1-3. The authors’ location of this joint is shown in 
Fig. 4 and is labelled D,-3. The error can be found in the calcula- 
tion of angle y for this linkage. 

In the authors’ example, linkage 1-4 will not operate over the 
required range and hence is unsatisfactory. The function gener- 
ated by linkage 2-3 departs considerably from the ideal function 
between precision points and linkage 2-3 would probably be an 
unsatisfactory solution, see Fig. 5. Linkages 1-2, 2-4, and 3-4 
generate functions which closely approximate the ideal function 
and probably would be satisfactory solutions. 


Authors’ Closure 


The authors are pleased that six of the leading kinematicians 
in this country have shown an interest in this paper. Since their 
discussions are complete in themselves, we will make only a few 
comments. 

Dr. Wolford has shown that there is an error in the calculation 
for linkage 1-3, Table 1. The correct values are: Y = —86. 
6637 deg, b = 0.9592, and c = 1.9111. He has also plotted the 
curves for the ideal function and for linkage 2-3, Fig. 5, and has 
pointed out that this linkage is probably unsatisfactory. It 
should certainly be discarded in favor of one of the other three 
linkages mentioned. However, the curves can be used to illus- 
trate an important concept in function generator theory, i.e., 
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when a fairly satisfactory linkage is obtained, the structural error 
can usually be reduced by respacing the precision points. After a 
solution is obtained, the curves are examined and the precision 
points respaced to force the generated function closer to the ideal 
function. In Fig. 5, it can be seen that points chosen at 5, 30, 
70, and 95 degrees, would probably force the generated curve 
much closer to the ideal. This theory has been fully discussed 
by Freudenstein.!* 

The paper by Shaffer and Cochin, (Reference [11]) is a major 
contribution to the field of kinematic synthesis. Not the least 


16 F, Freudenstein, ‘Structural Error Analysis in Plane Kinematic 
Synthesis,” Trans. ASME, Series B, JouRNAL oF ENGINEERING 
FoR INDusTRY, vol. 81, 1959, pp. 15-22. 
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value of the paper is that it shows conclusively that some func- 
tions cannot be generated exactly. For these cases the designer 
should concentrate his efforts on reducing the structural error, not 
in trying to eliminate it. 

Mr. Harding has first posed an interesting problem, and then 
has made a worth-while contribution since the function that he 
treats is very common. All discussions, in one form or another, 
have indicated that there is still much work to be done in the field 
of synthesis. 

In conclusion, the authors should point out that signs are often 
troublesome when trigonometric functions are used. This paper 
is no exception. A rough graphical solution can be helpful. At 
the upper left of Fig. 1, letters d and b should be interchanged. 
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_ | A Theory for the Effect of Mean Stress 


on Fatigue of Metals Under Combined 


W. N. FINDLEY 


Professor of Engineering, 
Brown University, Providence, R. I. 
Mem. ASME 


Torsion and Axial Load or Bending 


The concept that alternating shear stress is the primary cause of fatigue with the normal 
stress on the critical shear plane as an influencing factor has been developed for the case 


of mean (or static) stresses superimposed on combinations of torsion and axial load or 
bending. The influence of the maximum stress of the cycle of stress on the allowable 
alternating stress for a given number of cycles and on the orientation of the critical shear 
plane is explored. The predictions of the theory are consistent with the known trends of 
fatigue data both for ductile metals and cast irons. The theory explains the fact that the 
influence of mean stress is weak for torsion and stronger for bending of ductile metals, 
but strong for both torsion and bending of cast irons. As far as is known this is the 
first rational theory for the influence of mean stress. 


I. EARLIER papers [1-3!! the effect of the normal 
stress acting on the plane of failure in fatigue was discussed. It 
was postulated that the alternating shear stress was the primary 
cause of fatigue, but the ability of the material to resist alternating 
stress was influenced by the normal stress acting on the critical 
shear plane. This theory was examined for fatigue under com- 
pletely reversed cycles of stress, and suitable relationships were 
derived. These were applied in particular to combinations of 
bending and torsion. 

In the present paper, this theory is extended to include the in- 
fluence of static or mean stresses superimposed on the alternating 
stress in a cycle of fatigue. In making the extension to include the 
effect of mean stress the same principles are employed; namely, 
that the alternating shearing stress is the prime factor causing the 
phenomenon of fatigue fracture but that the resistance to this 
action is influenced by the normal stress occurring on the critical 
shear plane. The alternating shearing stress is, of course, not in- 
fluenced by the mean stress but the normal stress occurring on the 
critical shear plane is affected by mean stresses. In other words, 
the normal stress is the sum of normal stresses resulting from 
both the alternating and mean stress. Under these conditions, 
the following relationships may be derived. 


Effect of Mean Stress 


In previous papers [1, 2], it was suggested as a result of ex- 
amination of available test data that the normal stress o¢ on the 
critical shear plane might have a linear influence on the allowable 
alternating shearing stress 7.¢, where 9 was defined as the angle 
which the principal stress o; made with the plane of the shearing 


‘stress 7.9. In another paper [3], a parabolic form of nonlinear in- 


fluence of the normal stress was investigated. This investigation 


showed little improvement in agreement between theory and test 
-data as compared to the linear theory. Accordingly, the present 


1 Numbers in brackets designate References at end of paper. 

Contributed by the Machine Design Division and presented at 
the Annual Meeting, November 30-December 5, 1958, New York, 
N. Y., of Toe American Society or MECHANICAL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 


-of the Society. Manuscript received at ASME Headquarters, June 30, 


1958. Paper No. 58—A-61. 
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theory will employ the linear form, as follows: The critical shear 
stress 7.9 decreases with increase in a9 from a value f when o¢.is 
zero according to the relation 


Te =f — kog (1) 


where f and k are constants for a given number of cycles to failure. 

In the following, this relationship is applied to the problem of 
a combination of alternating axial load (or bending) and alternat- 
ing torsion with superimposed mean (static) stresses having the 
same proportions of bending and torsion as the alternating 
stresses. In applying these principles when static stresses are in- 
volved, it must be remembered that 7.9 is the allowable alternating 
stress—not the maximum stress, and g9 is the maximum nor- 
mal stress, i.e., the sum of the normal stresses resulting from both 
the alternating and mean stresses. In the following, additional 
subscripts will be employed to differentiate between alternating 
stresses A and maximum stresses M. Thus Equation (1) be- 
comes 


Tea = f — koom. (2) 


The shear stress 79 and the normal stress o¢ acting on any plane 
in a biaxial stress field making an angle 6 with the greatest prin- 
cipal stress are given by 


01 — G3 


sin 20 (3) 


O1 + 01 — G3 
= 4 
06 + cos 26 (4) 


according to elastic theory, as may be seen from Mohr’s circle 
diagram. 
For combined bending ¢ and torsion 7 stresses, the principal 
stresses are 
a1, = (0/2) [(o/2)?+ | 


= 0. 


Substituting Equations (5) in Equations (3 and 4), the following 
relations are obtained 


Te = [(o/2)? + sin 20 (6) 
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= (a/2) + [(o/2)? + cos 20. (7) 


Equation (6) thus gives the shearing stress acting on the plane 
making the angle @ with the greatest principal stress, and Equa- 
tion (7) gives the value of the normal stress acting on the same 
plane. 

The difference Ar between the limiting value of the alternating 
shearing stress on any plane 7,94 and the applied alternating shear 
stress 7¢4 on the same plane is given by Equation (8). 


Ar = Téa — Toa. (8) 
Substituting Equation (2) in Equation (8) yields 
Ar =f — koom — Toa. (9) 


From Equation (6) 794 may be obtained by using the alternating 
components of bending o, and torsion tT, stresses. Similarly 
oom may be obtained from Equation (7) by introducing the maxi- 
mum bending oy and torsion 7 stresses. 

It will be understood for the present derivation that o, and 
T,4 are synchronous and in phase for combined bending and tor- 
sion. Hence oy and Ty occur simultaneously. Also, the pro- 
portion of mean bending stress to mean torsion stress will be 
taken the same as the proportion of alternating bending to al- 
ternating torsion stresses. Thus the maximum bending and 
torsion stresses also occur in the same proportion, i.e., Ty/o4y = 
T4/04 = a. The values of 794 and oy obtained as described 
may be substituted in Equation (9) with the following result: 


Ar =f — (koy/2) — k[(o4/2)? + cos 20 
— [(o4/2)? + (74)?]'/2 sin 20. (10) 


The plane for which Av is a minimum is the plane on which 
failure is most likely to occur. Differentiating Av with respect to 
6 yields 


= 2k sin 20 [(o4/2)? + 


— 2 cos 20 [(o4/2)? + (1) 


Setting the derivative equal to zero defines the most likely plane 
of failure as follows: 


[(o4/2)? + 


tan 20 = ; 12 
k[(oy/2)? + 

For zero mean stress oy = 04 and Ty = T, 80 that 
tan 20 = 1/k (13) 


which is the same result as obtained previously for zero mean 
stress [2]. 

Equation (12) shows that the plane of failure for nonzero mean 
stress is a function of the alternating and maximum stresses in a 
stress cycle as well as the influence of the normal stress k; whereas 
for zero mean stress the plane of failure is a function of k only, as 
shown by Equation (13). A change in character of the fracture 
with different mean stresses has been observed [4], but changes 
in angle have not so far been reported. 

When the mean stress is not zero the right-hand side of Equa- 
tion (12) reduces to 74/k7y, for pure torsion and o4/koy for pure 
bending. 

From trigonometry the values of sin 20 and cos 206 correspond- 
ing to Equation (12) are 


(04/2)? + 742]'/2 
[k2(o + + (04/2)? + 


k[(o¢/2)? + 
[k2(o 4/2)? + + (4/2)? + 


sin 20 = 


(14) 


cos 26 = 


(15) 
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Failure will occur at a given number of cycles only if the applied 
shear stress on the plane of failure is equal to the critical shear 
stress corresponding to the given number of cycles. Hence Ar 
must be zero if failure is to occur at the given number of cycles, 
Thus setting Av equal to zero in Equation (8), substituting Equa- 
tion (6) (for alternating stresses), and making use of Equation 
(14) for sin 20, the following expression is obtained: 


((o4/2)? + 
+ + (04/2)? + 


TOA (16) 


Equation (7) (for maximum stresses) may be expressed as fol- 
lows by the use of Equation (15): 


k ((oy/2)? + Ty?) 
+ + (64/2)? + 742]? 


Substituting Equations (16) and (17) into Equation (2) and re- 
arranging terms, the following relationship governing fatigue 
failure under combined bending and torsion with superimposed 
static stresses is obtained: 


f = {(o4/2)? + + + }'/2 + (koy/2). (18) 


Torsion. When Equation (18) is specialized to the case of pure 
torsion, for which oy, = a4 = 0, the following results: 


f= (ta? + (19) 


or solving for 74 and dividing by f, Equation (19) may be ex- 
pressed in dimensionless form as follows: 


= — (20) 


Since f and & are constants for a given material and number 
of cycles to failure, the fatigue strength 7, at the given number of 
cycles decreases with increasing magnitude of maximum stress ac- 
cording to the elliptic function, Equation (20). 

Equation (20) has been plotted in Fig. 1 for values of k of 0.1, 
0.2, 0.3, and 1.1. If k is a sufficiently small number, the second 
term of Equation (20) may be neglected so that the fatigue 
strength in torsion is nearly independent of the mean stress, as 
shown in Fig. 1. This result has been observed experimentally 
by many investigators and summarized by Smith [5]. Equation 
(20) does not discriminate between negative and positive maxi- 
mum stresses, so the greatest magnitude of either controls. 

The diagonal lines in Fig. 1 represent conditions for which the 
alternating stress and maximum stress are equal, i.e., for which 
the mean stress is zero. The region between the two diagonal 
lines cannot be tested in a torsion fatigue test because if the 
positive maximum stress in the cycle is decreased below the al- 
ternating stress, the negative (minimum) stress becomes larger in 
magnitude so that fatigue failure occurs as a result of the negative 


Fem = (64/2) + (17) 
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Fig. 1 Effect of maximum stress 717 on fatigue strength (alternating stress 
7A) in torsion 
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minimum stress, that is, outside the region between the diagonal 
lines. 

An upper bound for the value of 7,,/f may be found from 
Equation (20) by setting 7, = 0. This yields 


Tu/f = 1/k. (21) 


This upper bound for 7 may also be expressed as a proportion to 
40, the allowable alternating stress at zero mean stress. If Ty = 
Ir4o is substituted in Equation (21) together with Equation (19) 
for zero mean stress, the upper bound of the ratio 1 may be ex- 
pressed as follows: 


L = (1 + k*)'/2/k. (22) 


Limiting values of J, the ratio of maximum stress to alternating 
stress at zero mean stress for torsion are 10.1, 5.1, and 1.35 for 
values of k of 0.1, 0.2, and 1.1, respectively. Equation (20) may 
also be expressed in terms of the mean stress Tmean a8 follows, 
since Ty = T4 + Tmean: 


(74/f)? + kl(ra/f) + (Tmean/f)]? = 1. (23) 


Axial Load or Bending. For axial load or pure bending Equation 
(18) reduces to the following since 7), = T, = 0: 


f = [(o4/2)? + + (koy/2). (24) 


Solving for a4 and dividing by 2f a dimensionless expression re- 
sults 
o4/2f = [1 — (25) 


Thus for bending the fatigue strength o4 at a given number of 
cycles to failure varies with the maximum stress according to a 
parabolic function, Equation (25). In Fig. 2, Equation (25) has 
been plotted for the same values of k as used in Fig. 1. According 
to Equation (25) the fatigue strength o, decreases with tensile 
(positive) maximum stress and increases with compressive (nega- 
tive) maximum stress, as shown in Fig. 2. This prediction is in 
accord with many observations for axial load tests and for tensile 
maximum stress in bending. For symmetrica] cross-sections of 
test specimens in bending equal absolute values of mean stress in 
tension and compression are realized in every cycle. Since the 
greatest maximum stress corresponds to the tensile mean stress, 
and since the greatest maximum stress results in the Jowest fatigue 
strength, the tensile mean stress will control failure in bending. 

Again the diagonal lines in Fig. 2 represent the conditions when 
the alternating and maximum stresses are equal (mean stress is 
zero). In bending, however, the region between the diagonal 
lines is available. The minimum stress of the stress cycle will 
always correspond to a higher alternating stress for failure than 
the maximum stress of the cycle, according to Fig. 2. Thus the 
maximum stress controls for any value of maximum stress. 

The following upper bound for oy, is found when o, = 0: 


255 25 | 
\ | Pd 
2 
= 4 S 
05 \ 
5 
ONS | oN 
| | \ 
-| -05 05 +5 2 25 
/2t 


Fig.2 Effect of maximum stress o 1; on fatigue strength (alternating stress 
o4) in bending 
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Oy /2f = 1/2k. (26) 


As in the case of torsion just considered, the upper bound may 
be expressed by l’ where oy, = l’o 49 and 4p is the limiting value 
of the alternating stress at zero mean stress. The upper bound 
of the ratio l’ is 


[k + (1 + k?)'/2]/2k. (27) 


The limiting values of |’, the ratio of maximum stress to alternat- 
ing stress at zero mean stress for bending, are 5.5, 3.1, and 1.18 for 
values of k, of 0.1, 0.2, and 1.1, respectively. 

Equation (25) may be expressed in terms of the mean stress 
Omean 28 follows: 


(4/2)? + + (Gmean/f)] = 1. (28) 


Combined Torsion and Axial Load or Bending. The restrictions 
placed on the present solution include the conditions that 7, = 
ag, andTy = ady, [see text preceding Equation (10)]. Substitut-. 
ing these identities in Equation (18) and dividing by f the follow- 
ing dimensionless form of the relation is found when both bending 
and torsion are present: 


o4/2f = — (2koy/2f) — + (29) 


This relation is illustrated in Fig. 3 for a = 0.868 and k = 0.1, 
0.2, 0.3, and 1.1. 

It should be noted that the parameters in Equation (29) repre- 
sent only the bending component of the alternating and maximum 
stresses. The shapes of the diagrams are, however, proportional 
to the total stresses. The different scales used in Figs. 1 to 3 were 
selected to permit visual comparison, assuming that shearing 
stress was the co-ordinating parameter. 


| | 0.50 = 4207 | 
Li | 


-1.00 -O75 -O50 -G25 0 O25 G30 O75 100 125 
/ 


Fig.3 Effect of maximum stress 7 on fatigue strength (alternating stress 
a4) in combined bending and torsion. 74 = 0.868 04,7u = 0.8680. 


An interesting prediction stemming from Fig. 3 is that for ma- 
terials having sufficiently large values of k, such as k = 1.1, the 
effect of decreasing the maximum stress will be to increase the 
allowable alternating stress up to a point. A further decrease in 
the maximum stress of the cycle results in a decrease in the al- 
lowable alternating stress. The decrease is caused by the 
minimum stress of the cycle intersecting the downward curving 
left-hand portion of the curve of limiting stress. Thus the alter- 
nating stress corresponding to the minimum stress of the cycle 
would control for still smaller values of maximum stress. 


Comparison With Experiments 


Comparing Equations (20) and (25) (Figs. 1 and 2), it was ob- 
served that the maximum stress had a greater influence in bending 
than in torsion even for small values of k. For k = 0.1, for ex- 
ample, the influence of the maximum stress in torsion may be 
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considered negligible but not in bending. A similar effect has 
been observed experimentally as reported in reference [5] for a 
number of different metals. While individual] tests showed a 
variety of behavior, the following trends were observed: It was 
shown that the fatigue strength of ductile metals in torsion was 
nearly independent of the mean stress, whereas the fatigue 
strength under axial load generally decreased with increasing 
tensile mean stress. However, a re-examination of the evidence 
shows that the conclusion regarding torsion was influenced by 
the method of plotting so that a reappraisal suggests that in gen- 
era] the fatigue strength in torsion tends to decrease slightly with 
increasing mean stress. 

Reference [5] also shows that the fatigue strength of cast irons 
in tension decreases rapidly with increasing tensile mean stress, 
and the fatigue strength in torsion also decreases rapidly with 
torsion mean stress. Subsequent tests of individual metals also 
show a variety of behavior: A 65-ton alloy steel [6] showed a 
small decrease in fatigue strength of about the same amount for 
both alternating bending with superimposed static bending and 
alternating torsion with superimposed static torsion. The bend- 
ing curve corresponds approximately to k = 0.1 and the torsion 
curve approximately to k = 0.3. 

A 76S-T61 aluminum alloy [4] showed curves of fatigue 
strength at 3 X 104 cycles versus mean stress which approximate 
those predicted in Figs. 1 and 2 for k = 0.3. The corresponding 
value of the ratio of the bending fatigue strength at zero mean 
stress 6 to the fatigue strength in torsion at zero mean stress ¢, 
b/t = 1.62, yields a value of k of 0.24 according to Equation 2.56 
of reference [2]: 


b/t = 2/{1 + [k/(1 + k?)/]}. (30) 


The data of fatigue strength versus mean stress for this material 
at 108 cycles agreed best with the curve of k = 0.1 for bending 
and 0.3 for torsion, whereas the value of b/t remained the same 
as for 3 X 10‘ cycles. It is also interesting to note the similarity 
between the curves extrapolated from the data of reference [4] 
into the unavailable regions and the curves shown for these regions 
in Figs. 1 and 2. 

The results of tests of a 25S-T aluminum alloy [7] show an in- 
fluence of mean stress in both bending and torsion; however, the 
fatigue strength in torsion decreases more rapidly with increase of 
mean stress than does the fatigue strength in bending. 

Results from an SAE 4340 steel [8] show a slight tendency for 
the fatigue strength in bending to increase with tensile mean 
stress, and a slight tendency for the fatigue strength in torsion to 
decrease with mean stress. 

Tests of 75S-T6 aluminum alloy which were designed to dif- 
ferentiate between the influences of yielding and mean stress in 
tests at high mean stress [9] showed a reduction of 43 per cent in 
axial tension and 7 per cent in torsion resulting from the high 
mean stress alone. 

Results of tests of an alloy steel under axial load [10] and under 
torsion of hollow specimens [11] show a small decrease in fatigue 
strength with mean stress, but one which is somewhat greater for 
torsion than for tension. 

Thus while the general trend seems to be consistent with the 
proposed theory, results of individual investigations show con- 
siderable dispersion. These may be due to some of the following 
factors: Anisotropy, yielding, and hence change in structure of 
the material under the action of the maximum stresses, slip dis- 
tribution as discussed in reference [12], and various factors in- 
volved in surface condition of the specimens and testing tech- 
niques. 


Anisotropy 


Anisotropy may be considered in the analysis presented in this 
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paper in the same manner as in previous papers [2, 3, 4] by in- 
troducing a constant d before each shearing stress term in Equa- 
tions (18, 19, 20) etc., where d is of the nature of a stress concen- 
tration factor resulting from elongated inclusions or similar aniso- 
tropic discontinuities. The factor d should be a constant when 
the orientation of the bending stress o and twisting stress 7 re- 
main constant with respect to the texture of the material. It 
will be observed that the factor d will not change the shape of the 
curves of Figs. 1 and 2. So this will not account for variations in 
the influence of mean stress between tension (or bending) and 
torsion unless a different factor (d,,) is required for the maximum 
stress term than for the alternating stress term (d,4). It is, of 
course, entirely possible that d,, would not equal d,. 

Anisotropy would affect the value of k computed from the ob- 
served value of b/t since, as shown by Equation 2.58 of reference 
[2], the relationship involves d as a factor. This suggests that it 
might be better to determine the value of k from data of fatigue 
strength versus maximum stress as in Figs. 1,2. Then if all other 
factors but anisotropy were eliminated from these observations 
and dy were known to equal d, the anisotropy d could be deter- 
mined from the values of k and b/t substituted in Equation 2.58 of 
reference [2]. 


Orientation of Critical Shear Plane 


Equation (12) expresses the orientation 0 of the critical shear 
plane as a function of the alternating and maximum stresses of 
the stress cycle. Since Equation (18) shows that the alternating 
stresses are a function of the maximum stresses, Equation (12) 
may be expressed as a function of the maximum stresses only. 
Since the present analysis has been restricted to the condition 
74/04 = Ty/Cy = a, Equation (18) may be solved for rT, as a 
function of 7, anda. For the same reasons Equation (12) may 
be reduced to a function of 7, and 7,,. Then by substituting 7, 
from Equation (18) into Equation (12) the following equation 
for the orientation 6 may be found for combined bending and 


torsion: | 
_ 
tan 20 = + (31) 


The corresponding relation for pure torsion is 


tan 20 = [(f?/k®ryy2) — 1]'/2. (32) 
And the reJationship for axial load (or pure bending) is 
tan 20 = — 4f/koy)'/*. (33) 


Equations(31)—(33) have been plotted in Fig. 4, where the factor 
a has been taken equal to 0.868. These curves show that the 
orientation 6 of the critical shear plane varies rather widely with 
maximum stress as well as with state of combined stress. How- 
ever, for a given material, the range of possible values is much 
smaller than indicated by the entire diagram. 


| | 
Bending, Eq. \ 


Combined Bending 


Torsion, Eq. 30. 


30 
4 
| Torsion, Eg. 3/— 10 
-25 -20 -15 «+40 -05 1.0 5 
or k 


Fig. 4 Effect of the maximum stress of the stress cycle on the orientation 
of the critical shear plane 
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For tests in torsion at zero mean stress the smaller of the two 
possible values of @ are 45, 42.2, 40, 21, and 0 deg for values of k 
of 0, 0.1, 0.2, 1.1, and ©, respectively. These are the upper 
bounds for the smaller of the two possible values of @, since in 
torsion values of 7, cannot be less than 7 4, the alternating stress 
for zero mean stress. It is evident that the critical shear plane for 
zero mean stress is oriented a few degrees from the principal shear 
plane for small values of k and a few degrees from the principal 
plane for large values of k. An upper bound for all values of k was 
found as described earlier. For the upper bound the value of 0 
was zero for all values of k. However, the ductile metals (for 
which k is small) will have an upper limit of ry, set by yielding. 
For such materials a reasonable upper limit for 7,, may be taken 
as2740,forexample. The values of # corresponding to this value 
of Ty are 40 and 33 deg for values of k of 0.1 and 0.2, respectively. 
As noted before, when k = 1.1, 7, cannot exceed 1.35 740. 

At zero mean stress in axial load (or bending) the values of the 
smaller of the two possible values of # are 45, 41.8, 39, 20.5, and 0 
deg for values of k of 0, 0.1, 0.2, 1.1, and ©, respectively. An 
upper bound for all values of k was found as described in previous 
pages. For this upper bound the value of 6 was zero for all values 
of k the same as for torsion. Also, as for torsion, when ductile 
metals are considered (for which k is small) a practical upper 
bound will be set by yielding. Again choosing as a reasonable 
upper limit oy = 2049 the following values of 9, 38 and 29 deg, 
are found for values of k of 0.1 and 0.2, respectively. As noted 
before, 74, cannot exceed 1.18049 for k = 1.1. 

According to the theory 04 may take any lower values. How- 
ever, this is again limited by yielding and by buckling. It is in- 
teresting to note that when the maximum stress of the cycle is 
zero, that is the entire stress cycle is in compression, the value of 8 
is 45 deg for al] values of k. Thus for materials having a large 
value of k, such as cast iron, the value of @ for zero mean stress 
may approach zero (the plane of the principal stress); and as the 
maximum stress of the cycle is depressed toward zero the value of 
6 approaches 45 deg (the plane of the principal shearing stress). 

In this regard it is interesting to recall the tests of cast iron, 
aluminum alloy, and steel under a range of stress from zero to 
compression [13]. Figs. 7b and 12, together with the accompany- 
ing text of reference [13], show the fatigue fractures to have 
occurred on planes of 45 deg to the axis of the specimen for cast 
iron and aluminum. This is in agreement with the present theory. 
However, the planes shown in the reference are planes of crack 
propagation and not necessarily the planes on which fatigue 
cracks initiated. The plane of crack propagation for steel under 
compression cycles of loading, however, has been observed to be 
transverse to the axis of the specimen [9, 13]. Whether the 
cracks initiated on transverse planes has not been determined. 


This theory is admittedly incomplete at present. It has so far 
been related only to the special cases (although technically the 
most important cases) of axial load (or bending), torsion, and 
combinations 9: the two with superimposed mean stresses in the 
same ¢on;oviu.tion. The general state of stress will be considered 
later. The influence of the orthogonal normal stress has not been 
considered here nor has the fact that the state of stress is different 
at the two extremes of the stress cycle for many conditions so that 
a different critical plane would be predicted for each extreme of 
the cycle. It is planned to consider these matters at a later date. 
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DISCUSSION 
F. R. Shanley? 


It is interesting to present some of the author’s results in the 
form of dimensionless interaction curves. For example, his Equa- 
tion (28), for combined axial load and bending becomes 


Rmean (m*R 4? + 2km R,) (34) 
where 
Rimean 
Tmo 


Omo = base mean stress, i.e., the critical mean stress when 
= 0 
o4 


Ra 
Tao 


base alternating stress; i.e., critical alternating stress 
when Omean = 0 

m = (1 + k?)'2 — k where k represents influence of 

normal stress in fatigue 


The interaction curves determined by Equation (34) are shown 


2 Professor, University of California, Los Angeles, Calif.; also, 
Consultant, The RAND Corporation, Santa Monica, Calif. 
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in Fig. 5. The influence of k is seen to be relatively small. When 
k = 0 the interaction curves reduce to a single parabolic curve. 
Some of the test results available indicate that a value of k = 0.3 
gives a good fit when oo is taken as the ultimate tensile stress. 
(The test points shown in Fig. 5 are from a text of the writer.*) 


1.6 
k=0.3 
1.4 
SSS Ra (Alternating) 
1.2 
k=0 
| 1.0 
1.0 o! 
(Parabola) 
8 t 
0.1 
6 0.2 
Ra 0.3 
a 4 
Test 


R mean 


Fig. 5 Interaction curves for fatigue failure under combined alternating 
and mean stress 


For combined bending and torsion the author’s Equation (29) 
becomes 


R;? = 1 — (m*R,? + 2km R,) (35) 
where 
To 
R, 
Tro 


This equetion, plotted in Fig. 6, reduces to a circular are when 
k =0. 

A parabolic equation for the effects of mean stress was origi- 
nally proposed by Gerber and the circular equation for com- 


3 F. R. Shanley, ‘‘Strength of Materials,”” McGraw-Hill Book Com- 
pany, Inc., New York, N. Y., 1957. 
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Fig. 6 Interaction curves for fatigue failure under combined torsion and 
axial or bending stress 


bined bending and torsion was proposed by Gough, many years 
ago. By expressing the results of the author’s theory in inter- 
action form we can see clearly how the effect of normal stress will 
modify these earlier theories. 


F. B. Stulen* 


Tn reviewing this paper in which the author extends his earlier 
study (reference [2] of the paper) of pure alternating stresses to the 
case of combined alternating and steady stresses, the writer was 
gratified to see that his results are in exact agreement with some 
unpublished studies that the latter had made several years ago. 
The notation as well as the step-by-step procedures used in 
these earlier studies were somewhat different from those of the 
author. 

As the author repeatedly shows throughout his paper, the 
various fatigue characteristics of a material depend strongly on 
the parameter k which in turn is directly related to the ratio of 
the pure bending fatigue strength to the pure torsion fatigue 
strength by equation (30). For this reason, the ratio of these two 
fatigue strengths may be a fundamental material constant. A low 
ratio of these fatigue strengths (approaching unity) is indicative 
of a poor-quality material having low dynamic ductility, whereas 
a high ratio implies a large dynamic ductility. 


Author’s Closure 


The interesting contributions of the discussers to this paper 
are very much appreciated. 


4 Assistant Chief Engineer—Analysis, Propeller Division, Curtiss- 
Wright Corporation, Caldwell, N. J. 
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An Analysis of Critical Stresses and 
Mode of Failure of a Wire Rope 


Many wire-rope manufacturers and machine designers are under the impression that 
the significant stress in a wire rope 1s the tensile stress, or possibly the stress due to ten- 
sion and bending. This paper proves by mathematical analysis that by far the greatest 
stress in a wire rope results from Hertz contact stresses at points of contact of wire-on- 
wire, and asserts that the usua! mode of failure of a wire rope is fretting-fatigue initiated 
at such points of contact. Design relationships based on these concepts should be of 


great value to designers who use wire rope. 


hie DESIGN of a machine part is usually based on 
equations which relate forces, stresses, and dimensions. Each 
such design equation may be associated with a discrete critical 
point at which the state of stress is suspected to be more severe 
than in the surrounding region of the part. The efficiency of the 
design is largely dependent upon the proper selection of critical 
points, the judicious application of the principles of elasticity to 
calculate the state of stress induced at each critical point by the 
forces applied to the part, and upon the correct selection of the 
mode of failure and failure data. 

Most wire-rope design relationships of the present are based on 
critical-point selections and force-stress relationships which in- 
clude only the effects of a direct tensile stress, and a bending 
stress if the rope is bent around a sheave or drum. Drucker and 
Tachau [1],! by analysis of experimental data, demonstrated that 
the critical point in a wire rope is actually the point of contact be- 
tween wires and that the most important force-stress relationship 
associated with this critical point is the Hertzian contact stress 
induced at this point. The present paper presents theoretical 
analyses of states of stress at dominant critical points of a wire 
rope and the results of experimental tests which simulate the 
fretting conditions in a wire rope. These analyses substantiate 
the conclusions of Drucker and Tachau. 

To provide a concrete illustration of the principles of stress 
ansiysis set forth in this paper, a 6 X 7 wire rope, with a hemp 


1 Numbers in brackets designate References at end of paper. 

Contributed by the Machine Design Division and presented at the 
Annual Meeting, November 30-December 5, 1958, New York, N. Y., 
of THe AMERICAN SocieTY OF MECHANICAL ENGINEERS. 

Norte: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, July 1, 
1958. Paper No. 58—A-63. 


core, loaded in simple tension, has been selected as an example. 
Although the stress calculations presented herein apply to this 
specific rope, the principles illustrated have general applicability 
to other wire-rope constructions and orientations. 

In the analyses which follow, two critical regions in a 6 X 7 wire 
rope are identified and two separate force-stress relationships, 
based on two selected theories of failure for multiaxial static 
states of stress, are derived for each critical point. 


Assumptions 


The following assumptions apply throughout this paper: 


(1) All stresses in the hemp core are assumed to be zero. 

(2) The effects of internal friction between wires are neglected. 

(3) In each strand, consisting of six outer wires and one inner 
wire, the inner wire is assumed to be either larger or smaller in 
diameter than the outer wires, thus inducing contact either tan- 
gentially among the outer wires or between each outer wire and 
the inner wire. Both modes of contact may not exist simul- 
taneously. 

(4) In the vicinity of contact between wires, the wires are as- 
sumed to be cylindrical in shape. 

(5) It is assumed that all stresses are in the elastic region be- 
low the yield stress. 

(6) It is assumed that no strain in the rope is large enough to 
cause sufficient shape change to result in a distortion of the stress 
pattern. 


Failure Theories 


Wire ropes are usually made of cold-worked, heat-treated, high- 
carbon steel. Under static uniaxial stress conditions such ma- 


Nomenclature: 
a = semimajor axis of ellipse of con- B = term dependent on geometry of B = lay angle of strands in rope, degrees 
tact, in. contacting bodies A term dependent on material and 
b = semiminor axis of ellipse of con- E = modulus of elasticity, psi geometry of contacting bodies 
tact, in. F normal force between contacting o, = normal stress in axial direction in 
d = wire diameter, in. bodies, lb/in. core wire of strand, psi 
r = cylinder radius, in. F, radial force between wires at a, normal stress in axial direction in 
x = axial co-ordinate, in. parallel-wire critical region, outer wires of strand, psi 
y = tangential co-ordinate, in. Ib/in. Ont normal stress due to direct tensile 
z = co-ordinate radially inward from - = radial force between wires at force only, psi 
cylindrical surface, in. crossed-wire critical region, lb Cer critical normal stress, psi 
A = term dependent on geometry of T tensile force per strand, lb Te critical shear stress, psi 
contacting bodies a lay angle of wires in strand, degrees mn Poisson’s ratio 
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terials are on the border between the ductile and the brittle 
state. Consequently, two failure theories for multiaxial states of 
static stress are used in this analysis: (1) The maximum-normal- 
stress theory, which is applicable for materials in the brittle state; 
and (2) the maximum-shear-stress theory, which is applicable for 
the ductile state. The actual behavior of most wire-rope ma- 
terials will probably fall somewhere between the predictions of 
these two theories. 


Critical Regions 


If it be asserted that the contact stresses between wires are the 
most significant stresses, there are two important critical regions 
in a 6 X 7 wire rope subjected to a tensile force [2]. The first of 
these, called the parallel-wire critical region, is located at the 
point of contact between the wires of a strand where the wires are 
in contact along a line, as two cylinders with parallel axes. 
Whether this critical region is between an outer wire and the core 
wire of a strand or between two outer wires is dependent on the 
diameter of the core wire, as indicated in Figs. 1(a, 6). The 
critical stresses are higher in the arrangement of Fig. 1(a) because 
the direct tensile stress is higher in the core wire than in the outer 
wires [2]. Therefore, the arrangement of Fig. 1(a) is the more 
serious of the two. 


& 
(a) (b) 


Fig. 1 Two possible conditions of contact of wires within a strand of 
wire rope 


Fig. 2 Two possible conditions of contact of adjacent strands of a wire 
rope 


The second critical region, called the crossed-wire critical re- 
gion, is located between wires of adjacent strands. Here point 
contact exists, since the wires simulate two cylinders in contact 
with axes askew. As indicated in Figs. 2(a, 6), two modes of con- 
tact may exist depending on the axial location. The arrangement 
depicted in Fig. 2(a) is the more serious of the two. 


Forces 


A tensile force on a wire wrapped in the form of a helix around 
a cylinder results in a radial force between the wire and the cylin- 
der of [2,3] 


F = (T/r) sin? a, (b/in.) (1) 


By using this relationship and the geometry of the 6 X 7 wire 
rope, the radial contact forces at the two critical regions can be 
written. 

At the parallel-wire critical region [2], Fig. 1(a), 


T cos? a sin? a@ 


d(1 + 6 cos? a) 


(b/in.) (2) 
At the crossed-wire critical region, Fig. 2(a), 
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_ aT sin? B 
tena’? (3) 


If the wire lay-angle a be assumed equal to the strand lay- 
angle 8, Equation (3) can be simplified. The lay-angles of many 
wire ropes were examined and this assumption proved reasonable. 
The largest difference between the lay-angle of the strands and the 
lay-angle of the wires was found to be one degree in the ropes 
studied. The simplified expression for F, then becomes 


F, = 0.1745 T sin a cos a (4) 


The two lay-angles will be assumed to be equal throughout the 
rest of this paper. 


Direct Tensile Stress 


The tensile stresses induced in the wires of a single strand of 
wire rope by a tensile force on the strand are given in Equations 
(5) and (6) [2]. These equations are based on the assumption 
that all wires undergo the same strain in the direction axial to the 
strand. 

In the core wire, axial to the wire, 


4T 


= md?(1 + 6 cos? a) 6) 
In the outer wires, axial to these wires, 
4T cos? a 

(6) 


eT md?(1 + 6 cos? @) 


Contact Stresses at the Parallel- Wire Critical Region 


The general equations for contact stresses between two parallel 
cylinders of equal diameter pressed together radially are con- 
veniently presented by Seely and Smith [4] as, 


/2 


where the quantity 6 is one-half the width of the rectangular 
area of contact between the two wires and may be found from 


the equation, 
(74) (10) 


The term A is a function of the geometry of the contact region 
and the wire material. For the parallel-wire critical region, A is 
given by, 


E 


(11) 


Contact Stresses at the Crossed-Wire Critical Region 
The general contact-stress formulas apply at the crossed-wire 


critical region. They can be simplified by recalling that the con- 
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gion 
A is 


(11) 


con- 


tacting wires are cylinders of equal radii, are made of the same 
material, and cross at twice the lay angle. Applying these facts, 
the expression for A can be simplified to 


(1 — 
sin? a@ 
) 


If the material of the wire be restricted to steel, with u = 0.26 
and E = 3 X 10’ psi, the variation in A with a is very small be- 
tween the practical limits, a = 12 toa@ = 20deg. With these re- 
strictions, A is closely approximated by 


A = 3.03 X 10-*d (13) 


A= (12) 


Critical Stresses at the Parallel-Wire Critical Region 


The critical normal stress at the parallel-wire critical region can 
be found by substituting Equations (2), (10), and (11) into 
Equation (9) and simplifying. This gives 

[ 2TE cos? a sin? a 
md?(1 + 6 a)(1 — 


(14) 


The critical shear stress at the parallel-wire critical region 
occurs beneath the surface of the wire [2] at z/b = 0.407. It is 


b b \? 
= 
Ter = 0.288 — + 4.95 X 10 (. >) (15) 


where 


2TE cos? a sin? a Va 
(16) 


b 
A E.. + 6cos*a)(1 — 


Critical Stresses at the Crossed-Wire Critical Region 


The critical normal stress at the crossed-wire critical region 
can be found from an equation given by Seely and Smith [4] 


Sor = —Co(b/A) (17) 


where 
b = (FA) (18) 


and A is given by Equation (12). 

The factors cg and c, are functions of the lay anglea. They are 
presented graphically on pages 356 and 357 of Seely and Smith 
[4]. To evaluate co and ¢, for any particular lay angle, a quan- 
tity B/A must first be determined from the equation 


( sin? a sin? a@ (19) 
1+ 3 ) - ( - 3 ) cos 2a 


Thus the critical normal stress may be calculated by combining 
Equations (12), (19), (18), and the graphs of Seely and Smith. 

The shear stress at the crossed-wire critical region is a complex 
function of the normal contact stresses and the direct tensile 
stress. The shear stress reaches a maximum, or critical value, at 
some distance below the wire surface. An analytical determina- 
tion of this critical stress is extremely difficult because the per- 
tinent variables are not explicit, but are contained implicitly 
within elliptic integrals. In addition, the presence of direct ten- 
sile stress hampers identification of the particular shear stress 
which becomes critical. 

The complexity of an analytical evaluation indicated that a 
graphical investigation of the critical shear stress at the crossed- 
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Fig.3 Plot of critical shear stress versus direct tensile stress at the crossed- 
wire critical region of a 6 X 7 wire rope under axial tension 
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Fig. 4 Semilog plot of critical shear stress versus direct tensile stress at 
the crossed-wire critical region of a 6 X 7 wire rope in axial tension 


wire critical region should be attempted. This was successfully 
done [6] and the results are shown in Figs. 3 and 4. The log-log 
plot, Fig. 4, is included because it enlarges the region where 
stresses are reasonable. It is also used later for curve-fitting. 
The critical shear stress is plotted against the direct tensile stress 
in the outer wires as calculated by Equation (6). The data for 
these curves were found using material presented by Thomas and 
Hoersch [5], together with independent calculations [6] involving 
evaluation of the elliptic integrals. These graphs are valid only 
for steel wire rope. 

The particular principal shear stress which is critical at the 
crossed-wire critical region changes at high stress levels due to the 
effect of the direct tensile stress being added to the contact stress. 
The point of change is indicated on the graphs of Figs. 3 and 4 by 
dotted lines. 


Simplified Equations 


If the 6 X 7 wire rope were made of steel, with a modulus of 
elasticity of 3 X 10’ psi, and a Poisson’s ratio of 0.26, simple ex- 
pressions can be written relating the critical stresses and the 
direct tensile stress. These expressions, shown in Table 1, were 
derived analytically for all of the critical stresses except the 
critical shear stress at the crossed-wire critical region. For the 
latter case, the equations were constructed by curve-fitting the 
log-log plot of critical-shear-stress versus direct-tensile-stress, 
Fig. 4. 


Magnitude of Contact Stresses 


The critical stresses found are the highest stresses in the 
wire rope. The effects of both the contact stresses and the direct 
tensile stress were considered in calculating these combined 
stresses. The relative importance of the critical stresses can be 
determined by comparing them with the direct tensile stress 
This comparison is presented graphically in Figs. 5 through 8. 
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Fig. 5 Plot of ratio of critical normal stress/direct tensile stress versus 
direct tensile stress at the parallel-wire critical region for a 6 X 7 steel 
wire rope 
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Fig. 7 Plot of ratio of critical normal stress/direct tensile stress versus 
direct tensile stress at the crossed-wire critical region for a 6 X 7 wire 
rope 
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Fig. 6 Plot of ratio of critical shear stress/(direct tensile stress/2) versus 
direct tensile stress/2 at the parallel-wire critical region for a 6 X 7 wire 
rope 
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Fig. 8 Plot of ratio of critical shear stress/(direct tensile stress/2) versus 
direct tensile stress/2 at the crossed-wire critical region for a 6 X 7 wire 


rope 


Table 1 Relationships between critical stresses and direct tensile stress in a 6 X 7 steel wire rope under simple tensile loading 


- Helix angles of wires in strands and strands in rope, deg — 


Critical stress-——-—— a = 12 a= 15 a=17 a = 20 
Critical 
normal Cor = —815(o21)'/2 Car = —1000(o2:)'/2 Cor = —1120(o21)'/2 Cor = —1290(o21)'/2 
stress 
Parallel-wire critical 
region 
Critical 
shear Ter = 235(ozt)!/2 Ter = 291(o2t)'/2 Ter = 324(ozt)'/2 Ter = 327(ozt)'/2 
stress + 0.562: + 0.502: + 0.502 
Critical 
normal ger = oer = oer = ocr = 
stress 
Crossed-wire critical 
region 
Critical 
shear Ter =8270(o21)'/2 Ter = 9560(o24)'/8 Ter = ter = 11,500(o21)'/2 
stress 


In these graphs, the ratio of the critical-normal-stress to the 
direct-tensile-stress is plotted against the direct-tensile-stress for 
the maximum-normal-stress theory of failure. The ratio of the 
critical-shear-stress to one half the direct-tensile-stress is plotted 
versus one half the direct-tensile-stress for the maximum-shear- 
stress theory of failure. 

This analysis shows that the critical stresses in a wire rope are 
due largely to contact stresses and are much higher than the 
direct-tensile-stress, the stress which is often used as the primary 
criterion of failure for a wire rope. This study predicts, for ex- 
ample, that if a 1.5-in. diameter, 6 X 7 wire rope were subjected 
to an axial tensile force of 70 pounds, the stresses at the critical 
points identified in this analysis should induce local yielding. 
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Fretting in Wire Rope 


The equations derived in the previous sections of this paper in- 
dicate that very high contact stresses may develop in regions of 
point or line contact within a wire rope. An additional phe- 
nomenon, called fretting, probably also occurs at such regions and 
may contribute significantly to the initiation of fatigue cracks. 

Experiments were performed to investigate this fretting phe- 
nomenon. Fretting, under point-contact conditions similar to 
those found in a wire rope, was induced in steel specimens [7]. 
Three different degrees-of-severity of fretting were used in the 
test program. Prot failure tests were applied subsequent to the 
fretting treatments to determine whether or not fretting had re- 
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Table 2 Summary of results of fretting tests 


Mean prot Unbiased 
failure standard 


me Sample stress, deviation, Low range 

Test condition size psi psi value 
Nonfretted, pol- 

ished, steel 15 78,200 5456 70,000 
Mildly fretted, 

polished, steel 15 77, 280 4155 70,000 
Medium fretted, 

polished, steel 15 71,850 5492 60,000 
Severely fretted, 

polished, steel 15 67,700 6532 51,000 


duced the endurance limits. The results of these tests on steel 
specimens are given in Table 2. 

It is evident that fretting did cause fatigue damage. Subse- 
quent microscopic observations revealed that fretting had in fact 
initiated fatigue cracks. 


Probable Mode of Failure of a Wire Rope 


The probable mode of failure of a wire rope may now be pro- 
posed. High contact stresses, especially at regions where wires 
cross to give point-contact, cause stresses far greater than the 
yield stress. Hence yielding and flattening occurs, and subse- 
quent stresses are of yield-stress value. Relative motion of wire- 
on-wire occurs due to strains in the wires, particularly strains due 
to rope bending. These relative motions and contact pressures 
cause fretting, which induces fatigue cracks. The fatigue cracks 
first propagate through the highly stressed regions near the points 
of contact. By the time the cracks propagate out of the highly 
stressed contact regions, they are sufficiently large that they will 
continue to propagate under the less severe tensile and bending 
stresses which prevail elsewhere throughout the wires. 

Thus fretting-induced fatigue cracks propagate to cause com- 
plete failure of individual wires. After several wires have frac- 
tured, the rope is so weakened that it is dangerous and must be 
replaced. 
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W. J. Kaufman? 


The authors are to be commended for neatly packaging the in- 
formation pertinent to calculation of stresses in wire ropes, and 
for pointing out the importance of the effects of Hertzian contact 
stresses in producing wire-rope failures. In naval arresting-gear 
applications, where wire rope is used as a machine element which 


_ is of primary importance to operational reliability, it has become 


quite evident to us over the years that wire-rope manufacture is 
an art, not a science. Any scientific effort which enhances the 
art and aids in transforming it into a science accordingly fills a 
sore need, and this paper does exactly that. 

The authors show that calculated crossed-wire stress magni- 
tudes are intolerably high under their assumption that “all 
stresses in the hemp core are assumed to be zero.’”’ It should be 
emphasized that in any practical wire rope this assumption does 
not hold. In fact, the wire-rope designer tries to choose a hemp 
core of such size and physical characteristics that it acts to sup- 
port the strands fully and preclude completely the existence of 
any crossed-wire stresses. In the practical case, he rarely suc- 
ceeds completely, so the rope fails at a slightly lower load as a 
result. It follows, then, that control of the physical characteristics 
(the size, compressibility, frangibility, and like properties) of the 
hemp center core is one of the most important of the many wire- 
rope design factors, and is accordingly a fertile field for detailed 
scientific study. 

The writer’s colleague, Dr. Ringleb, has shown that the parallel- 
wire contact stress effects, calculated by the authors’ method, 
may well explain most of what the wire-rope designer now dis- 
misses as “rope efficiency,’ even though this calculation takes 
account of only the inner forces within a strand, arising from the 
helical construction of an individual strand. 

If now, we consider that in the practical case the hemp center 
is actually a kind of elastic foundation for each strand, and as 
such is capable of applying a radial compression force directly to 
half the outer wires of each strand, then it is apparent that the 
hemp center acts to increase the parallel-wire contact stresses in 
the strand. The radial loads applied to the strand are a function 
of rope lay and axial tension, and are superposed on the authors’ 
stresses derived from strand lay and axial tension. The effect 
must be to raise the critical’ parallel-wire stress level to a value 
more nearly comparable to the one derived from tests of individual 
wires in the simplest stress state. This probably accounts for the 
remainder of the difference discussed by Dr. Ringleb. 

One other observation seems pertinent. The mode of failure 
of the individual wires in a wire rope was noted by the authors as 
inherently borderline. In the Navy’s arresting-gear applications, 
where over-all factors of safety for wire rope are of the order 1.5, 
and where fatigue is not a factor because of short service life 
from other considerations, the mode of wire-rope failure is always 
noted to be typical of ductile rather than brittle material—a 
point that might not be the case in other applications. 

The writer again commends the authors. There is serious need 
for a scientific approach to the wire-rope-design problem, of the 
sort exemplified by the present paper. Much remains to be done. 


F. 0. Ringleb* 
Among all mechanical machinery using wire ropes for the 


2Chief Project Engineer, Naval Air Engineering Facility (SI), 
Naval Air Material Center, Philadelphia, Pa. 
3 Woodbury Heights, N. J. 
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transmission of forces the aircraft arresting gear requires the 
highest performance of the cable by which the landing airplane 
has to be decelerated. At present landing speeds, these per- 
formance requirements are so high that even the best available 
wire ropes begin to fail. Therefore, any investigation with the 
purpose of explaining the stress mechanism in a wire rope, and to 
show ways of better cable construction, is of the greatest value. 
The present paper is a valuable contribution to the stress 
analysis of a wire rope and when the writer offers in the following 
some critical remarks, it is with the intention only of suggesting to 
the authors a continuation of their work in a specific direction. 
The authors point out that the highest stresses in a wire rope 
with hemp core are Hertz contact stresses between the wires. The 
contact stresses between parallel wires can be computed from 
simple formulas which have been given in the paper. Depending 
on the lay angle these stresses are a few times higher than the 
direct tensile stresses. In the case of crossed wires, however, the 
computation is complicated and the resulting stresses are tre- 
mendously high. In the writer’s opinion, this case is of no prac- 
tical importance. It can occur only if the strands are in contact. 
But in a good cable such contact is avoided by using a hemp core 
which is somewhat larger than the inner space between the 
strands, so that a small gap separates them. In this way the 
hemp core can act as an elastic support for the wire strands. The 
size of the gap has to be chosen so that they just close at maximum 
working load of the cable. The cross-wire case cited in the paper 
appears as a consequence of the assumption that all stresses in 
the hemp core are zero, which is not true in practice. The writer 
suggests basing the stress analysis of the wire rope on the contact 
stresses of parallel wires only; h r idering the hemp core 
as an elastic support for the strands. 
The writer believes that from the contact-stress theory of 
parallel wires as presented in the paper an interesting conclusion 
can be drawn. He has found for a6 X 7 cable, using the authors’ 
graphs, that the average stress in a wire produced by the parallel 
contact stresses is equal approximately to the breaking stress of 
the steel if the cable is loaded with its breaking tension. Thus 
from the paper can be explained quantitatively why a cable 
breaks at a considerably lower stress than the breaking stress of 
the wire steel. Actually the breaking stress of any cable should 
be computable in this way if its geometry and the elastic data of 
its material are known. The writer repeats that this paper is a 
valuable contribution to the stress analysis of a wire rope. 


C. M. Zerr* 


Wire rope has long been considered a piece of machinery by 
most producers and many users. Those who fail to recognize it as 
such, fail to design or use it properly. Wire rope is made up of 
many moving parts. These parts include the core which is made 
up of either a specially constructed hard fiber or an independent 
wire rope core, commonly known as an IWRC, and the various 
steel wires which will number from 42 in the 6 X 7 construction 
to 294 in a6 X 49. The wire sizes are drawn to very close toler- 
ances, usually to less than 0.001 in. and with anywhere from 2 to 6 
different sizes in the strands. 

As with any piece of machinery, the chemistry of the material is 
important. Carbon steel in the 1000 series is the most popular, 
because of its workability and economics. Special steels are used 
upon occasion. The amount of cold-working which is reflected in 
tensile strength and torsions of the wire is determined by the end 
use. No heat-treating has ever been done successfully on finished 
drawn rope wire as a means of improving the wire characteristics. 

In addition to the closely controlled physical properties of the 
component parts, great care must be exercised in the design and 
actual assembly of these parts. The lays or length of helix of the 


4 Chief Engineer, Union Wire Rope Corporation, Kansas City, Mo. 
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wires in the strand and the strand in the rope are closely related to 
each other and also to the end use. 

Finally, there is that item known as “the art” of closing which 
cannot be measured and comes only through experience. 

In the design of a wire rope, we must take into consideration 
the end use. This will determine the specifications which must 
be set. up for the core wires, lays, and constructions. Tensile and 
bending stresses are naturally of primary consideration but, as 
pointed out by the authors, critical stresses are set up within the 
rope which must be considered. 

In their initial approach to the problem, the authors indicate 
that stresses are higher in a strand with a larger center wire and 
with the 6 cover wires supported by this center wire than the con- 
dition where the 6 outer wires are in contact with each other. In 
actual] practice, it has been found necessary to support the cover 
wires on the center wire properly; this whether in a standing or 
operating rope. It must be remembered that friction is always 
present and its effect contributes materially to the physical prop- 
erties of the strand. 

Wire rope, being a machine, has moving, working, and stressed 
members. Wear takes place between these working parts as in all 
machinery. The bending of a rope over a sheave causes the wires 
making up the rope to slide back and forth against each other. 
A similar motion takes place between the strands. Continued 
flexing of the rope over the sheaves and repeated loadings wear 
the core and wires along the points and lines of contact until the 
rope ultimately will not perform safely as designed. 


Fig.9 Center wire showing fluted condition 


Fig. 9 shows a center wire which has been in actual service. 
The wire, which was originally round, is now fluted. This con- 
dition was induced by the movement of the cover wires along the 
surface of this center wire. The friction developed between wires 
in a strand, which has been the cause of the wear, fluting, or 
galling, is great enough that, within 2-3 ft of a broken end, the 
wire will pick up its proportionate share of the total load. 

Another example of contact stresses at the parallel-wire critical 
region is found between the outer wires of the strand and the 
IWRC. Fig. 10 shows the typical surface appearance and Fig. 11 
the cross section. 

Contact stresses at the crossed-wire critical region, which is 
commonly referred to as strand nicking, are illustrated in Figs. 12 
and 13 which show the surface appearance and cross-sectional 
view. In both Figs. 11 and 13 can be observed the effect of the 
parallel-wire (fluting) wear as observed in Fig. 9. They are 
marked as points A. 

There is one other area which has to be given serious considera- 
tion in properly designing wire rope. This area is that portion of 
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Fig. 10 Surface appearance resulting from contact stresses between 
center wires and IWRC 


Fig. 11 Cross section of wire in Fig. 10 


Fig. 12 Result of contact stresses at crossed-wire critical region, com- 
monly called strand nicking ; 


Fig. 13 Cross-sectional view of wire in Fig. 12 


the rope which comes in contact with the external surroundings. 
Figs. 14 and 15 shows a typical wire failure across the crown of a 
strand. Note the actual deformation of the wire as it has flowed 
around the supporting wires inside the strand. The flowing of 
the metal is the result of pressure between a sheave and the wire. 
In this condition there is generally no fretting of the adjacent 
Wires as the unit pressures have been reduced to the point that the 
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Fig. 14 Typical wire failure across crown of strand 


Fig. 15 Cross section of wire in Fig. 14 


Fig. 16 Martensitic structure at surface of a wire; crack acting as a 
notch, ultimately causing fracture : 


lubricating film of oil does not break down. In addition, surface 
abrasion between the wire and the surrounding elements takes a 
toll on rope life. A fifth condition which often is the cause of early 
retirement is the result of excessive bending. Under these con- 
ditions, all the wires in each strand will break along the crown of 
the strands. 

Friction, which the authors have neglected in their assump- 
tions (see item 2), is an ever-present culprit in reducing the ef- 
ficiency and useful life of wire rope. When two wires rub to- 
gether or the rope slides over some object, the momentary heating 
due to the friction and sudden quench changes the iron-carbide 
formation to martensite. Fig. 16 is a typical example of the ex- 
ceptionally hard and brittle martensite along the surface of a wire 
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with a typical crack acting as a notch which ultimately will cause 
a wire fracture. Owing to the hard enamel-like characteristics 
of this formation, it easily chips off before the sample reaches 
the laboratory and thus is overlooked as a possible source of the 
wire-rope failure. 

Fig. 17 shows two strands and an IWRC which illustrate the 
parallel and crossed-wire condition as they appear in a used rope 
The broken wires in the IWRC are the result of the compressive 
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Fig. 17 Two strands and an IWRC illustrating parallel and crossed-wire conditions 
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Fig. 18(b) 


and restrictive forces imposed on it by the main strands of the 
rope. 

The value of 30,000,000 psi, which has been assumed for the 
modulus of elasticity for a drawn rope wire, is of questionable 
value. Fig. 18 is a graph of the stress-strain curve for an 0.060 and 
0.120-in. wire which is commonly used in the average wire rope. 
It will be noted that the wire forms a continuous curve. As the 
wire or the rope does not follow the usual known curve from which 
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is determined the value of Z, it thus becomes impractical, if not 
impossible, to establish a valid value. In addition, repeated load- 
ings as illustrated in Fig. 19 for these same two wires throws ad- 
ditional light on the constant change taking place as the rope is 
used. 

The diagram in Fig. 20 shows a typical arrangement of the 
wires and strands in a6 X 7 wirerope. The two illustrations used 
by the authors are limiting conditions rarely found in actual prac- 
tice. The strands and rope are designed to the circumscribed 
circle. 

The various illustrations used in this discussion were taken from 
a wire rope which had given excellent results in the field. The 
fact that there was parallel-wire wear and crossed-wire nicking 
present and yet they had not resulted in wire failure at these 
points seems to indicate that the manufacturers and designers 
have learned to live with a condition which is not ideal but which 
can be reduced in its over-all effect. There are times, however, 
when conditions are such that these critical areas become the 
direct cause of the rope failure and in such cases, a re-evaluation 
of the specifications becomes necessary. 

The mode of failure of a wire rope at points of “high contact 
stresses,’ as suggested by the authors, while valid as to an analyti- 
cal analysis, must be qualified in practical application. Factors 
contributing to the rope design include: 

1 Physical properties of the wire. 

2 Strand construction. 

3 Rope fabrication. 


4 Inner strand and rope stresses from parallel and cross-wire 
action. 
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Fig. 20 Typical arrangement of wires and strands in a 6 X 7 wire rope 
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5 Type of iupricants used. 
6 Conditions to which the finished product will be subjected. 


These typical factors are all variables within themselves and in 
relation to each other. Their individual and collective influence 
on the ultimate rope service is therefore very difficult to ascertain 
mathematically. 
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The study of critical-stress areas in a wire rope presented in the 
paper are extremely interesting, informative, and will no doubt 
become a part of the general science of machine design. Other 
critical factors that have been pointed out in this discussion indi- 
cate that the authors have just started in their task of analyzing 
the compounded stresses that are inherent within the static or 
dynamically loaded machine we call a wire rope. 
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An Application of Statistics to 


M. F. SPOTTS 


Professor of Mechanical 
Engineering, Northwestern 
University, Evanston, Ill. 
Mem. ASME 


the Dimensioning of Machine Parts 


An outstanding industrial achievement during recent years has been the development of 
quality control and acceptance sampling methods for the purpose of maintaining the 


high standards of manufactured products. It is equally important to make use of 
statistics in the engineering office in connection with the dimensioning of machine 
parts. This paper will be concerned with the dimensioning of assemblies, and will 
demonstrate how tt is possible to maintain tolerances for proper engineering function- 
ang, and at the same time widen the tolerances so far as production is concerned. 


A PRODUCTION process is said to be in statistical con- 
trol when all assignable causes of error are located and removed, 
and the variation in dimensions is due solely to chance or random 
causes. These are inherent in the process and are present even 
though all conditions are held as constant as humanly possible. 

When parts are made in large quantities, and the production 
system is in statistical control, small errors occur more frequently 
than large errors. If a plot is made of the frequency of occurrence 
of an error versus its magnitude, the bell-shaped or normal curve 
of Fig. 1 results. This is the probability curve or Gaussian law of 
error, and the derivation of the equation can be found in books on 
statistics or advanced calculus. The curve is symmetrical about 
the mean or arithmetic average of the measured dimension. The 
area under the curve represents the number of piece parts. The 
standard error, usually called sigma, a, is equal to the root mean 
square of the errors of the parts. The errors for the dimension 
are squared and added together, the sum is divided by the number 
of parts, and the square root taken. 

Dimension % + wu is shown in Fig. 1(c). Here a is the mean or 
average dimension, and +uw is the tolerance. Unilateral toler- 
ances have many advantages for shop drawings. However, when 


Fre 


Error 
“a tu 
olerance 
u 
<< — 
Fig. 1 
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Note: Statements and opinions advanced in papers are to be 
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statistical aspects are to be considered, bilateral tolerances are 
usually more convenient. In Fig. 1(a), the production process is 
such that three standard deviations are equal to the tolerance w. 
This is the so-called ‘‘natural tolerance” for a controlled produc- 
tion process. The area under the curve then represents 99.73 per 
cent of the entire production. 

In Fig. 1, and throughout the paper, it is assumed that the 
mean of the production process coincides with the mean value of 
the dimension. 

In Fig. 1(b), the errors of production are larger. The area 
under the curve to the left of a — u represents the proportion of 
parts with a value for the dimension less than 7 — u. Likewise 
the area under the curve to the right of 7 + u represents the pro- 
portion of parts larger than @ + wu. ‘ 

Thus the shape of the probability curve can vary widely from 
high and narrow to low and flat. Areas under the curve can 
be found in tables in handbooks. In order that a single table can 
suffice for all shapes of curves, the horizontal scale is so made that 
the area under the curve is unity. This is done by using the non- 
dimensional variable ¢ defined as 


t= = (1) 
o 
The value of ¢ thus represents the number of standard errors in a 
given error z. On the ¢-axis the value of @ is equal to unity. 
Tables for area A under the normal curve can be arranged in ~ 
various ways. For example, Fig. 1(b) shows the area from — © to 
the desired value of t. Numerical values for this area can be 
found in Table 1. Other tables give Ao‘ or the area from the origin 
on the ¢-scale to the desired value of ¢. One type of table is about 
as convenient as another when it is remembered that 


+ A- = 0.5 (2) 
Still other tables give A_,' where 
= (3) 


Standard Deviation for Assembly 


When a manufacturing process employs 100 per cent inspec- 
tion that is 100 per cent efficient, the critical dimensions will all lie 
within the limits of the drawings, and the assembly clearances 
will be as anticipated by the designer. The situation is dif- 
ferent when the parts are made in large quantities, and the quality 
is maintained by a system of quality control or acceptance 
sampling. The entire production is now expected to be assembled 
into the finished product despite the obvious fact that a certain 
small percentage of the parts may not be able individually to 
meet the blueprint specifications. 
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The detail dimensions for two parts that are to be assembled 
end to end are shown in Fig. 2. The maximum possible length of 
the assembly is m% + wu + wd + uw. The minimum possible 
length of the assembly is aw — m + a — uw. The maximum 
variation in the length of the assembly is 2u: + 2we. 

Such extreme combinations are of very rare occurrence for 
large-volume manufacture when the production processes are in 
statistical control, when the assembly is random, and when the 
lengths of the individual parts are in no way dependent on each 
other. The theory of probability gives information on the manner 
in which the lengths of the assemblies are more likely to vary. 
This is given by the following equation which is easily proved: 


a, = (012 + (4) 


where o; and 2 are the standard deviations for parts 1 and 2, re- 
spectively, and a, is the standard deviation for the assembly when 
they are placed end to end. An equation similar to the foregoing 
can be written in terms of tolerances if uw is proportional to 
and uw is proportional to oz but not otherwise. If more than two 
bodies are in the assembly, additional terms will appear on the 
right side of Equation (4). 

If two bodies are mated side by side, it can be proved that 
Equation (4) is also applicable to the difference between the 
lengths of the two bodies. 

It is a fundamental theorem of statistics that if each of two 
variables is normally and independently distributed, then the 
sum or difference of these variables will be normally distributed 
with a mean that is the sum or difference of the component 
means, and with a standard deviation given by the foregoing 
Pythagorean equation. 


Dimensioning of Clearance Fit 


The assembly of a hole and shaft can be taken as an example 
of a clearance fit. As shown in Fig. 3, let the diameter of the 
shaft be specified as a + wi, and let the diameter of the hole be 
specified as % + ue The allowance or minimum clearance is a, 
and the maximum clearance is a.. Let the manufacturing proc- 
esses be in statistical control with normal distributions and 
standard deviations o; and oz. The frequency curves for the 
natural tolerances of +30 can be drawn as shown in Fig. 3(a). 

The shaded areas in the tails of the curves represent the pro- 
portions of shafts and holes that lie outside the specified tol- 
erances. For 100 per cent inspection, no assembly will have a 
clearance less than a,, and no assembly will have a clearance 
greater than a2. 

However, let us assume that the entire production of the 
quality indicated is to be assembled. It must be conceded that 
some assemblies will have a clearance less than ai, or will even 
interfere, and that some assemblies will have a clearance greater 
than a2. Methods for determining the proportion of such misfits 
will now be given. 

Let a new variable w be defined as the clearance or the difference 


between the diameters for a random hole-and-shaft assembly. 


The distribution for w, as shown in Fig. 3(b), will be normal with a 
mean # equal to the difference between the means for the hole 
and shaft 


D=h-th (5) 
By Equation (4), the standard deviation o,, for the assembly is 
given by 
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The distribution curve for w can be drawn by choosing a suita- 
ble scale for o,. The clearance between a hole and shaft is 
represented by a point w,; on the w-axis. With w known, points 
such as the origin, ai and az can be located. For example, for 
w; = ai, the assembly is at the minimum clearance specified by 
the drawing. The area under the w-curve to the left of a; gives 
the proportion of assemblies with clearances less than a;. For w; 
= 0, there is no clearance in the assembly. The area of the w- 
curve to the left of w; = 0 gives the proportion of assemblies that 
will interfere. 

Equation (1) can be adjusted for the w-distribution as follows: 

— w 


(7) 


Example 1. As shown in Fig. 3, let the shaft tolerance be 
+0.0015 in., with 7, = 0.0010 in. Hole tolerance is +0.0020 in. 
with o2 = 0.0016 in. Difference of means, w is 0.0045 in. If the 
entire production is accepted for assembly, find the expected pro- 
portion of assemblies with clearance less than the allowance, a; = 
0.0010 in., and the proportion of assemblies that can be ex- 
pected to interfere. 

Solution: 


o,, = (0.00102 + 0.00162)'/* = 0.001887 in. 
By Equation (7) 


_ 0.0010 — 0.0045 _ 


be = 001887 — 1.855 


From Table 1 


= 0.0318 


Hence 3.2 per cent of the assemblies will have clearances less than 
the specified allowance for 0.0010 in. 
For interference 


w; = 0 


By Equation (7) 
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Table 1 Areas under the normal curve 
(Proportion of total area between limits of — © and f;) 


i ty A_ow% & A_ow'i t A_ow & A_a% 
0.00 0.5000 —0.50 0.3085 —1.00 0.1587 —1.50 0.0668 —2.00 0.0228 —2.50 0.0062 
—0.01 0.4960 —0.51 0.3050 —1.01 0.1562 —1.51 0.0655 —2.01 0.0222 —2.51 0.0060 
-0.02 0.4920 —0.52 0.3015 —1.02 0.1539 —-1.52 0.0643 —2.02 0.0217 —2.52 0.0059 
-0.03 0.4880 —0.53 0.2981 —-1.03 0.1515 —-1.53 0.0630 —2.03 0.0212 —2.53 0.0057 
—0.04 0.4840 —0.54 0.2946 —1.04 0.1492 —-1.54 0.0618 —2.04 0.0207 —2.54 0.0055 
-0.05 0.4801 —0.55 0.2912 —-1.05 0.1469 -1.55 0.0606 —2.05 0.0202 —2.55 0.0054 
-0.06 0.4761 —0.56 0.2877 —1.06 0.1446 —-1.56 0.0594 —2.06 0.0197 —2.56 0.0052 
-0.07 0.4721 —0.57 0.2843 -1.07 0.1423 —1.57 0.0582 —2.07 0.0192 —2.57 0.0051 
—0.08 0.4681 —0.58 0.2810 —1.08 0.1401 —1.58 0.0571 —2.08 0.0188 —2.58 0.0049 
-0.09 0.4641 -—0.59 0.2776 —-1.09 0.1379 —1.59 0.0559 —2.09 0.0183 —2.59 0.0048 
-0.10 0.4602 —0.60 0.2743 —-1.10 0.1357 —1.60 0.0548 —2.10 0.0179 —2.60 0.0047 
-0.11 0.4562 —0.61 0.2709 -1.11 0.1335 —1.61 0.0537 —2.11 0.0174 —2.61 0.0045 
-0.12 0.4522 —0.62 0.2676 —1.12 0.1314 —1.62 0.0526 —2.12 0.0170 —2.62 0.0044 
-0.13 0.4483 —0.63 0.2643 —-1.13 0.1292 —1.63 0.0516 —2.13 0.0166 —2.63 0.0043 
-0.14 0.4443 —0.64 0.2611 —1.14 0.1271 —1.64 0.0505 —2.14 0.0162 —2.64 0.0041 
-0.15 0.4404 —0.65 0.2578 —1.15 0.1251 —1.65 0.0495 —2.15 0.0158 —2.65 0.0040 
-0.16 0.4364 —0.66 0.2546 -1.16 0.1230 —1.66 0.0485 —2.16 0.0154 —2.66 0.0039 
-0.17 0.4325 —0.67 0.2514 -1.17 0.1210 —1.67 0.0475 —2.17 0.0150 —2.67 0.0038 
-0.18 0.4286 —0.68 0.2483 —-1.18 0.1190 —1.68 0.0465 —2.18 0.0146 —2.68 0.0037 
-0.19 0.4247 —0.69 0.2451 —-1.19 0.1170 —1.69 0.0455 —2.19 0.0143 —2.69 0.0036 
-0.20 0.4207 —0.70 0.2420 -1.20 0.1151 —-1.70 0.0446 —2.20 0.0139 —2.70 0.0035 
-0.21 0.4168 —0.71 0.2389 —1.21 0.1131 -1.71 0.0436 —2.21 0.0136 —2.71 0.0034 
-0.22 0.4129 —0.72 0.2358 —1.22 0.1112 —1.72 0.0427 —2.22 0.0132 —2.72 0.0033 
-—0.23 0.4090 —0.73 0.2327 —1.23 0.1093 —1.73 0.0418 —2.23 0.0129 —2.73 0.0032 
-0.24 0.4052 —0.74 0.2297 —-1.24 0.1075 —1.74 —2.24 0.0125 —2.74 0.0031 
-0.25 0.4013 —0.75 0.2266 —-1.25 0.1057 —-1.75 0.0401 —2.25 0.0122 —2.75 0.0030 
-0.26 0.3974 —0.76 0.2236 —1.26 0.1038 —1.76 0.0392 —2.26 0.0119 —2.76 0029 
-0.27 0.3936 —0.77 0.2207 0.1020 -1.77 0.0384 —2.27 0.0116 —2.77 0.9028 
-0.28 0.3897 —0.78 0.2177 —1.28 0.1003 —-1.78 0.0375 —2.28 0.0113 —2.78 0.0027 
-0.29 0.3859 —0.79 0.2148 —1.29 0.0985 —-1.79 0.0367 —2.29 0.0110 —2.79 0.0026 
-0.30 0.3821 —0.80 0.2119 —1.30 0.0968 —1.80 0.0359 —2.30 0.0107 —2.80 0.0026 
—0.31 0.3783 —0.81 0.2090 —1.31 0.0951 —1.81 0.0351 —2.31 0.0104 —2.81 0.0025 
-0.32 0.3745 —0.82 0.2061 —-1.32 0.0934 —1.82 0.0344 —2.32 0.0102 —2.82 0.0024 
-0.33 0.3707 —0.83 0 2033 -—1.33 0.0918 —1.83 0.0336 —2.33 0.0099 —2.83 0.0023 
—0.34 0.3669 —0.84 0.2005 —1.34 0.0901 —1.84 0.0329 —2.34 0.0096 —2.84 0.0023 
—0.35 0.3632 —0.85 0.1977 —1.35 0.0885 —1.85 0.0322 —2.35 0.0094 —2.85 0.0022 
—0.36 0.3594 —0.86 0.1949 —1.36 0.0869 —1.86 0.0314 —2.36 0.0091 —2.86 0.0021 
—0.37 0.3557 —0.87 0.1922 —1.37 0.0853 —1.87 0.0307 —2.37 0.0089 —2.87 0.0021 
—0.38 0.3520 —0.88 0.1894 —1.38 0.0838 —1.88 0.0301 —2.38 0.0087 —2.88 0.0020 
—0.39 0.3483 —0.89 0.1867 —1.39 0.0823 —1.89 0.0294 —2.39 0.0084 —2.89 0.0019 
-—0.40 0.3446 —0.90 0.1841 —1.40 0.0808 —1.90 0.0287 —2.40 0.0082 —2.90 0.0019 
—0.41 0.3409 —0.91 0.1814 —1.41 0.0793 —1.91 0.0281 —2.41 0.0080 —2.91 0.0018 
—0.42 0.3372 —0.92 0.1788 —1.42 0.0778 —1.92 0.6274 —2.42 0.0078 —2.92 0.0017 
—0.438 0.3336 —0.93 0.1762 —1.43 0.0764 —1.93 0.0268 —2.43 0.0075 —2.93 0.0017 
-—0.44 0.3300 —0.94 0.1736 —1.44 0.0749 —1.94 0.0262 —2.44 0.0073 —2.94 0 0016 
—0.45 0.3264 -—0.95 0.1711 —-1.45 0.0735 —-1.95 0.0256 —2.45 0.0071 —2.95 0.0016 
—0.46 0.3228 —0.96 0.1685 -1.46 0.0721 —1.96 0.0250 —2.46 0.0069 —2.96 0.0015 
—0.47 0.3192 —0.97 0.1660 —1.47 0.0708 —1.97 0.0244 —2.47 0.0068 —2.97 0.0015 
-0.48 0.3156 —0.98 0.1635 —1.48 0.0694 —1.98 0.0239 —2.48 0.0066 —2.98 0.0014 
-0.49 0.3121 —0.99 0.1611 —1.49 0.0681 —1.99 0.0233 —2.49 0.0064 —2.99 0.0014 
—0.50 0.3085 —-1.00 0.1587 —-1.50 0.0668 —2.00 0.0228 —2.50 0.0062 —3.00 0.00135 
—0.0045 2 where o,, is the larger of the standard deviations, and a, is the 
~~ 0.001887. 385 smaller. Equation (7) then becomes 
w 
‘From Table 1 o,= (9) 
= 0.0086 ty ( 
o, 
‘Hence 0.86 per cent of the assemblies can be expected to interfere. ; 
Equation (6) can be written Dimensioning of Assembly When Standard Deviations Are 


= (7,2 + = 6, ( a) 8) 
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The allowance or minimum clearance, and the tolerances for a 
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clearance fit are specified by the designer at values that will give 
the best engineering functioning of the device. These may be 
modified to some extent by whatever information is available 
about the ease or difficulty of holding a particular tolerance. In- 
formation regarding the standard deviations resulting from 
various manufacturing operations is badly needed by engineering 
offices. However, when such information is lacking, it is still 
possible for the designer to specify approximately the dimen- 
sional quality that must be obtained in the different piece parts. 

It is assumed that the production processes are in statistical 
control, the dimensions are in no way dependent on each other, 
and the entire production is to be utilized by random assembly. 

When values for wi, uz, and a; are chosen, the value for w is 
determined. 

Utilization of the entire production necessarily infers that some 
few parts will be outside the blueprint tolerances so that a small 
percentage of the assemblies will have clearances less than the 
specified value of w; = a. The designer must form an estimate 
of the percentage of such assemblies that can be permitted. The 
table for areas under the probability curve gives the correspond- 
ing value of t,,. 

In other cases the critical condition may be interference upon 
assembly, w; = 0. Again the designer must form a satisfactory 
estimate. 

Although nothing may be known about the individual values of 
the o’s of the two parts, an estimate can be made for the value of 
the ratio a,/o,. Should information be totally lacking, for a 
first approximation, the value of the ratio can be taken as unity. 

With the foregoing information at hand, Equation (9) gives 
the value of o,, the required quality of the production. In gen- 
eral, instead of specifying the o, it is more appropriate to give the 
percentage of parts that must meet blueprint specifications. The 
process is illustrated by the following example. 

Example 2. Let the shaft tolerance wu; be +0.0005 in., and the 
hole tolerance u, be +0 0006 in. Allowance a; is 0.0005 in. If 
2 per cent of the random assemblies can be permitted to have 
clearances less than a, find the percentage of shafts that must 
meet the specified tolerance, and also the percentage of holes. 
Assume ¢,/0, = 02/0, = 1.5. 

Solution: w; = a; = 0.0005 in. 


=u + a = 0.0005 + 0.0006 + 0.0005 = 0.0016 in. 


For A_,,‘” = 0.02, probability table gives t,, = —2.054 
By Equation (9) 


0.0005 — 0.0016 


o2 = 0, = 1.5 X 0.000297 = 0.000446 
0.000500 
0.000297 


By probability table 
= 0.0462 
and 


A-~1.683 *1-83 = 2(0.5 — 0.0462) = 0.9076 


By probability table 


A_,, 1347 = 0.0891 
and 


A_4.347'-347 = 2(0.5 — 0.0891) = 0.8218 
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Note that the production process for shafts need only meet a 
tolerance of +30, or + 0.0009 in. instead of +0.0005 in. as called 
for by the engineering tolerances. For holes, the production 
tolerance is +30, or +0.0013 instead of +0.0006 in. 

If the shaft diameter is the only critical dimension on the draw- 
ing, it should be marked, “‘At least 91 per cent of shafts must meet 
specified tolerances.’’ Similarly, the drawing for the holes should 
be marked, “At least 82 per cent of holes must meet specified 
tolerances.” The use of such markings was suggested years ago 
in an article by Rice.! 

Hence by this method of dimensioning, the engineer retains his 
desired tolerances on the parts, and on the assembly, and at the 
same time an increased tolerance is permitted the production de- 
partments. Accuracy of the results of course will depend on the 
degree of precision with which the value of ratio o,/¢, is known. 


Special Case for o,/o, = uz/u, 


When the values of the standard deviations of the production 
processes are sufficiently well known, the tolerances can be taken 
proportional to them if so desired; that is 


(10) 


where wu, is the larger tolerance and u, is the smaller. 
Let the clearance w,; for the assembly be at the specified mini- 
mum valuea;. Then 


2) (11) 
s 
This should be substituted into Equation (9). 
Us (: 
Us 
ty ) 
u,? 


After transposing o,, we can substitute ¢, for u,/o,. Then 


2\'/ 
Us 
(138) 


Us 


Example 3. Let wz/u, = o,/o0, = 1.5. If 2 per cent of the 
assemblies can be permitted to have clearances less than a speci- 
fied value for a, find the percentage of parts that must meet the 
tolerances of the drawings. 

Solution. For A_,“” = 0.02, probability table gives t, = 
—2.054. 

By Eguation (13) 


(1 + 1.52)'/2 
= war 2.054 = —1.481 


By probability table 
= 0.0693 
A_+.48:!-481 = 2(0.5 — 0.0693) = 0.8614 


Hence 86 per cent of the parts must meet the tolerances. 

The curves in Fig. 4 have been plotted by making computa- 
tions similar to the foregoing. These are convenient and give 
quick results. 


1B. Rice, ‘Setting Tolerances Scientifically,” Mechanical Engi- 
neering, vol. 66, 1944, p. 801. 
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It should be noted in Equation (13) that ¢, depends only on 4, 
and the ratio u,/u,, both of which are independent of a. For 
an assembly with a; equal to zero, these curves then give the 
per cent of interferences. Equation (10), however, must be ap- 
plicable. 


Overlapping Tolerances 


Figs. 3 and 5 illustrate the fact that the normal curves for the 
two members of an assembly can have a considerable overlap with 
a percentage of interferences much less than popularly supposed.? 
By permitting the curves to overlap and accepting the possibility 
of a small percentage of interferences, the values of the o can be 
increased considerably, with a resulting reduction in the cost of 
production. Interference depends on the overlap of the normal 
curves, and is not dependent on any overlap which arbitrarily 
specified tolerance zones may have. 

Example 4. Suppose the engineering requirements indicate 


that the tolerances should be taken as in Fig. 6(a). For o,/o, = ~ 


02/0, = 1.2, and 1 per cent of permissible interferences, find the 
percentage of parts that must meet the tolerances. 

Solution: For interference, w; = 0. 

In Fig. 6(a) 


w = 0.006 in. 
By probability table, for A_,” = 0.01 

= —2.3267. 
By Equation (9) 


2R. W. Hanna and E. C. Varnum, “Interference Risk When Nor- 
mal Distributions Overlap,’ Industrial Quality Control, vol. 7, Sep- 
tember, 1950, p. 26. 
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Interference 
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(a) Overlapping folerances 
w:.006 


ASS LL 

Uz 002 Y-002| | uz=.004 || 
(b) Tolerances revised togive zero 

overlap and zero allowance 


Fig. 6 


0.006 
2.3267(1 + 1.22)'/2 


= 0.00165 in. 


o, = 1.2 X 0.00165 = 0.00198 


It should be noted that the overlap of the tolerances had no effect 
on the values of 0; and 02. 


By probability table 
= 0.0346, 


and 
A-1.81;1'87 = 2(0.5 — 0.0346) = 0.9308 
= = 2.524 
= 0.0058, 

and 


A-~2.54?574 = 2(0.5 — 0.0058) = 0.9884 


Hence 93 per cent of the shafts and 99 per cent of the holes must 
meet the specifications. 

It is questionable, however, whether the detail drawings should 
show an overlap in the tolerances. Such drawings are likely to be 
checked by persons not familiar with the implications of statistical 
theory. They would likely conclude that a large number of inter- 
ferences would occur, and disagreement and confusion result. 
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When the engineering requirements indicate a need for an over- 
lap, the tolerances can be revised in such a way that no overlap 
appears on the drawings. A suitable modification must be made 
in the percentages of parts that must meet the new tolerances. 
The process is illustrated by the following example. 

Example 5. Suppose the tolerances of Fig. 6(a) are revised to 
those shown in Fig. 6(b) where both overlap and allowance are 
zero. If the previous values for W, 0), and o2 are retained, find 
the percentage of parts that must meet the dimensions of the 
revision. 


Solution: 
0.00200 
= 1.2114 
0.00165 
= 0.1130 
and 
A-1.a114!'7114 = 2(0.5 — 0.1130) = 0.7740 
U2 0.00400 
=— = = 2.019 
02 0.00198 
and 


A-s.o0% = 2(0.5 — 0.0217) = 0.9566" 


Hence for the revised tolerances, 77 per cent of the shafts and 96 
per cent of the holes must meet specifications. 


A discussion of truncated normal curves has been given by 
Acton and Olds.® 


DISCUSSION 
L. S. Linderoth, Jr.4 


The author’s paper is not only an interesting but also a very in- 
formative one which is particularly helpful to the individual not 
familiar with the statistical approach to dimensioning machine 
parts. During the presentation of the paper the writer raised the 
question as to whether this particular approach was equally 
applicable where more than two dimensions were involved in a 
given tolerance and finds after reading the paper again that this is 


3F. S. Acton and E. G. Olds, ‘‘Tolerances—Additive or Pythago- 
rean?” Industrial Quality Control, vol. 5, November, 1948, p. 6. 

4 Associate Director of Research, Engineering Department, Con- 
tinental Can Company, Inc., Chicago, II]. Mem. ASME. 
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ably covered and is a very workable system. Every young engi- 
neer entering the machine-design field should have an oppor- 
tunity to read this very excellent paper. 


Lawrence Slote® 


The paper presented is primarily expository in nature. While 
nteresting as a lecture to undergraduates taking a course in 
imachine design, there is nothing new in the presentation which 
can be considered as a material contribution to machine design. 
This paper merely reiterates what can be found in any elementary 
text on statistics or statistical quality control. Examples of this 
may be found in the textbooks by E. L. Grant,® and C. Eisenhart, 
M. W. Hastay, and W. A. Wallis.” 

The basis of the present application to the dimensioning of 
machine parts is the fact that the frequency of occurrence of an 
error versus its magnitude follows the bell-shaped or normal 
curve. In practice a great many manufacturing distributions do 
not follow the normal distribution curve. Some examples are 
radial play, eccentricities, and so on. This being the case, the 
methods shown in this paper are not applicable. Other methods 
can be used such as that discussed by R. L. Thoen.® 

Another application of statistical techniques in machine design 
is well presented by K. H. Moltrecht and R. M. Caddell.® 

What is pointed out in the paper has been clearly stated by 
Prof. Bernard Ostle:° 

“Tf a basic course in statistics is introduced at an early stage 
into all undergraduate engineering curricula, the engineer will at 
least become acquainted with the statistical approach as well as 
with statistical terminology.” 


5 Research Scientist, New York University, College of Engineer- 
ing, Research Division, New York, N. Y. Assoc. Mem. ASME. 

6 E. L. Grant, ‘‘Some Statistical Aspects of Tolerances in Statistical 
Quality Control,” second edition, McGraw-Hill Book Company, Inc., 
New York, N. Y., 1952, chapt. 12. 

7C, Eisenhart, M. W. Hastay, and W. A. Wallis, ‘“Tolerance Limits 
for Normal Distributions in Selected Techniques of Statistical Analy- 
sis for Scientific and Industrial Research and Production and Manage- 
ment Engineering,’’ McGraw-Hill Book Company, Inc., New York, 
N. Y., 1947, chap. 2. 

8R. L. Thoen, ‘‘An Easy Approach to Statistical Tolerancing With 
Punched-Card Computers,’’ Machine Design, vol. 29, 1957, pp. 121- 
123. 

9 K. H. Moltrecht and R. M. Caddell, ‘‘How to Determine Produc- 
tion Tolerances,’ Tool Engineer, vol. 39, October, 1957, pp. 81-85, 
November, 1957, pp. 85-89. 

10 Bernard Ostle, ‘“‘Statistics in Engineering,” Journal of Engineering 
Education, vol. 64, 1957, pp. 410-414. 
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The Dynamic Analysis and Design of 
Relatively Flexible Cam Mechanisms 


R. C. JOHNSON 


Assistant Professor of Mechanical 
Engineering, School of Engineering, 
Yale University, New Haven, Conn. 

Assoc. Mem. ASME 


Having More Than One Degree of Freedom 


The author presents a simple numerical method for analyzing and designing cam sys- 
tems which have more than one degree of freedom. [Illustrative examples are given for 
follower systems having two degrees of freedom. In addition to the general analysis 


and design methods which are derived, the examples also illustrate the importance of 
significant figures in the calculating procedure and the importance of the consideration 


of the distribution of mass throughout elastic regions. 


The general methed of analysis 


and design can be extended to follower systems having more degrees of freedom than 


two. 


Pitenas we have considered the application of 
finite differences to the dynamic analysis and the design of cam 
mechanisms which have single-degree-of-freedom follower sys- 
tems.! Also, it has been pointed out that certain follower 
systems will exist in practice which have more than one significant 
degree of freedom, and by rights they should be treated as such 
in analysis and design work. The basic methods which have been 
developed for the single-degree-of-freedom systems can be ex- 
tended easily to studies of multidegree-of-freedom systems. Asa 
typical example, let us consider a follower system which consists 


Contributed by the Machine Design Division and presented at 
the Annual Meeting, New York, N. Y., November 30—-December 5, 
1958, of Tae American Society or MECHANICAL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, July 10, 
1958. Paper No. 58—A-192. 


Nomenclature. 


of two relatively rigid masses which are bounded by linear elastic 
regions, as depicted in Fig. 1. 


Typical Follower System Having Two Degrees of Freedom 


In Fig. 1 a cam mechanism is shown having a follower system 
consisting of two relatively rigid masses, m, and m,, bounded by 
linear elastic regions. The elastic region which exists between 
the cam surface and the mass m, has a gradient denoted by k,,. 
Another elastic region separates mass m, from mass m,, and its 


1R. C. Johnson, ‘‘The Extension of Finite Differences to the Dy- 
namic Analysis and Design of Relatively Flexible Cam Mecha- 
nisms,”” to be presented by author at the Fifth Conference on 
Mechanisms at Purdue University, October 13-14, 1958; see equa- 
tion (4). Preprints are planned. Will be published in the Trans- 
actions of the Fifth Conference on Mechanisms shortly after the Con- 
ference. 
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acceleration; in/sec? if translation 
motion or radians/sec? if rota- 
tional motion. If a is used as a 
subscript it refers to massa. If 
a is used as a superscript, it 
denotes a constant power which 
helps to define a nonlinear elas- 
tic region.! 

C = a constant which helps to define a 
nonlinear elastic region.’ If is 
used as a subscript it refers to 
the mass c. 

f, = natural frequency at which a fol- 
lower system undergoes free vi- 
brations, cycles per sec. For a 
follower system having two de- 
grees of freedom, f,1 denotes the 
lower mode of free vibration 
whereas f,2 denotes the higher 
mode of free vibration. 

force, lb 

= spring constant for a linear elastic 

region, lb per in. 

m = mass; lb-sec?/in. 

8 = follower displacement which is 


| 
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Aé 


actual 


machined into cam _ surface 
under existence of zero force in 
follower system, inches for trans- 
lation or radians for rotation 


velocity, in/sec for translation or 


rad per sec for rotation 

follower displacement 
which results under dynamic 
operation of actual cam mecha- 
nism, inches for translation or 
radians for rotation 


cam angle interval at which dis- 


placements are specified or cal- 
culated, deg 


angle through which cam rotates 


during one complete cycle of 
free vibration of follower system, 
deg. For a follower system 
having two degrees of freedom 
6,1 denotes the lower mode of 
free vibration whereas 6,2 de- 
notes the higher mode of free 
vibration. 


constant velocity at which the 


cam rotates, rpm 


Subscripts 


a, b, c, ete., refer to masses of follower sys- 
tems having more than one degree of 
freedom. 

0 refers to point under consideration with 
respect to velocities, accelerations, or 
displacements. 

1 refers to point immediately previous to 
the one under consideration, AO deg 
away. 

2 refers to point immediately following the 
one under consideration, A@ deg away. 

s used as a subscript refers to region im- 
mediate to cam surface. 

L refers to a load spring acting on a fol- 
lower system. 

ni and n2 refer to natural frequencies as 
defined previously. 


For multidegree-of-freedom systems, 
follower mass displacements will be de- 
noted by two subscripts; the letter refers 
to the mass under consideration while the 
numeral 0, 1, or 2 refers to the angle as de- 
fined in the foregoing. 
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gradient is denoted by k,,. A load spring might act on top of 
mass m, and its gradient has a value of kz. The follower system 
shown has two natural frequencies for free vibrations and both 
can be calculated quite easily by common methods.? The angles 
through which the cam will rotate during one complete cycle of 
free vibration can be calculated for each mode by applying equa- 
tion (10) of the author’s previous paper! 


(1a) 


(1b) 


where the frequency f,; is for the lower mode of free vibration, 
f,2 is for the higher mode of free vibration, 6, is the angle through 
which the cam will rotate during one complete cycle of free 
vibration of the follower system at the lower mode, and 6p:2 is 
the angle through which the cam will rotate during one complete 
cycle of free vibration of the follower system at the higher mode. 


x 
Xao |” | 


Fig. 1 Schematic representation of a typical cam mechanism with a 
follower system having two degrees of freedom and linear elastic 
regions. Only vertical motion of follower masses permitted. 


Since fn2 is the higher natural frequency at which either of the 
two masses will have a free vibratory component of motion, the 
derivation of equation (12)! for one-degree-of-freedom systems 
can be extended to the same problem of determining cam-angle- 
interval size for follower systems having two degrees of freedom. 
Hence by mathematical induction, the desirable size of cam-angle 
interval (A9) in order to obtain accuracy in the finite-difference 
design and analysis methods is according to equation (2), as 
follows: 


(On2) 


(AA) < 


(2) 


If the cam-angle interval A@ satisfies equation (2), then it will 
generally be relatively small when compared to the variation of 
the total motions for either of the masses, m, or m,. As before, 
in addition to satisfaction of equation (2), accuracy with finite 
differences can only be obtained if the displacements are carried 
to enough significant figures. 


2 J, P. DenHartog, ‘Mechanical Vibrations,’’ McGraw-Hill Book 
Company, Inc., New York, N. Y., 1947. 
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Before proceeding any further, it would be appropriate to de. 
fine clearly the notation which will be used in the subsequent 
derivations and equations. This notation can be extended easily 
to follower systems having more degrees of freedom than two, 
For any displacement of a mass, two subscripts will be used. The 
first subscript will refer to the particular mass under considera- 
tion. The second subscript will designate the time relationship or 
cam-angle relationship for the displacement. Hence the subscript 
0 will denote the instant under consideration, the subscript 1 will 
denote the adjacent displacement value previous to the one under 
consideration A@ deg away, and the subscript 2 will denote the 
adjacent displacement value following the one under considera- 
tion A@ deg away. Therefore, for the mass m,, its displacement 
at the instant under consideration is Zao and at that instant its 
acceleration factor is (%a1 + %a2 — 2%a0). Likewise, for the mass 
mp, its displacement at the instant under consideration is x and 
at that instant its acceleration factor is (21 + 22 — 22). The 
cam displacements will be denoted as before, where so is the cam 
displacement at the instant under consideration and the cam- 
acceleration factor at that instant will be (s: + s2 + 280). There- 
fore equations (2) and (3) of Johnson! can be used directly for 
calculating accelerations and velocities of the masses m, and m,. 


tao = 3 (+) Xai) (3a) 
2 
dao = 36 on + Zaz — 2Ta0) (3b) 


Equations (3a) and (3b) give the velocity and acceleration for 
the mass m, at the instant under consideration. 

Likewise, the following equations give the velocity and ac- 
celeration for the mass m, at the instant under consideration: 


Yy = 3 () — 21) (4a) 


2 
ay = 36 (+) + — 220) (4b) 


The velocities and accelerations for the cam command can be cal- 
culated directly from equations (2) and (3) of Johnson.! For all 
of the equations which have been covered and for those which will 
follow, the positive direction for displacements, velocities, and 
accelerations will be in the direction away from the cam surface. 

Let us now consider a free-body diagram for mass m, from 
Fig. 1. This free-body diagram is shown in Fig. 2. At this time 
it will be appropriate to apply Newton’s second law to the free- 
body diagram in Fig. 2. 


Fig. 2 Free-body diagram 


Fig.3 Free-body diagram 
of mass m from Fig. 1 


of mass mz from Fig. 1 


Transactions of the ASME 


cons 


The 


Pas = 
Ther 
= 

30 E 
Kap 

favo Apt 
ae A 
Sur 
mo 
| 
Lik 
ma 
be 
an 
the 
Fa = Kr 
= Kab + th 
an 
Kis 
01 
Fs Fa m 
Je 


to de. 
quent 
easily 
n two, 
. The 
‘idera- 
hip or 
script 
1 will 

under 
te the 
‘idera- 
ement 
int its 
Mass 
and 
The 

e cam 
cam- 

lhere- 
tly for 
id ™,. 


(3a) 
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DEF = mearo 


F, Fa = mf 30( + — 20m) | 


But from Fig. 1, F, = k,,(so — 20) and from the same figure, 
Fa = kav( — Zao). 
Therefore, 


2 
= mM E (21 + — | (5) 


In a similar manner, a free-body diagram of mass m, could be 
considered as shown in Fig. 3. Applying Newton’s second law 


yor = 


Fa — Fr = madao 


Therefore 


2 
— Zao) — Fy, = Ma E (+) + La: — 


(6) 


Equations (5) and (6) will be used directly to derive the analy- 
sis and design equations which follow. 


Application of Equations (5) and (6) to Follower Analysis 


Assume the existence of a follower system as shown in Fig. 1. 
Also assume that the cam displacements s are given values. 
Suppose that the problem is to determine the responding 
motion for both of the follower masses, mq and m,. 

Solving equation (5) for 2,2 results in the following: 


(7) 


Likewise, solving for xa2 in equation (6) results in the following: 


Ta = 
(8) 
| 


Equations (7) and (8) can be applied directly in a sequential 
manner for analysis problems. The cam-angle interval A@ should 
be small enough for accuracy as suggested by equation (2). The 
analysis should start at a dwell region or any region where both of 
the follower displacements are known at two adjacent points; 
this is equivalent to saying that the initial velocities and displace- 
ments of the follower masses must be known in order to analyze 
the ensuing follower motion, and this requirement is true for any 
analysis method. Hence the displacements 8, 240, 21, Zao, and 
qi Will be known values for the first calculation. Being a simple 
linear function of the displacement x40, the load spring force F, 
can be determined directly each time since Zao is a known value 
for each calculation. Of course, if no load spring exists, Fz is 
merely zero in equation (8). 

For the first follower-response calculation, the displacements 
%2 and 2x2 are determined directly using equations (7) and (8). 
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Of course, once these displacements are found, they become known 
values and we can proceed to calculation of the follower-mass dis- 
placements for the next cam-angle position. This procedure is 
repeated as far as desired, using equations (7) and (8) for each 
calculation; and in this manner the follower displacements 2, 
and z, are determined for each cam angle. After the determina- 
tion of the follower displacements, the acceleration factors for the 
two masses can be found by the normal procedure, (a1 + 22 — 
2250) and (2a; + — 2%a0). In this manner the responding 
motions for both of the follower masses can be found for a given 
cam command. 

A discussion already has been given on the limitations of analy- 
sis and significant figures, for a single-degree-of-freedom follower 
system.! The basic thoughts behind this discussion can be ex- 
tended quite easily to follower systems having more than one 
degree of freedom. Hence, to save space, the discussions will not 
be repeated. Incidentally, the specification of smallness for the 
cam-angle interval A@ given initially by equation (12)! has been 
extended already to follower systems having two degrees of free- 
dom by equation (2) in this paper. 

The analysis method for systems having two degrees of freedom 
will be illustrated in an example which follows. However, before 
this let us consider the application of equations (5) and (6) to cam 
design. 


Application of Equations (5) and (6) to Cam Design 


Assume the existence of a follower system having two degrees 
of freedom, as shown in Fig. 1. Also assume that the desired dis- 
placements x, for the mass m, are given values at equal cam- 
angle intervals, A9. Suppose that it is desired to find the cam 
displacements s which are required in order to obtain the desired 
response for the mass m, under dynamic operation. 

Solving equation (6) for 2 results in the following: 


32 
Ma E 200) + (9) 
+ 
al 


Solving equation (5) for s results in the following: 


2 
mp [0( + te — |+ — Xao) 


8 = + 
kos 


(10a) 


Combining equation (9) with equation (10a) results in equation 
(10d). Equation (106) is sometimes easier to use than equa- 
tion (10a) since the acceleration factors (a1 + a2 + 2%a0) are gen- 
erally known values at the start, as are the load spring-force 
values, since the displacements x, are given values. Therefore 


Fy, 36 \? 
ie ( + Laz — 2La0) 


+ man + tee — + 


(106) 


J 


Using equations (9) and (10a) or (106), the problem of cam design 
becomes a straightforward procedure. The first step is to deter- 
mine the required motion for mass m, by substituting directly in 
equation (9). This is generally very easy since the motion for 
mass m, has been given initially, and therefore the displacements 
x, and the acceleration factors (ta. + Zaz — 2%a0) are known 
values at the start. Also, the load-spring force F', can be calcu- 
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lated at each point, since it is linearly related to the given dis- 
placements zz. The second step is to calculate the required cam 
displacements 8 by direct use of either equation (10a) or (106), 
whichever is more convenient. This is a straightforward step 
also, since the displacements 2, have been calculated already in 
the first step. 

A discussion already has been made on limitations and sig- 
nificant figures for a single-degree-of-freedom follower system.! 
The basic thoughts in this discussion can be extended directly to 
follower systems having two or more degrees of freedom. As an 
example, for design control of vibrations it is assumed that the 
cam-angle interval A@ is small enough to satisfy equation (2) and 
if so, that the required displacements s can be held in manufactur- 
ing to a sufficient degree of accuracy to obtain motion control. 
In order to save space, all of the discussions in the earlier paper! on 
limitations and significant figures will not be repeated here, and 
the extension to follower systems having two or more degrees of 
freedom will be left to the reader. The design method for systems 
having two degrees of freedom wil] be illustrated in an example 
which follows. 

Let us now consider some results from several illustrative ex- 
amples which were solved. 


Illustrative Examples 


Example 1—Calculating Follower Response for a Follower System 
Having Two Degrees of Freedom. For this example assume that a 
follower system exists which has two degrees of freedom and that 
the specified cam command is that of cycloidal motion, DRD 
type. The problem is to determine both the acceleration and dis- 
placement response of the follower masses m, and m,. As de- 
picted in Fig. 4, the follower system has no load spring, the fol- 
lower masses ma and m, have a weight of 0.25 lb, the elastic re- 
gion between m, and m, has a linear force gradient of ka, = 8000 
Ib per in. and the elastic region near the cam surface has a linear 
force gradient of k,, = 8000 lb per in. The cam is to rotate at 
1200 rpm and the cam displacements will be specified according 
to a cycloidal motion, the total throw being 0.2000 in. during 30 
deg of cam rotation. 

To begin with, we will calculate the two natural frequencies of 
the follower system, using common methods.? Hence 


1 1/2 1/2 


Ma 0.25/386 
Likewise, 
1 (kas + 
f= 2( ms ) ~ 0.25/386 
And 


= 561 eps. 


Therefore the lower natural frequency for the follower system is 
determined as follows: 


far fon 7942) [(™ : sor] 


= 347 eps 


Similarly, the higher natural frequency is as follows: 


2 2 2 1/o) 1/2 


= 908 cps 
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Using equation (1a), the lower mode of free vibration will occur 
during a cam-angle rotation On, 


1200 
= = 20.8 deg = 


Using equation (1b), the higher mode of free vibration wil] 
occur during a cam-angle rotation On2, 


1200 
= = 7.93 deg = Onz 


Therefore, for accuracy of the finite-difference equations, the cam- 
angle interval A@ should be at least as small as what is suggested 
by equation (2). Hence, using equation (2), 

Ono 7.93 


= 1.26 deg 


Therefore we will choose a cam-angle interval of A@ = 1 deg, and 
our analysis work should be quite accurate, provided of course 
that the displacements are carried to enough significant figures 
for acceptable accuracy of the acceleration factors. 

It is interesting to note at this point that the two-mass follower 
system for this problem, Fig. 4, has exactly the same total re- 
sultant mass and stiffness characteristics as did the single-degree- 
of-freedom follower system which was used in the previous ex- 
amples.!' For both systems, the total mass is that of a 0.50-lb 
weight and the total elastic stiffness has a force gradient of 4000 
lb per in. for each case. The distribution of mass with elastic 
properties is the important difference between the two systems, 
The natural frequencies for the two-mass system of this example 
are 347 cps and 908 cps, whereas the natural frequency for the 
one-mass system of the previous paper! is 279 cps. Hence it is 
obvious that the manner in which the mass and elastic properties 
are distributed throughout a follower system will have a pro- 
nounced influence on the dynamic behavior of the system. 

In equations (7) and (8) we will now place the following values: 
kts = kay = 8000 Ib/in., w = 1200 rpm, AO = 1 deg, m, = ma = 
0.25/g = 0.25/386. Therefore equations (7) and (8) become as 
follows: 


zie = [80 + Lao — (0.2383) + (2x0 — | (1a) 


re: = (zi — %¢0)(0.2383) + — | (110) 


For this specific example, equations (11a) and (11b) will be used 
to calculate the response for the follower masses mg and mp» a8 
previously explained. 

As a first step, the cam displacements were calculated by exact 
mathematical methods according to the cycloidal motion pre 
viously discussed, and these displacement values are tabulated in 
Table 1. From these displacements, the cam-acceleration factors 
were calculated at each point (s: + se — 289). Plotting these 
acceleration factors resulted in the cycloidal motion between 0 

‘0 deg for the cam command in Fig. 4. Let us now use equa- 
livas (11a) and (116) to calculate the follower response. 

Assume that the cycloidal command starts after a zero velocity 
dwell which exists between the 359-deg point and the 360-deg 
point. Therefore, for both of these points, s = rq = % = 
0.00000 in., since the cam does not command the follower masses 
to move until slightly after the 0-deg point. Hence for the first 
pair of calculations we will consider the 2,2 and the 2q2-values a8 
existing at the 1-deg position of the cam. Therefore, for this first 
calculation, so = Zao = 250 = 0.00000 and 2 = 2%, = 0.00000. 
Therefore, from equation (lla), 22 = [0.00000 + 0.00000 — 
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Table 1 Follower-response calculations for Example 1 


(Refer to equations (11a) and (116) for response calculations) 


Cycloidal : 
motion 
0 given 
deg 8 Xb Lo (tar + a2 — (201 + 2a — 22X00) 
359 0.00000 0.00000 0.00000 .00000 0.00000 
0 0.00000 0.00000 0.00000 0.00000 0.00000 
1 0.00005 0.00000 0.00000 0.00000 +0.00001 
Z 0.00039 0.00001 0.00000 0.00000 +0.00009 
3 0.00129 0.00011 0.00000 +0.00003 +0.00026 
4 0.00301 0.00047 0.00003 +0.00011 +0.00050 
5 0.00577 0.00133 0.00017 +0.00028 +0.00078 
6 0.00973 0.00297 0.00059 +0.00057 +0.00104 
7 0.01501 0.00565 0.00158 +0.00097 +0.00126 
8 0.02168 0.00959 0.00354 +0.00144 +0.00144 
9 0.02973 0.01497 0.00694 +0.00191 +0.00160 
10 0.03910 0.02195 0.01225 +0.00231 +0.00178 
11 0.04968 0.03071 0.01987 +0.00258 +0.00194 
12 0.06129 0.04141 0.03007 +0.00270 +0.00204 
13 0.07372 0.05415 0.04297 +0.00266 +0.00200 
14 0.08672 0.06889 0.05853 +0.00247 +0.00178 
15 0.10000 0.08541 0.07656 +0.00211 +0.00137 
16 0.11328 0.10330 0.09670 +0.00157 +0.00081 
17 0.12628 0.12200 0.11841 +0.00086 +0.00016 
18 0.13871 0.14086 0.14098 —0.00003 —0.00048 
19 0.15032 0.15924 0.16352 —0.00104 —0.00111 
20 0.16090 0.17651 0.18502 —0.00203 —0.00169 
21 0.17027 0.19209 0.20449 —0.00296 —0.00225 
22 0.17832 0.20542 0.22100 —0.00371 —0.00275 
23 0.18499 0.21600 0.23380 —0.00425 —0.00315 
24 0.19027 0.22343 0.24235 —0.00451 —0.00340 
25 0.19423 0.22746 0.24639 —0.00452 —0.00341 
26 0.19699 0.22808 0.24591 —0.00425 —0.00316 
27 0.19871 0.22554 0.24118 —0.00373 —0.00266 
28 0.19961 0.22034 0.23272 —0.00295 —0.00199 
29 0.19995 0.21315 0.22134 \ —0.06195 —0.00120 
30 0.20000 0.20476‘ 0.20795 —0.00076 —0.00037 
31 0.20000 0.19600 ~“ 0.19383 « +0.00052 +0..00044 
32 0.20000 0.18768 0.18023 
w 
3 10 | 
| 
5 
< FL =0 | 
< COMMAND 
M,= 0.287 386 
Kay = 6000 «RESPONSE 
20 b= 0.257 386 
5 Kus = 8000 
3 200 RPM A 
} 4 8 e 0 20 2 2 32 
5 4 i > + 24 26 32 Fig. 5 C d displacement curves for a cam-driven two- 


1 20 
Cam (DEGREES) 


Fig. 4 C d-resp acceleration curves for a cam-driven two- 
mass system, Example 1. Illustrates a response from a cycloidal com- 
mand. 


2 X 0.00000] (0.2383) + (2 x 0.00000 — 0.00000) = 0.00000 = 
a. Likewise, from equation (11b), 242 = 0.00000. 

For the second pair of calculations, consider the x2 and the 
tervalues as existing for the 2-deg position of the cam. There- 
fore, for the 1-deg position, 8) = 0.00005, 20 = 20 = 0.00000, 
and at the 0-deg position 2; = za: = 0.00000. Therefore, from 
equation (11a), x42 = [0.00005 + 0.00000 — 2 x 0.00000] (0.2383) 
+ (2 x 0.00000 — 0.00000) = 0.00001 in. = 2. Likewise, 
using equation (11b), x42 = 0.00000. 

For the third pair of calculations, consider the 22 and the a:- 
values as existing for the 3-deg position of the cam. Therefore, 
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mass system, Example 1. Illustrates a response from a cycloidal 
command. 


at the 2-deg position s) = 0.00039, 2 = 0.00001, and ra) = 
0.00000, and for the 1-deg position z = 0.00000 = 2a. 
Hence, from equation (11a), 


Ze = [0.00039 + 0.00000 — 2 X 0.00001](0.2383) 
+ (2 X 0.00001 — 0.00000) = 0.000088 + 0.00002 = 0.00011 


Likewise, using equation (11b), z,2 = (0.00001 — 0.00000) 
(0.2383) + 2 x 0.00000 — 0.00000) = 0.00000. 


As has now been illustrated, calculations are made in pairs 
using equations (11a) and (11b). The illustrated procedure is re- 
peated in the sequential manner until the 32-deg point is reached, 
and the resultant response is tabulated in Table 1. Of course, 
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the analysis could be extended beyond the 32-deg point if so de- 
sired. It is suggested that the reader refer to Table 1 in checking 
the analysis procedure. Being a very systematic procedure, the 
method of analysis using Table 1 becomes almost automatic 
after but a few calculations and, therefore, the response calcula- 
tions can be made very rapidly using a desk calculator. As a 
final step, the acceleration factors for bodies m, and m, were 
calculated directly from the calculated displacements; i.e., (721 + 
%a2 — 2240) and (Ze 220), respectively 

From Fig. 4 it is obvious that the follower mass m, follows the 
cam command a little closer than does the follower mass ma, and 
this is probably what would be expected from a consideration of 
the follower system as depicted in Fig. 4. The important thing 
to notice is that the response differs far from the command, even 
though the command is the dynamically accepted cycloidal 
motion, A plot of displacements for this problem is given in Fig. 
5, and displacementwise the follower masses undergo something 
quite different from cycloidal motion. It is obvious that after 
the throw is complete, vibrations of large amplitude will result 
for both of the follower masses, and this can be seen readily from a 
study of Figs. 4 and 5. These ensuing vibrations after the throw 
period of the cam could have been determined specifically by ex- 
tending the analysis beyond the 32-deg point. 

An interesting point to notice from Figs. 4 and 5 is that the 
predominant mode of vibration seems to be that of the lower 
mode. This is generally the case for a very smoothly varying 
command, as we have with the cycloidal cam. 

Example 2—Designing a Cam to Give a Desired Response for a Fol- 
lower System Having Two Degrees of Freedom. For this example 
assume the existence of a follower system which is identical to 
that of Example 1. However, assume that the problem is to de- 
sign a cam which will give the desired cycloidal motion for mass 


m,. Hence this will be an example of tuned cam design for 
follower systems having more degrees of freedom than one. 

The natural frequencies for the follower system can be found 
in exactly the same way as they were for Example 1. Hence, from 
Example 1, fu = 347 cps and f,z = 908 cps. Therefore, as before, 
one cycle of the lower mode of free vibration will occur during a 
cam-angle rotation 0, = 20.8 deg by equation (1a) and one cycle 
of the higher mode of free vibration will occur during a cam-angle 
rotation A. = 7.93 deg by equation (1b). Hence, for vibrational 
control in design of the cam, equation (2) shows that A@ should be 
at least less than 1.26 deg. Therefore, as before, we choose 
= 1 deg. 

For the cam design in this example we will apply equations (9) 
and (10b) directly. As in Example 1, ms = m, = 0.25/386, 
kab = ky, = 8000 lb per in., w = 1200 rpm, and Aé = 1 deg, 
By using these values, equations (9) and (10b) become as follows 
for this particular example: we ia 


tho = (4.197)(ra1 + Lae — 2240) + Lao (12a) 


8 = (4.197) [(2a1 + Laz — 2240) + 


(126) 
(xo + — 2%0)] + 


Instead of applying equations (12a) and (126) in pairs, as we 
did with equations (11a) and (11b) for the analysis Example No. 
1, we will first apply equation (12a) to all of the points giving us 
all of the z,-values, since all of the z,-values are given. Then 
we will calculate the acceleration factors for the mass m, by 
direct usage of the calculated displacements, (2: + 22 — 22). 
As a final step we will calculate the required cam displacements s 


Table 2 Cam-displacement calculations for Example 2 
(Refer to equations (12a) and (12b) for design calculations) 


6, Desired 
deg Le + — (xe + — 2200) 8 (81 + 82 — 28» 
358 0.000000 0.000000 0.00000 0.00000 0.00000 0.00088 
359 0.000000 0.000000 0.00000 +0.00021 0.00088 +0.00219 
0 0.000000 +0.000049 0.00021 +0.00084 0.00395 —0.00265 
1 0.000049 +0.000289 0.00126 +0.00045 0.00437 +0.00233 
2 0.000387 +0.000565 0.00276 +0.00047 0.00712 +0.00103 
3 0.001290 +0.000819 0.00473 +0.00065 0.01090 +0.00056 
4 0.003012 +0.001033 0.00735 +0.00085 0.01524 +0.00047 
5 0.005767 +0.001205 0.01082 +0.00099 0.02005 +0.00049 
6 0.009727 +0.001323 0.01528 +0.00108 0.02535 +0.00066 
a 0.015010 +0.001384 0.02082 +0.00112 0.03131 +0.00079 
8 0.021677 +0.001383 0.02748 +0.00114 0.03806 -+0.00054 
9 0.029727 +0.001323 0.03528 +0.00108 0.04535 +0.00071 
10 0.039100 +0.001205 0.04416 +0.00098 0.05335 +0.00052 
ii 0.049678 +0.001034 0.05402 +0.00084 0.06187 +0.00058 
12 0.061290 +0.000818 0.06472 +0.00067 0.07097 +0.00039 
13 0.073720 +0.000565 0.07609 +0.00047 0.08046 +0.00016 
14 0.086715 +0.000290 0.08793 +0.00023 0.09011 +0.00024 
15 0.100000 0.000000 0.10000 0.00000 0.10000 0.00000 
16 0.113285 —0.000290 0.11207 —0.00023 0.10989 —0.00024 
17 0.126280 —0.000565 0.12391 —0.00047 0.11954 —0.00016 
18 0.138710 —0.000818 0.13528 —0.00067 0.12903 —0.00039 
19 0.150322 —0.001034 0.14598 —0.00084 0.13813 —0.00058 
20 0.160900 —0.001205 0.15584 —0.00098 0.14665 —0.00052 
21 0.170273 —0.001323 0.16472 —0.00108 0.15465 —0.00071 
22 0.178323 —0.001383 0.17252 —0.00114 0.16194 —0.00054 
23 0.184990 —0.001384 0.17918 —0.00112 0.16869 —0.00079 
24 0.190273 —0.001323 0.18472 —0.00108 0.17465 —0.00066 
25 0.194233 —0.001205 0.18918 —0.00099 0.17995 —0.00049 
26 0.196988 —0.001033 0.19265 —0.00085 0.18476 —0.00047 
27 0.198710 —0.000819 0.19527 —0.00065 0.18910 —0.00057 
28 0.199613 —0.000565 0.19724 —0.00047 0.19287 —0.00101 
29 0.199951 —0.000289 0.19874 —0.00045 0.19563 —0.00234 
30 0.200000 —0.000049 0.19979 —0.00084 0.19605 +0.00265 
31 0.200000 0.000000 0.20000 —0.00021 0.19912 —0.00219 
32 0.200000 0.000000 0.20000 0.00000 0.20000 —0.00088 
33 0.20:000 0.000000 0.20000 0.00000 0.20000 0.00000 
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by using equation (125) directly. The cam-acceleration factors 
(s: + 82 — 280) can then be found from these cam displacements 
by straightforward calculation. The design method is so straight- 
forward that a specific numerical example will not be given. 
However, the complete results are given in Table 2 for those who 
might be interested and all of the acceleration factors are plotted 
in Fig. 6. It is interesting to note that the response of follower 
mass mq is cycloidal as desired, and that the accelerations for 
follower mass m, are of low magnitude also. The only accelera- 
tion which varies to a high amplitude is that of the cam command, 
and that is unimportant for this example since negligible mass 
exists in the elastic contact region. Even at that, the accelera- 
tion for the tuned design cam command in Fig. 6 does not vary to 
as high a magnitude as does the untuned design response of fol- 
lower mass m, in Fig. 4. The displacement variations for the 
tuned design are shown in Fig. 7, where the follower response of 
m, Will be that of the desired cycloidal motion. 

Before concluding this example, the importance of significant 
figures in design probably should be illustrated for follower sys- 
tems having two degrees of freedom. We will assume that cam 
displacements s can be machined to the nearest 0.0001 in. Now, 
suppose that the designer calculated the desired response values 
z, to four significant figures. Then each z,-value might be off as 
much as 0.00005 in. Therefore, from equation (12a), each of the 
calculated 2-values might be off as much as (4.197)[+0.00005 + 
0.00005 — (+2 X 0.00005)] + 0.00005 = +0.00079 in. There- 
fore, from equation (125), each of the calculated s:-values might be 
off as much as (4.197) [(+0.00005 + 0.00005 + 2 x 0.00005) + 
(¥0.00079 0.00079 2 0.00079)] + 0.00079 = 0.01163 
in. In summary, therefore, if the desired z,-values are initially 
calculated to only four places, the required s-values will only be 
accurate to the nearest +0.01163 in. Of course, the probability 
of unfavorable errors occurring for all points in a calculation at 
one time is very low, but the possibility still exists. Hence, if 
we calculate cam-displacement values s which are only accurate to 
+0.01163 in., we are certainly not taking full advantage of the 
machining capabilities which are available. 

On the other hand, suppose the desired z,-values are initially 
calculated accurate to six decimals, Hence any one z,-value 
might be off as much as +0.0000005 in. Therefore, from equa- 
tion (12a), each of the determined z,-values will be accurate to 
(4.197)( +0.0000005 + 0.0000005 + 2 x 0.0000005)  0.0000005 
= +0.0000079 in. Asa next step, suppose that we round off these 
%;-values to 5 decimals. Then any one 2,-value might be off as 
much as +0.0000079 + 0.000005 = +0.0000129 in. Therefore, 
using equation (12b), we find that the so-values might be off as 
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much as (4.197)[(+0.0000005 + 0.0000005 + 2:0.0000005) + 
(F0.0000129 + 0.0000129 + 2 x 0.0000129)} + 0.0000129 = 
0.000195 ~ 0.0002 in. Therefore, under this assumption, 
the cam-displacement values calculated in design will be accurate 
to +0,0002 in. at the worst possible condition. Therefore, if the 
initial desired follower displacements z, are calculated to six 
decimals we will be taking more advantage of our available 
machining capabilities, since then our calculated cam displace- 
ments will be near in accuracy to the limitations of machining 
capabilities. Of course, if greater accuracy for the cam displace- 
ments was desired, the displacements x, could be specified to 
seven or more places but a certain amount of effort would then be 
wasted since the required cam displacements would be apprecia- 
bly more accurate than what could be obtained by common 
machining methods. 

Example 3—Cam Design—Showing Effects of Distribution of Both 
Mass and Elastic Regions. As a final example assume that a two- 
mass follower system actually exists as depicted in Fig. 8. For 
this follower system F, = 0, m, = m, = 0.25/386, k,, = ky, = 
8000 Ib per in., and the cam rotates at 1200 rpm. Hence, as in 
Examples 1 and 2, the two natural frequencies for free vibrations 
are again fn = 347 cps and fre = 908 cps. Therefore, as cal- 
culated in Examples 1 and 2, one complete cycle of the lower 
mode for free vibrations would occur during a cam-angle rotation 
6,1 = 20.8 deg, and one complete cycle of the higher mode for 
free vibrations would occur during a cam-angle rotation 0,2. = 
7.93 deg. Hence, from equation (2) we again choose A@ = 1 deg 
which should give us good accuracy for both analysis and design, 
providing of course that the displacements are carried to enough 
significant figures. Let us assume that the designer wishes to 
obtain a linear type of acceleration, and that three different 
methods of solution will be used. The results from each method 
of solution will be illustrated in this example. 

First, assume that the designer neglects follower deflections 
and therefore the cam command is that of the desired linear 
acceleration as shown in Fig. 8. From the given linear-accelera- 
tion cam displacements, the follower response was determined for 
both masses m, and m, using equations (7) and (8), the method 
having been illustrated already in Example 1. The results of 
the analysis are shown in Figs. 8 and 9. From these figures it 
is obvious that the response of mass m, is far from that which was 
desired, both accelerationwise and displacementwise. From 
the acceleration plot for mass m, in Fig. 8 it seems that the 
deviation from command is primarily influenced by the lower — 
mode of free vibration, although a component at the higher mode 
is apparently having some effect also. The important thing to 
notice is that neglect of significant flexibilities in the cam design, 
i.e., untuned design, has resulted in an actual follower response 
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mass system, Example 3. Illustrates typical response from constant- 
pulse command. 
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mass system, Example 3. Cam was designed to give constant-pulse 
responses under the false assumption that a single concentrated mass 
follower existed. 
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Fig. 12 C d-resp acceleration curves for a cam-driven two- 
mass system, Example 3. Cam properly designed to give constant- 
pulse-follower response. 


which differs appreciably from the desired motion, both with re- 
gard to acceleration and displacement. 

Secondly, assume that the engineer designs his cam to obtain 
the desired follower response by the methods of the earlier paper,! 
but that he assumes the follower system to have only one degree of 
freedom consisting of a 0.5-lb weight having a 4000 lb per in. elas- 
tic region of contact. Assume that he obtained the 0.5-lb weight 
value by weighing the actual follower system and that he obtained 
the 4000-lb-per-in. elastic-region value by measuring the total de- 
flection of the follower system for various static contact forces. 
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pulse-follower response. 


These values could have been obtained from the actual two-mass 
system by static-condition measurements. In other words, the 
designer has made an error, in that he has failed to recognize the 
proper distribution of mass and elasticity in the follower system. 
Assume, therefore, that the supposedly tuned cam design was 
made using equation (9b)! under the belief of a single-mass follower 
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system, and the cam displacements were obtained in that manner. 
From these cam displacements, the actual response of the existing 
two-mass system was calculated using equations (7) and (8) of 
this paper, and the results are plotted in Figs. 10 and 11. 

From Fig. 10 notice that the actual acceleration response de- 
viates appreciably from the desired response, the oscillations 
being rich in the higher mode for free vibrations. However, from 
Fig. 11 it is seen that the actual displacement of follower mass m, 
deviates only slightly from the desired displacement response of 


follower mass m,. In other words, conducting a tuned cam de- 


sign under the false belief of having a lumped follower mass re- 
sulted in a much better actual follower-mass displacement re- 
sponse than neglecting flexibilities entirely in an untuned design. 
However, it is questionable as to which of the methods con- 
sidered so far gives the better acceleration response. 


Thirdly, the cam design was carried out properly using equa- 
tions (9) and (10) of this paper, resulting in the exact desired 
response for m, as shown in Figs. 12 and 13. 

In conclusion, we might make the observation that most flexible 
follower systems can be approximated by lumping one, two, or 
three masses bounded by elastic regions. However, in reality 
some systems will have appreciable mass distributed throughout 
an elastic region, resulting in higher modes of free vibration than 
for the approximated case. Therefore, if a tuned design is made 
under the approximated assumption of lumped masses, high-fre- 
quency appreciable acceleration oscillations probably will result. 
However, the follower displacements generally will be very close 
to that designed for and oftentimes follower displacements are of 
primary importance, providing that accelerations do not deviate 
to an unreasonable extent. Hence it is oftentimes much better 
to conduct a tuned cam design under the approximate assumption 
of lumped masses than it is to conduct an untuned design com- 
pletely neglecting follower-system flexibilities, 

Other Multidegree-of-Freedom Follower Systems. Up to this point 
we have only considered follower systems having one degree of 
freedom or at the most two degrees of freedom. For the latter 
we have only considered the existence of linear elastic regions 
bounding the two masses. The basic method of analysis and de- 
sign using finite differences is definitely not limited to such situa- 
tions only. In fact, the possible combinations of conditions be- 
come too numerous for complete coverage in this paper. Hence 
the engineer should be thoroughly familiar with the basic method 
of derivation, which has so far been illustrated in deriving equa- 
tions (7), (8), and (9) of the earlier paper,! for the single-mass 
systems and equations (5), (6); (7), (8), (9), and (10) of this 
paper for the two-mass systems. 


For follower systems having more degrees of freedom than two, 
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merely extend the two-mass-system derivation method; that is, 
extend the same notation for masses and displacements by using 
more letters for subscripts and draw a free-body diagram for 
each mass. Apply Newton’s second law to each free-body dia- 
gram, 2F = ma. Using finite differences, express all forces, 
accelerations, and velocities in terms of displacements. Hence 
basic equations which are comparable to equations (5) and (6) 
will be obtained. For analysis problems, where the cam dis- 
placements s are given, solve these basic equations for 2,2, 12, 2,2, 
and so on, and apply these sequentially as before to determine the 
follower response. Of course it is assumed that the cam-angle in- 
terval A@ satisfies equation (2), except 0, should be the angle 
through which the cam rotates for one cycle of free vibration at 
the highest natural frequency of the follower system. For design 
problems where the z,-values are given, solve the basic equations 
which are comparable to equations (5) and (6) for x40, 2,0, . . .) 80. 
Using these derived equations, calculate all of the 2,-values 
first, all of the z,-values next, and so on, and then all of the cam- 
displacement values s last. 

For follower systems having a nonlinear region of contact, 
merely repeat the previous steps except in the proper free-body 
diagram use equation (4) from the earlier paper! for expressing 
the elastic region force between the cam surface and its adjacent 
mass. In this equation use the proper letter subscript for the 
%-term, such as x. If a nonlinear elastic region exists between 
any other two bodies in the follower system merely include the 
effect by applying equation (4)! to express the nonlinear forces 
in the proper free-body diagrams. Thus, as an example, if the 
elastic region between masses m, and m, in Fig. 1 were nonlinear, 
the force in the free-body diagram of mass m,, Fig. 3, would be 
Fa, = C(x — 2a0)*. Similarly, the force in the free-body dia- 
gram of mass m,, Fig. 2, would be expressed by the same relation- 
ship. Hence in applying Newton’s second law, the nonlinear 
forces would be expressed in such a manner. 

For any of the extended cases to be considered, use the same 
general notation as before. For instance, consider that 89 = tao = 
Zo = Xo, and so on, only for the case where zero force exists in 
the follower system; i.e., zero deflections. As before, the method 
for including load-spring effects is obvious. Initially, during a 
short dwell say, the x,, x, and x,-values will differ from the s- 
values due to the static elastic-region deflections induced by the 
load spring. Under dynamic operation, the load-spring force is 
generally a linear function of the displacements 2a. The 2a0- 
values are generally either given values in design problems or 
else known values when solving for zq:-values in analysis prob- 
lems. Hence the load-spring force can be taken into account 
easily during the cam design or follower analysis procedure as it 
was before. 
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Power Heat Balance Considerations in 
Design and Operation of Industrial Plants 


Relations between design of process equipment and power-service facilities for industrial 
plants and the effect of each on over-all investment and operating costs are shown by 
typical examples. A brief discussion is given of a method of controlling manufactur. 
ing costs by application of power-performance standards once a plant is in opera- 


tion. 


PERFORMANCE of a manufacturing enterprise 
often is measured in terms of return on investment. Many fac- 
tors have a bearing on performance but two of the most significant 
ones are investment cost and operating cost. In this paper the 
design and operation of “‘power-service’’ facilities for industrial 
plants will be considered in relation to the problem of minimizing 
over-all plant investment and operating cost. 

The term power services refers to steam, electricity, refrigera- 
tion, compressed air, water of different qualities, and other 
utilities supplied by modern industrial power plants. 

Often it is profitable to accept an increase in investment or 
manufacturing cost in one area to obtain a cost reduction or an 
advantage of greater magnitude in another area. For instance, 
the higher cost per Btu of electricity as compared with other heat 
sources can be justified in certain industrial-furnace applications 
because of lower investment cost and improvements in product 
quality. 


Design Objectives 


In designing an industrial plant for maximum economy or re- 
turn on investment, it is logical to start with a layout which will 
obtain minimum investment consistent with the basic require- 
ments of safety, working conditions, and product quality. Any 
deviation from this base that will increase the over-all return on 
investment without sacrificing standards of safety, working 
conditions, or quality, may be acceptable for incorporation in the 
final project design. An example of this type would be an ex- 
penditure of additional capital for automatic control systems to 
secure a substantial reduction in operating labor cost. 

After a decision has been made to build a plant, it is usually 
good business to complete the job as quickly as is practical since 
no return on the investment can be realized until the manu- 
factured product is on the market. Often, as in the case of manu- 
facturing a new product, added impetus is given by the prospect 
of becoming established in the market before competition enters 
the field. 

With these considerations in mind, project schedules usually 
are shortened toward the practical limit which requires that much 
of the basic data development, actual design, procurement, and 
construction be carried on simultaneously. A time schedule for a 
typical new plant project is shown on Fig. 1. 

During the early stages of basic data development and design, 
numerous economic evaluations are needed to establish optimuin 
conditions for highest return on investment. In consideration of 

Contributed by the Process Industries Division and presented at 
a joint session with the Power Division at the Annual Meeting, 
New York, N. Y., November 30—-December 5, 1958. of THE 
AmeRICAN Society of MECHANICAL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 


the Society. Manuscript received at ASME Headquarters, Septem- 
ber 12,1958. Paper No. 58—A-166. 
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Fig. 1 Typical project schedule 


the advantage of completing a project in a minimum of time, 
minor sacrifices in accuracy and thoroughness of the economic 
evaluations can be tolerated. 


Economic Evaluations 


Now let us turn to several examples of economic evaluations 
involving power facilities. 

Steam for process and building heating normally is supplied 
from a central power plant through a distribution system. 
An alternative system might include semiautomatic package- 
type boilers located near the points of heat consumption, with a 
minimum investment in steam-distribution lines. Still another 
alternative might involve use of direct-fired heaters, burning 
gas or oil, with no intermediate fluid, for transmitting the heat 
energy. Under certain circumstances it is even economical to 
use electricity as a source of or as a means of distributing heat 
energy. 

Nylon Processing Plant. A nylon yarn plant requires heat at a 
temperature of about 500 F for melting the nylon polymer. Dow- 
therm A, an eutectic mixture of diphenyl and dipheny] oxide, is 
used as the medium for transmitting this heat, since its tempera- 
ture-pressure relation is well suited to this application. In the 
first commercial nylon plant, Dowtherm vaporizers were centrally 
located and heated with electricity. Because of the relatively 
extensive Dowtherm distribution system and the high operating 
temperature, radiation losses amounted to about 80 per cent of 
the total heat load of the system, only about 20 per cent being 
used for melting the nylon polymer. In subsequent installations 
a number of small Dowtherm vaporizers were located at the 
points of use, each being supplied with electricity. As a result, 
substantial reductions in both investment and operating cost were 
accomplished through elimination of Dowtherm distribution pip- 
ing with its heavy radiation loss. 

In later installations a change was made back to a central sys- 
tem with an oil-fired Dowtherm vaporizer and a distribution sys- 
tem, using No. 6 fuel oil as the heat source instead of electricity. 
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The justification for this reversal was the much lower operating 
cost and lower allocated investment for power facilities with direct 
oil firing, in spite of high radiation loss from the distribution 
system. 

The economic studies which provided the bases for decisions in 
these cases involved many considerations. In the first case, a 
new, spectacular product was about to be manufactured commer- 
cially for the first time. The pioneering work, done on a rela- 
tively small scale, had utilized electricity as the heat source, and 
many difficult problems of control had been well worked out. 
The incentive for getting a large-scale plant into production as 
quickly as possible was strong enough to justify proceeding with 
the initial arrangement without risking the long delay which 
might have been experienced in developing a more economical 
system. 

Process and Space Heating. Even problems which look simple 
sometimes have a suprising number of angles. For a plant with 
on-site electric power generation, selection of optimum pressure 
levels of steam for process and space heating requires careful 
analysis. One of the important considerations is the size and 
physical layout of the distribution system. When the steam is ex- 
tracted or exhausted from turbines or engines, the lowest exhaust 
pressure level usually results in minimum operating cost. Witha 
low-pressure system, however, friction drop in the distribution 
piping must be held low and pipeline sizes, therefore, must be 


To illustrate these relations, let us assume that we need to dis- 
tribute 120,000 lb per hr of dry saturated steam for a distance of 
1000 ft and that three pressure levels, 10, 20, and 50 psig are being 
considered. Reasonable allowable pressure drops would be 2 psi 
for the 10-psi system, 5 psi for the 20-psi system, and 10 psi for the 
50-psi system. For these assumed conditions the 10-psi system 
would require a 24-in. pipe, the 20-psi system an 18-in. pipe, and 
the 50-psi system a 14-in. pipe. Schedule 10 pipe would be used 
in all cases. 

Another important item for consideration is that process heat- 
exchanger surface must be greater with a lower steam pressure be- 
cause of the lower delta ¢. For an economic analysis of which 
pressure to use, 10, 20, or 50 psig, the following items would be 
significant. 

Advantages of Lower Pressure. Lower operating cost. This in- 
volves consideration of the over-all power heat balance and the 
over-all effect on boiler steam load. Corollary considerations 
may include: 


Lower allocated investment costs for: 


(a) Smaller boiler capacity. 

(b) Smaller high-pressure steam leads. 

(c) Smaller boiler auxiliaries. 

(ad) Smaller turbine condensers, due to reduction in amount of 
steam condensed. 

(e) Smaller condenser circulating system, including pumps, 
motors, piping system. 

(f) Smaller auxiliary loads, with lower over-all load on elec- 
tric generation and distribution system. 


The amount of consideration which should be given these 
corollary items depends on the size of the plant and how closely 
the factors must be weighed in order that a clear decision can be 
made. 

Disadvantages of Lower Pressure. Higher cost of steam distribution 
system, including piping, valves, fittings, and flow controllers. 

Solvent Recovery. Another example of economic evaluation in- 
volving power-service facilities deals with recovery of solvent from 
@ process air stream. As shown in the flow diagram for one ar- — 
rangement, Fig. 2, air containing an expensive solvent is drawn 
from process by steam-turbine-driven blowers which discharge into 
water-absorption towers at about 2 psig pressure. Air from the 
towers, essentially free of solvent, exhausts to atmosphere. The 
absorption medium is cool water which enters the towers at about 
60 F. From the towers the water and solvent are pumped 
through feed preheaters to distillation columns where direct steam 
is added and the solvent is vaporized. Water and condensate 
pass back through the feed preheaters and then to the sewer or to 
refrigeration units to be cooled and recirculated. The solvent 
vapor is condensed and returned to the system. 

An alternative arrangement for recovering the solvent, as 
shown in Fig. 3, involves operating the absorption towers at about 
75 psig pressure and use of a gas-turbine expander for recovery of 
the power of compression. As shown in the diagram, air contain- 
ing solvent would be drawn from process by a compressor which 
would be driven by the gas-turbine expander. The air-solvent 
mixture would be discharged from the compressor, passed through 
a heat exchanger and to absorption towers operating at about 75 
psig, instead of 2 psig pressure. From the towers the air would 
be passed through a regenerator to recover heat from the turbine- 
expander exhaust, thence to a combustion chamber to be heated 
directly by combustion of gas or fuel oil, and then to the turbine 

expander. The gas-turbine expander would recover a substantial 
portion of the power absorbed by the compressor and the de- 
ficiency would be supplied by the gas or fuel oil. Water at about 
85 F would be the absorption medium. Other components of the 
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Fig. 2 Solvent-recovery system initial arrangement 
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Fig. 4 Process waste-heat boiler and power-plant flow diagram 


system would be essentially the same as in the previous example. 
The following advantages would be realized with the higher 
pressure system: 


1 Elimination of the need for cool water, using water at 85 F 
instead of 60 F. 

2 Reduction in physical size and cost of the absorption 
towers. 

3 Reduction in water flow to the absorption towers. 

4 Reduction in steam consumption of the distillation column. 


A complete evaluation of these two arrangements indicates the 
possibility of reducing investment approximately 4 per cent and 
operating cost approximately 2 per cent through use of the higher 
pressure system. These advantages, however, are hardly large 
enough to justify the additional operating complications which 
would’ be involved, and it is likely that final design would incor- 
porate the simpler, low-pressure system. 

Waste-Heat Boiler Application. One more example of economic 
comparison concerns steam generated in a process waste-heat 
boiler, augmented by steam from coal-fired boilers. As shown in 
Fig. 4, air for combustion of a process material is supplied to a fur- 
nace at a pressure of approximately 200 in. water gage. The 
products of combustion pass from the furnace through a waste- 
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heat boiler, where heat is recovered by generation of steam, and 
thence to other equipment where the remainder of the process is 
completed. Steam from the waste-heat boiler is expanded 
through a back-pressure turbine to drive a blower which supplies 
air for the combustion process. 

Elsewhere in the factory, steam is required for space heating, 
water heating, vacuum jets, and miscellaneous services. In the 
power house, pulverized coal-fired boilers supplement the steam 
produced in the process waste-heat boiler. Plant electrical re- 
quirements are supplied partly from on-site turbine-driven gen- 
erators and also from the local public utility. About half of the 
on-site electric-generating capacity is from back-pressure or 
“topping” units and the remainder is from condensing units 
which receive steam exhausted from the topping units. The 
varying steam demand for feedwater heating is handled by proper 
selection and loading of steam and electric-driven auxiliaries. 

In designing the waste-heat boiler, one of the important de- 
cisions to be made is that of establishing the optimum pressure 
levels at the superheater outlet and at the blower-turbine exhaust. 
Another problem is involved in the selection of a turbine for driv- 
ing the blower, with regard to whether it should be a single-stage 
unit, or a multistage unit with high efficiency which would entail 
a substantial increase in cost. 
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STEAM FLOW RATE 
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Fig. 5 Steam generation by waste-heat boiler and steam required by 
blower turbine 


These relations can be visualized easily by reference to Fig. 5, 
which shows steam generation and steam requirements for the 
blower turbine at various production rates, and under the pro- 
posed operating conditions. 

Steam generation in the waste-heat boiler is proportional to 
production rate, irrespective of boiler pressure, since the waste- 
heat boiler is actually a process cooler in which heat removal is 
regulated to match production requirements. Power required by 
the blower varies as the cube of the production rate, since no 
throttling of air is used and the blower speed is regulated to match 
the system resistance. After allowing for expected turbine per- 
formance at varying speed and horsepower, and under the pro- 
posed operating conditions, the over-all relations are as shown. 

Production is limited at the point where the curve of steam 
generation intersects the curve of turbine steam consumption. 
The lower limit is computed for a single-stage turbine, with steam 
conditions of 400 psig, 600 F, at the throttle and an exhaust pres- 
sure of 150 psig. The next limit is computed for the same turbine 
with steam conditions of 500 psig, 700 F, to 140 psig. Thetwo 
highest limits are for a more efficient turbine operating under the 
same two sets of steam conditions. 

These relations show an advantage in production capacity for 
the more efficient turbine and the more favorable steam condi- 
tions. To evaluate this advantage fully, however, a number of 
factors need to be considered. Aside from the consideration of 
productive capability, there is also the matter of operating cost 
when running at less than full capacity since this represents a more 
normal production level. This involves consideration of certain 
factors in the power-plant operation, as well as in the process area. 

The relations between the process waste-heat boiler, the blower 
turbine, and the power plant can be visualized readily by refer- 
ting again to Fig. 4. It will be seen that the process waste-heat 
boiler and the power-plant boilers operate in parallel, supplying 
steam to a common system at approximately 500 psig and 700 F. 
This system supplies throttle steam to the blower turbine and also 
the topping turbines for electric generators, all of which exhaust 
at an intermediate pressure of approximately 140 psig. Certain 
boiler auxiliary drives also receive steam from the 500-psig sys- 
tem, but exhaust at 5 psig for feedwater heating. The 140-psig 
system supplies steam to process, to condensing turbines for addi- 
tional electric power generation, and to other power-house aux- 
iliaries which exhaust to stage feedwater heaters. 

The electric generators supply only a portion of the plant re- 
quirements, depending on steam-balance conditions, and the re- 
mainder is supplied by the local public utility. 

Under most operating conditions more steam is generated in the 
waste-heat boiler than is required by the blower turbine, with the 
surplus being available for electric power generation in the top- 
ping units. The coal-fired power-house boilers supply any addi- 
tional steam needed for process use and also for generation of extra 
electric power if and when it is profitable to do so. 

Evaluation of Operating Factors. To evaluate thoroughly the effects 


Journal of Engineering for Industry 


of different levels of steam pressure for the process waste-heat 
boiler or of blower-turbine efficiency, it is necessary to consider 
process production rates, process steam requirements, fuel cost, 
purchased electricity cost, and numerous details of less signifi- 
cance. With the results of such evaluations at hand and with 
estimates of investment cost for the alternatives being considered, 
it is possible to select the over-all optimum arrangement. 

The foregoing discussion has dealt with some of the factors in- 
volved in design and selection of process and power facilities for 
typical industrial plants. It seems appropriate at this point to 
add a few remarks about control of operating costs, once the facili- 
ties have been installed. 


Control of Operating Costs 


Whereas the over-all performance of an electric utility plant can 
be computed with remarkable accuracy and expressed in simple 
terms, such as net heat rate, the problem on the average industrial _ 
power plant is much more complicated. There are many good 
systems in use for measuring power-plant performance and con- 
trolling operating cost, and the one we use in our company has 
been found very effective. The objectives of our system are: 


1 To provide management with a measure of over-all perform- 
ance so that any change for better or worse can be detected 
easily and quickly. 

2 Tospotlight preventable losses and evaluate them. 

3 To show the causes of preventable losses and thus aid in 
correcting them. 


Most industrial power plants export several commodities, such 
as steam, electricity, refrigeration, compressed air, and water. 
Often steam is supplied at several pressure levels, refrigeration at 
several temperature levels, and water at several levels of quality. 
The performance standards used on du Pont power plants provide 
a measure of over-all efficiency, which is expressed as the ratio of 
standard total dollar cost to actual total dollar cost. The items 
included in each factor of ‘total dollar cost’? are fuel and any 
utilities, such as electricity, which may be purchased from outside 
sources. Operating labor, maintenance labor and material, mis- 
cellaneous supplies, and fixed charges are handled separately. 

Standard total dollar cost represents the minimum combined 
cost for fuel and purchased utilities which should be incurred for 
a given set of exported loads and uncontrollable operating condi- 
tions. In setting up the standard, allowance is made for scheduled 
maintenance which occasionally takes the most efficient units out 
of service. 

The standards are based on the best performance attainable, 
with each piece of equipment in first-class condition. In many 
cases, this corresponds with the performance observed during ac- 
ceptance tests. In other cases, the manufacturer’s predicted per- 
formance is used. In still other cases, standard performance is 
computed using accepted engineering procedures. 

After establishing standard performance data for each piece of 
equipment through the range of operating conditions to be en- 
countered, the information is put together in a manner which 
facilitates the calculation of optimum loading of all equipment in 
the entire power plant. Then, starting with actual exported 
loads and other uncontrollable factors such as river-water tem- 
perature, the results engineer computes the standard or minimum 
over-all cost for the period of time under consideration. This 
standard cost is compared with actual cost for the same period 
and serves as a measure of over-all performance. 

Location and evaluation of individual losses require detailed 
information regarding actual loading and performance of indi- 
vidual pieces of equipment. Care must be exercised, of course, in 
assigning properly the causes of loss. For instance, excess sieam 
consumption of a turbine-driven induced-draft fan may be due to 
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excess boiler load, improper fuel-air ratio for combustion, excess 
draft loss, leaks in the air heater, low fan efficiency, low turbine 
efficiency, or nonstandard conditions of steam pressure and tem- 
perature at the turbine throttle or exhaust. 

In a plant such as the one with the waste-heat boiler, described 
previously, it is usually necessary to give more attention, from 
day to day, to equipment loading than to the actual performance 
of individual units. 

Referring again to Fig. 4, it is more economical for this plant to 
generate than to purchase electricity under the following condi- 
tions: 


1 Whenever there is a surplus of steam from the waste-heat 
boiler which otherwise would be wasted. 

2 When the coal-fired boilers are needed for process use but are 
able to pick up additional load; provided, however, that: 

(a) Steam can be passed through both the topping turbines 
and the condensing turbines, picking up electric load in each. 

(b) Load carried on the condensing units is greater than the 
“break-even’’ point as compared with purchased power cost. 


Aside from considerations within the power plant itself, certain 
operating conditions in the manufacturing process can be altered 
so as to affect the generation of electricity in the power plant. For 
instance, some process changes will alter the steam required by the 
waste-heat-boiler blower turbine. These changes are accom- 
panied by changes in process manufacturing cost and they also 
affect the amount of steam available for use in the turbine genera- 
tors. In general, an advantage in one area is accompanied by a 
penalty in the other area. All of these relations are held under 
scrutiny continuously, and changes are made as required, to hold 
the over-all manufacturing cost at a minimum. 

Automatic Data Systems. Applications of automatic data logging 
and monitoring systems, together with electronic computers and 
controllers, hold great promise in operations of this sort. Some of 
these components are being used now with success on certain du 
Pont chemical processes, although they have not yet been applied 
to operations as involved and interrelated as those in the foregoing 
example. The speed and extent to which automatic systems of 
this type will replace human operators depend largely on econom- 
ics. Here again, the principal consideration will be safety, re- 
liability, and return on investment. 


Conclusion 


The foregoing discussion deals with problems which are typical 
of those encountered in the design and operation of power facilities 
for present-day industrial plants. In some cases, decisions regard- 
ing design of such facilities are based on judgment and experience, 
with a minimum amount of economic analysis. On the other 
hand, it is usually worth while to evaluate thoroughly all angles 
which have a bearing on over-all investment and operating cost. 
Even a well-designed plant, however, can be operated poorly and 
can experience large preventable losses. The system of power- 
plant performance standards, discussed briefly in this paper, is a 
valuable tool which can be used effectively in holding such losses 
to a minimum. 


DISCUSSION 
D. C. Calhoun? 


The author has covered his subject very well. He, like the 
writer, is engaged in heat balances in a chemical plant. Those 
who are familiar with this type of utility engineering will agree 
that it is most challenging. The utilization of low-level heat from 


1 Union Carbide Chemicals Company, South Charleston, W. Va. 
Mem. ASME. 
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chemical units and the balancing of the low-pressure steam sys- 
tem appear to be two of the most important problems outside 
of maintaining a constant supply of adequate heat. 

No one will deny that judgment and experience are important 
factors in determining what steps should be taken to improve the 
over-all heat balance of a plant. After a system is carefully de- 
signed and installed it may be found that problems develop 
which were not taken into consideration and, to a greater or lesser 
degree, will change the economical picture. 

Such a case which was experienced at first hand will serve as an 
example: 

One of our boiler plants has almost 100 per cent make-up to 
the feedwater. Since the feedwater comes from the river through 
a treating plant and since the steam plant produces steam for 
process in a chemical plant, there are many opportunities for the 
recovery of low level] heat into the feedwater. A low level heat 
recovery system was carefully engineered and finally installed. 
A long pipeline was required to carry the feedwater to and from 
the point of heat pickup. After considerable operating time it 
was found that the liberated oxygen in the warm return line was 
attacking the cast-iron pipe and causing a build-up of iron in the 
boilers. To prevent this condition from continuing, threshold 
phosphate treatment of the boiler water was inaugurated only to 
find that it not only coated the pipeline but also the heat-transfer 
surfaces in the high-level tubular feedwater heaters in the boiler 
plant. It became necessary to clean these heaters chemically at 
regular intervals to maintain their efficiencies. Thus we see that 
unforeseen difficulties can arise and materially affect the economy 
of even a carefully engineered system. 

The writer firmly believes that the basic data development 
should be completed prior to request for authorization on any 
project. Actual design, procurement, and construction can then 
be carried on simultaneously. Some changes in basic design may 
be necessary to take advantage of equipment design known only 
as procurement progresses. Field-construction changes are ex- 
pensive and should be kept to a minimum. 

The author’s two examples of specific problem conditions are 
good and only add impetus to the preparation and constant re- 
view of the internal steam and power balance. 

In most industrial areas, electric power can be purchased from 
a public utility at a lower average rate than can be generated ona 
condensing cycle by the industrial power plant. Most industrial 
plants require steam at various pressure levels and rates of con- 
sumption. Mechanical power, when produced by steam turbines 
under such conditions, will usually prove to be more economical 
than when produced by electric motors. It should, how*v<r, be 
remembered that the piping for an adequate steam distribution 
system may prove to be both involved and expensive. The in- 
creased use of high-speed centrifugal compressors has made the 
use of steam turbines a definite advantage over motor drives. 
The ideal situation is to have the steam turbine in the area re- 
quiring sufficient quantities of low-pressure steam, and thus only a 
high-pressure steam main need be installed to the turbine with 
appropriate low-pressure steam lines leading from the turbine to 
the nearby low-pressure steam load. 

Since industrial plants often produced a variety of products 
and since plant sites are usually not abandoned because of the 
failure of any particular product manufactured at that plant, it is 
reasonable to say that the utilities system can be modified to per- 
mit its effective use for other products and should therefore not 
be required to show as high return on investment as is required of 
other manufactured products. 

Major consideration in the design of industrial utilities involves 
(a) maximum continuous utilization, (b) protection against un- 
scheduled outages, (c) minimum operating labor consistent with 
safe practice, and (d) good steam balance with minimum vent to 
atmosphere. 
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The author stresses the fact that an economic evaluation must 
be made each time a major plant project is considered. The 
steam and power balance must be reviewed constantly to de- 
termine when, and if, utilities changes will result in better uti- 
lization of existing facilities and to predict the need for changes 
or additions to the system. 


L. F. Deming? 


The author is to be commended for bringing to our attention 
some of the many considerations which determine the design and 
economy of the plants which provide part or all of the utility 
services required by industrial processes. 

The writer is impelled to mention certain considerations not 
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cited by the author which generally affect the planning for utility 
services and upon which the author may desire to comment fur- 
ther. It is presumed that the basic data development has in- 
cluded an evaluation of the cost and availability of raw materials 
and services required before selection of a site. Since fuel is one 
of the raw materials generally accounting for one half to two 
thirds of the cost of the generated utility services, it would be in- 
teresting to know if, in the selection of fuel, the trends in fuel cost 
are projected into the future so that planning comprehends the 
use of the most economical fuel during the service life of the plant. 
If so, the basis of that projection would be of interest. An al- 
ternate solution, and the one generally followed by the Navy in 
areas where more than one fuel is either now or potentially avail- 
able, is to design the plant so that it is either equipped or readily 
adaptable to the economical use of an alternate fuel. 

Utility services must of course measure up to the requirements 
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of the plant served, both in capability and integrity of service. 
It has been found that by analysis of load components certain 
loads can be served in off-peak periods. In addition, base loading 
of the most efficient equipment can effect substantial economies. 
The Navy has found that by a Utilities Conservation Program 
substantial savings in utilities cost can be effected. 

The steam-load duration curve, Fig. 6, for optimum loading at 
the U. 8. Naval Gun Factory, Washington, D. C., for fiscal year 
1954, illustrates a study of steam loads. By restricting the waste- 
ful use of steam and by relating the loads shown to the boiler 
capabilities and loading the most economical units to their 
capacity, the peak loads can be carried by the older and less ef- 
ficient units. Under such conditions the short duty cycle mini- 
mizes the effect of the lower efficiency of the older units. A 
similar graph, Fig. 7, shows the utilization of the electrical plant 
capacity under the conditions of steam loading previously shown. 
Notice that the No. 6 unit operated throughout the year and at 
rated capacity for about 4000 hr per year. The less efficient No. 5 
unit, on the other hand, operated at capacity slightly over 2000 
hr or about 8 hr per working day. Purchased power is limited in 
demand and to about 700-hr utilization per year to provide the 
most economical combination of purchased power, electric power 
produced by extraction generation to the maximum extent possi- 
ble with the export steam load, and electric power generated by 
condensing operation. The analysis only is typical since the re- 
sults of course change with fuel costs, power rates, and so on. 

The conditions of urgency which generate the programming and 
scheduling illustrated by the author’s Fig. 1 are serious indeed to 
the private owner and not contributory to optimum economy. 
The Navy meets this emergency by utilization of equipment in its 
mobile power reserve which it owns and maintains to meet 
strategic requirements. 

In conclusion, it is desired to commend the author for inviting 
our attention to the multitude of factors which affect the com- 
parative economy both of operation and of capital invested for 
the industrial plant. 


Author’s Closure 


The author wishes to express his appreciation for the written 
discussions by Messrs. Calhoun and Deming and the oral com- 
ments given at the Annual Meeting. 
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As pointed out by Mr. Calhoun, the economic picture often 
changes from the original concept of a project, thus emphasizing 
the need for flexibility in design. It is, of course, desirable to 
complete basic data development prior to final design of facili- 
ties, but compromises often are justifiable or even necessary for 
expediency. Attractive profits obtainable from a new product 
often shrink significantly after competition enters the field. 

Back-pressure turbine drives are used in many of our process 
applications, located near the areas which can utilize the exhaust 
steam. The economics are affected by many factors, including 
local or state regulations which often dictate that prime movers 
be attended by licensed operators. No such requirements apply 


to electric motor drives. 


Mr. Deming raises the question as to whether trends in fuel 
cost are projected into the future as a basis for design of fuel 
burning facilities and planning operating economy during the 
life of our industrial plants. Such a procedure is a very desirable 
goal which is only partially attained. Trends in fuel prices often 
are upset by unpredictable factors, such as lobbying by various 
private interests and political maneuvering in Washington. 
Most of our prognostications are based on generally known trends, 
tailored to fit local conditions. In marginal cases, where coal 
and oil or gas are competitive, we sometimes design for oil or 
gas but leave space for later installation of coal handling, prep- 
aration, and burning facilities. A number of our plants are 
equipped to burn either type of fuel and some are equipped to 
burn two or three fuels simultaneously. In the operation of all 
our plants we follow what we have termed the “utility” values of 
various available fuels in order to obtain the lowest over-all 
operating cost. This involves many factors and requires con- 
tinual surveillance, but the rewards are large. Further discussion 
of this subject is given in the following ASME papers: 

L. F. Deming, ‘Navy Practice of Designing for Dual Fuel 
Firing in Shore Station Steam Plants,” Paper No. 55—Fu-4. 

W. H. Attwill, C. E. Day, and A. J. Johnson, “Evaluating Coal 
by Utilization Cost,”’ Paper No. 56—A-211. 

R. M. Schahfer, F. H. Strout, and G. W. Land, ‘Application of 
Bonus-Penalty Method to Coals Being Purchased or Evaluated,” 
Paper No. 57—A-115, 

E. D. Holdup, “Sound Methods of Solid Fuel Evaluation for 
Use in Thermal Power Stations,” Paper No. 57—Fu-]. 
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Influence of Hot and Cold Storage- 
Loads on the Refrigerated-Space 


CARL F. KAYAN 


Professor, Department of 
Mechanical Engineering, 
Columbia University, 

New York, N. Y. Mem. ASME 


Temperature of a Cooling Complex 


Of particular interest in cold-storage warehouse operations 1s the influence on the space- 
temperature-time history of a “‘hot’’ load of goods suddenly introduced into the refrig- 
erated zone of an insulated structure. 
occupied by warehouse goods already down in temperature. 


The zone may be completely empty, or partially 
Of further interest is the 


ultimate effect of the transient load on the performance of other components of the system, 
; such as the connected space-cooling heat exchanger. In previous work in this realm 
by the author on the behavior of the over-all complex, the studies have included the various | 


heat-exchange and machine energy-transport equipments of the system [1]}.! 


Whereas 


in the earlier paper the transient effect was explored via calculation methods, in the 
study of the problem in this paper experimental simulation methods are employed 
entirely, using pilot process-prototype transient-elements. 


Is THE present study of cold-storage warehouse 
operations, attention is focused on the behavior of the structure, 
the space, and the loads under given temperature conditions, as 
coupled with the space heat-exchanger equipment. Warm goods 
brought into an empty refrigerated space represents one type of 
transient cooling problem. However, the presence of existing 
stored goods at refrigerated temperature adds to the problem of 
predicting net heat-flow rates and heat-exchanger performance. 
Finally, the behavior of the structure wall, considered with real 
heat capacity versus one with negligible heat capacity, further 
complicates the analysis. Once again, orthodox calculation pro- 


1 Numbers in brackets designate References at end of paper. 

Contributed by the Process Industries Division and presented at 
the Annual Meeting, New York, N. Y., November 30—December 5, 
1958, of THe AMERICAN SocrETy OF MECHANICAL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, Sep- 
tember 12, 1958. Paper No. 58—A-165. 


cedures prove inadequate, and analogy methods using rheoelec- 
tric simulation techniques with pilot prototypes prove worth 
while. 

Using specialized electrical-potential transfer procedures along 
with resistance-capacitance representation for the different heat- 
capacity components, the different rates of energy flow are fol- 
lowed time-wise. The net energy flow is sensed in the heat-ex- 
changer simulation circuit; this is set up for a given value of 
brine coolant and circulating room air flow and brine-supply tem- 
perature. For the integrated system, resultant space and goods 
temperature, heat-exchanger energy-flow rates, as well as fluid 
temperatures, are predicted as a function of elasped time for the 
different cases covered in the program of experiments. 

Studies via rheoelectric analogy, as extensions of the re- 
sistance concept of heat transfer, have covered various applica- 
tions; in heat conduction, both the geometrical-sheet and net- 
work methods have been used by the author [2, 3]. Also to be 
mentioned are heat-exchanger analysis [4], energy-transport [5, 
6], and pipe-network fluid-flow distribution [7]. Rheoelectrical 


Nomenclature 
A = surface area, sq ft q = heat flow rate, Btu/hr t, = surrounding air temperature, deg 
a = k/(pc) = thermal diffusivity R’' thermal resistance, referred to 1 sq F 
c = specific heat, Btu/(lb-deg F) ft, deg F/[Btu/(hr-sq ft)] t,, = external atmospheric temperature, 
C, = specific heat of air, Btu/(Ib-deg F) R,’ thermal resistance of surface air, deg F 
¢, = specific heat of brine, Btu/(Ib-deg deg F/[(Btu/(hr-sq ft)] t) = original initial solid temperature, 
F) R thermal resistance, deg F/(Btu/hr) deg F 
C, = electrical capacitance (farad, mi- R, = heat-exchange resistance of air ts = surface temperature, deg F 
crofarad ) stream, = 1/(W,c,) t, = temperature in solid at any point 
C,’ = thermal capacity per unit of R, heat-exchange resistance of brine x, deg F 
frontal area, = cpY stream, = 1/(W,¢,) At = temperature difference, deg F 
fs = Schack heat-transfer function, Ry, = exchanger heat-transfer resistance, W, = weight flow rate for air, lb/hr 
temperature-difference ratio = 1/(UA) W, = weight flow rate for brine, lb/hr 
E = electrical voltage for analog R,' thermal resistance of solid path Y U = exchanger over-all coefficient of 
F = dimensionless — temperature-dif- in length, = Y/k, per sq ft heat transfer, Btu/(hr-sq ft-deg 
ference or resistance ratio R,’ = thermal resistance due to fluid F) 
h = heat-transfer coefficient, Btu/(hr- conductance h, per sq ft, = 1/h Y = half-slab thickness, ft 
sq ft-deg F) r, = electrical resistance, ohms, propor- p = density, lb/cu ft 
I = electrical current for analog tional to R,’ Tt = elapsed time, hours 
k = thermal conductivity of solid ma- =r, = electrical resistance, ohms, pro- Subscripts to denote special applications 


terial, Btu/(hr-sq ft-deg F/ft) 
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portional to R,’ 


are used where required. 
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analysis is carried along further in this paper, in the development 
of the transient refrigeration-load problem studied earlier [1]. 
In this work, pilot resistance-capacitance units are used to monitor 
the changing heat flows resulting from the varying temperature- 
difference between the solid-goods surface and the cooled-space 
air temperature. 

The system studied is shown in Fig. 1. Brine from a separate 
refrigeration plant is supplied at constant temperature. The heat 
load of the refrigerated space is absorbed in the movement of the 
space air, cooled through fan circulation in a countercurrent type 
of heat exchanger supplied with the cold brine. 

In the general development of the problem, many simplifying 
assumptions are found to be expedient: For example, the struc- 
ture (in one part of this study), the heat exchanger, and the heat- 
exchanger fluids are assumed to have negligible heat capacity; 
the heat-transfer coefficient for the body surface is assumed, and 
taken as constant in the range of operation; and the same is the 
case for the heat-exchanger fluid flows and their over-all coefficient 
of heat transfer. For the hot and cold transient loads, values of 
the thermophysical properties have been assumed, and are used 
as constant. The external air temperature also remains constant. 


Heat-Load Relationships 


Transient Load. For a solid slab-form body of given thermal 
properties, initially at a uniform temperature t throughout and 
then suddenly subjected to exposed-surface heat transfer with a 
surrounding fluid, the rate of heat exchange depends not only on 
the exposed area of surface and its heat-transfer coefficient, but 
also obviously on the realized temperature difference existent 
at any time between the surface and the fluid (such as the space 
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air in the present case). The time-varying surface temperature t, 


depends not only on its previous history and surface conditions ag _ 


well as on internal physical properties, but also on the value of the 
half-thickness dimension Y of the body. 

Thus referring to Fig. 2(a) to depict the solid section in ques- 
tion, in terms of its physical properties, and considering unit 
frontal surface area (1 sq ft) exposed to the action of the air co- 
efficient h,; with R,’ representing the thermal resistance of the 
half section Y in depth, R,’ representing the surface resistance 
due to coefficient h,, and with a as the thermal diffusivity: 


For sustained temperature difference between initial and uni- 
form solid temperature f and the surrounding air ¢,, the dimen- 
sionless ratio F, or specifically F, for the temperature-difference 


relationship at any point x for the solid (from surface to center- 


line core), is important: 


For the surface 


(4) 


In a previous study [1], reference was made to the “‘transient- 
body”’ resistance Ry’. This may be defined as an equivalent 
time-wise varying resistance which, when divided into the over- 
all temperature difference (t) — ¢,) determines the rate of heat 
flow per sq ft, thus q,, for the surface-area A,: 


— ta 
qu = A, (“=*) (5) 
Specifically 
R,’ = =— 6 
tb (6) 


The function f; is the specialized one for this type of case, cover- 
ing the transient surface temperature, defined by Schack [8] and 
originally used in the calculation of the transient-body resistance 
{1]. It represents the surface temperature-difference ratio ef- 
fect, that is, (t, — t,)/(to — t,), which herein is also defined as F,. 
F, may be used either as a calculated value, or as one otherwise 
established such as shown subsequently, through electrical simu- 
lation technique. 


It is of interest to note a rearrangement from Equation (6): 


The heat-flow g, between a transient-body load and the sur- 
rounding fluid such as air, may, of course, be stated in terms of 


the corresponding temperature difference and the applicable re- 
sistance for the load: 


t, — tg 3 1 

= ( R, ), with Ry = (8) 
Structure. With fresh ventilation-air (infiltration) cooling ef- 
fect included, for steady state, the enclosure hourly heat-flow 
rate g,, in accordance with the resistance concept of heat transfer, 
may be represented by an equivalent wall resistance R,, and the 
temperature difference At,, between the external outdoor and the 

internal space-air temperature (which may be varying): 


Vo = At,/Ry (9) 
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Ilnjer steady-state conditions, the heat capacity of the wall 
~~qoes not enter into the heat-transfer effect. However, should the 
face temperatures of the wall be altered by transient heat-flow 
effects, then the storage or heat-capacity property of the solid 
material comes into play, in a somewhat similar manner to that 
involved in the slab-form body discussed previously, under which 
circumstances the wall heat flow may be greater or less than that 
for steady state. Such a case is treated in the experimental pro- 
gram to be discussed subsequently. 


Heat-Exchanger Relationships 


The air-to-brine heat exchanger may be considered on a coun- 
tercurrent basis. As developed in earlier papers [1, 4], its per- 
formance may be explored and predicted through an extension of 
the resistance concept. Here the exchanger heat-flow gq, is re- 
lated to the over-all exchanger temperature difference At, , be- 
tween the entering air and the entering brine (f, — ¢,), through 
the total over-all thermal resistance R,, ,: 


Ru. = | ite — 


(10) 


(11) 


with R, and R, representing the heat-exchange resistance of the 
air and brine streams, and Ry, representing the over-all heat- 
transfer resistance. 

As will be apparent in the specific example to be cited, propor- 
tional distance along the total resistance path R,, , will represent a 
measure of the temperature change for the different fluids work- 
ing in the exchanger. 


Rheoelectrical Simulation Analogy 


Inasmuch as the solid body undergoing transient change com- 
prises both thermal resistance and thermal storage (heat-capacity 
effects), its internal thermal behavior may be simulated elec- 
trically through electrical resistance and electrical capacity. 
For the structure, as noted before, under steady-state conditions, 
only resistance comes into play, but with transient condi- 
tions, heat capacity. (Hence under the latter conditions, elec- 
trically, here likewise both electrical resistance and capacitance 
come into play, as for solid loads, and represented on a proto- 
type basis.) 

Fig. 2(b) thus represents the counterpart of Fig. 2(a) with re- 
spect to a solid slab, with the surface-to-center solid section sub- 
divided into a number of slices. Each slice mid-point is repre- 
sented by its condenser, and slice-to-slice connections are intro- 
duced via appropriate network resistances, as conventionally 
done in this field [9]. Similar provision may be made for the 
surface heat-transfer effect. 

As a first step in the present approach, the electrical values 
may be considered, individually, directly proportional to the 
thermal values, the solid body here being taken as the half 
thickness of a slab of infinite extent (edge effects neglected). 
Thus Equations (3) and (4) have their electrical counterpart, 
with temperatures replaced by their voltage counterparts. It is 
to be noted, therefore, that in the simulation process the dimen- 
sionless ratio F may be determined electrically: 


E, E, 


AE), 
E, = E, + Fo — E,) 


(12) 


Thus in effect 
(13) 
Stated otherwise 
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The dimensionless ratio F thus shows as a function of the specified 


It is significant that, since the electrical (rC) value is involved 
as a product, the result is independent of the value of the separate 
factors, and is so generally regarded. Therefore, for electrical 
simulation, the setting up of the (rC) product can be accomplished 
in terms of convenient values and chosen quality of condenser 
capacitance C,, with the resistance r, appropriately evaluated. 
Thus in terms of the value for r, so established, and through 
the required value for the ratio r,/r,, as established thermally, r, 
may be valued appropriately. In effect, this characteristic also 
permits a convenient choice of values for a desirable time relation- 
ship for the “pilot” resistance-capacitance unit, used (voltage- 
wise) to monitor the main circuit as a segregated prototype. 

The behavior of this pilot unit is then one of producing a voltage 
value for the “surface temperature,” which may then be applied 
separately to the monitoring of the main-circuit load-resistance 
component. 

The application of this transient-body simulation approach 
will be demonstrated best through the subsequent illustrative ex- 
ample. Here, in addition to the transient problem for storage 
loads and in part also for structure, the concurrent behavior of the 
heat exchanger is incorporated into the system. 


(14) 


(15) 
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Fig. 3 represents the general rheoelectric simulation circuit for 
the study of the subsequently cited examples. In addition to the 
resistive components for the structuré (in one part) and for the 
heat exchanger, load-transfer slide wires also are included. Via 
the pilot units, surface-temperature voltages result from the 
operation of the (rC) pilot networks; their value is reproduced at 
the load slide wires, so that in this way it is possible to track the 
real load change for the main circuit (as a function of time). 
Thus the flow of heat off the exposed surface of the transient 
loads the surrounding air may be simulated. 


lilustrative Example 


Time-wise performance prediction is desired for the refrigerated 
space and heat-exchanger system, represented by Fig. 1, under the 
following conditions, general: External temperature, 100 F, and 
brine-supply temperature (constant), 12 F. Two series of ex- 
periments are reported, Series I, using a wall of negligible heat 
capacity, and Series II, a concrete-cork wall, with real capacity. 

Series I: Design space temperature, at steady state, 32 F: (a) 
Load X, hot, plus load Y, hot (i.e., both uniformly and initially 
at 100 F), together charged into the empty 32-F space. (b) Load 
X already established in 32-F space, and uniformly at this tem- 
perature; load Y, hot, uniformly and initially at 100 F, suddenly 
introduced into the partially occupied space. 

Series II: Design space temperature, at steady state, 18 F: (a) 
Hot-wall condition, uniformly at 100 F, and cooling effect with 
12-F brine on heat exchanger, suddenly started. (b) Concrete- 
cork wall at steady-state conditions, 18-F space, and empty, then 
with sudden introduction here of load X with 625 sq ft of exposed 
surface, hot at 100 F initially. 

For these cases the temperature-time history, as well as heat 
flow-time history is desired. 

Thermal-Load Data. Structure heat load (ventilation-infiltration 
air included) for Series I wall of negligible heat capacity = 1260 
Btu/(hr-deg F); for Series II, 8-in. (0.667 ft) concrete outer wall, 
with 6-in. (0.50 ft) cork—presumed to represent the entire 
envelope of 7200 sq ft. Heat-exchanger: Brine flow rate W, = 
8000 lb/hr, c¢, = 0.70 Btu/(lb-deg F); air flow W, = 38,500 
lb/hr, c, = 0.241 Btu/(lb-deg F); UA = 11,850 Btu/(hr-deg F). 

Load X: Tightly packed vertical stack of goods, 2 ft thick, of 
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500 sq ft total exposed area; goods thermal conductivity k = 
0.70 Btu/(hr-sq ft-deg F/ft); density p = 94 lb/cu ft; ¢ = 0.22 
Btu/(lb-deg F); equivalent thermal diffusivity, a = 0.034. Load 
Y: Tightly packed vertical stack of goods, 2 ft thick, of 625 sq ft 
total exposed area; thermal conductivity, k = 0.35; density, p 
= 94 lb/cu ft, c = 0.22; thermal diffusivity, a = 0.017. The 
characteristics of load Y were so chosen, to have a thermal 
diffusivity of one half of load X. For both load X and load Y, 
the surface conductance h, is assumed to have the general value 
of 2 Btu/(hr-sq ft-deg F). Concrete: k= 0.50; ¢ = 0.27; p = 133; 
a = 0.0139. Cork: k = 0.025; c = 0.50; p = 10; a = 0.005; 
h, = 1.65 as a general value for both exposed sides of the wall. 

It is to be noted that both load X and load Y, under the dif- 
ferent programs of operation, will be subjected to heat transfer 
under time-wise variable space temperature. Accordingly, for 
reference purposes, the general characteristics of temperature 
variation, for both X and Y, under sustained heat-transfer fluid 
conditions, will be of interest. These characteristics have been 
determined separately through electric-analogy, resistance- 
capacitance studies of the ‘pilot loads.’’ 

Figs. 4 and 5 refer to load X, with its defined thermal character- 
istics. Fig. 4 shows the measured variation of the dimensionless 
ratio F for both the exposed surface and the core, under the de- 
fined surface-conductance condition versus elapsed time, and 
under constant-level ‘‘cooling’’ conditions. Fig. 5, for the same 
load, cooling, shows the variation of F versus internal position, 
for various values of elapsed time in hours. 

Figs. 6 and 7 refer toload Y. Fig. 6 is the counterpart of Fig. 4 
which applied to load X, and Fig. 7 is the counterpart of the load 
X, Fig. 5. 

For these studies, and for the subsequent ones on over-all per- 
formance of the space system under load conditions, a time ratio 
of 60:1 was established for the thermal versus electrical behavior. 
Thus in effect, for load X, for a value of 7/(R,'C,') = at/Y? = 
1.00 = 0.0347, giving r a value of 29.4 hr, this represents 29.4 
min = 1760 sec electrical. Similarly for load Y, 58.8 br thermal 
was represented by 58.8 min = 3520 sec electrical. Electrical 
circuit values were translated from the thermal numbers on a 
proportional basis throughout, for voltage versus temperature, 
for current versus heat-flow rate, and for ohms electrical re- 
sistance versus thermal resistance. 
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In the voltage ‘transfer’ technique used, great flexibility is 
possible since, in the pilot procedure, only the voltage pattern 
needs to be obtained. Thereby the voltage values at the load 
slide wires are monitored to govern the current flow which passes 
through the load surface resistance to the space-voltage bus. As 
stated before, for the pilot system, the electrical value of r, 
was set up via r, of the (r,Ce) value, in accordance with the re- 
quired thermal ratio of R,’/R,’. 

Performance Predictions. Series I: Figs. 8 and 9 cover the stipu- 
lated Case (a) of both load X and load Y, initially and uniformly 
hot at 100 F, suddenly moved into the refrigerated space. Fig. 8 
shows the temperature-time history of the two loads, reported as 
to core and surface temperature, the heat-exchanger fluid tem- 
peratures (brine in at 12 F constant), and lastly but of great 
significance, the variation of the space temperature. The effect 
of suddenly charging the cold zone with masses of hot goods is 
sharply apparent. Fig. 9 reports various energy flow rates, for 
the structure, the two loads X and Y, and the ultimate heat ex- 
changer. 

Figs. 10 and 11 report the results for Case (b), with load Y, hot, 
suddenly charged into the 32-F space already occupied by load X 
at its uniform temperature of 32 F. Fig. 9 shows the tempera- 
ture-time history with its different detail, with special attention 
due the cold load X as to its behavior. Fig. 11 again reports the 
different energy flows. 

The Series II results concern the system with a wall structure 
of finite heat capacity. For this arrangement, a general study of 
the wall characteristics, to parallel the load studies depicted in 
Figs. 4 and 5, for load X, and in Figs. 6 and 7, for load Y, presents 
some interesting results. Fig. 12 shows the value of F as a func- 
tion of elapsed time, for different position points within the wall, 

ncluding the concrete-cork interface. Fig. 13 shows F versus 
distance inward from the outer concrete-wall face, reporting 
values for different elapsed-time periods. 
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Case (a), Series II, concerns the starting-up operation for the 
concrete-cork structure, initially uniformly at the 100 F tempera- 
ture. Fig. 14 shows the temperature-time history, emphasizing 
the sizable time required to approach the steady-state conditions 
of 18-F design-space temperature. Fig. 15 shows the tremendous 
variation of the heat load on the heat exchanger under the start-up 
conditions, with the brine temperature fixed at the 12-F value. 

Figs. 16 and 17 cover Case (b): The empty structure, at steady- 
state condition (18 F) suddenly charged with a hot load of X com- 
position, of larger size than usual, to wit, of 625 sq ft exposed area. 
This entails disturbance of the steady-state condition of the wall 
and of course involves an outright transient behavior on the part 
of the hot load. Severe disturbance in many places is to be 
noted. Fig. 16 presents the temperature-time history, and Fig. 17, 
the energy flows. 

It is again to be emphasized that, in addition to the use of many 
simplifying assumptions for reasons of expediency, the brine tem- 
perature throughout was maintained at 12 F, and the rates of 
flow for both the air and the brine remain unchanged. It is of 
course recognized that, in given plant operations, these conditions 
might well vary. Particularly to be noted is that the brine-supply 
temperature, often dependent on refrigerating-machine per- 
formance, might well rise due to the added heat load. Under 
these circumstances, the space temperature would be disturbed 
further, with extended effects on the different heat loads. . (fa 
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Vibratory Feeders and Unit Vibrators 


Highlights in the development of an entirely new electropermanent magnetic motor 


drive and a new disk-type springing system through the use of a laminated fiberglas- 
epoxy resin product, resulting in an outstanding line of vibratory feeders and unit or 
bin vibrators. Also some of the many uses for these products in controlled materials 
handling, as well as engineering information which should be considered in the selection 
and application of these versatile performers. 


Mace HAS BEEN SAID and written regarding the sub- 
ject of controlling material flow by means of vibration. Many 
types of vibratory feeders and unit vibrators are available for 
this purpose and special skills have been developed in their design 
and application to problems of materials handling and control. 

Many engineers, and industry in general, are familiar with the 
principle of operation of vibratory feeders, which in essence con- 
sists merely of tossing the particles upward and forward on a 
feeder tray in a series of short controlled hops. The amplitude 
and frequency of the tray movement determine the rate of flow 
of the material. Because the amplitude of each movement is 
small and the frequency is so high, usually 60 cycles per second, 
casual observation of the material being conveyed could create 
the impression that something mysterious is going on. The 
material seems to flow forward or even up inclines or around 
curves as smoothly as a flow of liquid. 

There are many successful methods of accomplishing the cor- 
rect tray movement and there are arguments relative to the ad- 
vantages and disadvantages of each. No single method provides 
a satisfactory solution to the innumerable problems that are 
normally encountered in the field of vibratory feeding. 

Surveys indicate a definite trend toward selection of electro- 
magnetic feeders over mechanical, eccentric weight or air-operated 
types because of the simple, accurate, and versatile controls 
available. Multiple controls can be easily installed in central or 
remote locations. 

The type of vibratory feeder (Fig. 1) and unit vibrator (Fig. 2) 
with which this paper is concerned is a combination electromag- 
netic and permanent magnetic type. It performs the same func- 
tions as ovher vibratory products except that the design is such 
that many new features, including an entirely new type of motor 
drive, have been incorporated. These features contribute to a 
more stable, durable, and efficient product. 


Objectives 


At the outset it was recognized that development of a new 
product superior to proven electromagnetic feeders and vibrators 
would be difficult. Any work performed had to be preceded by a 
thorough study and analysis of existing products. This study 
was undertaken, and disclosed that existing products, although 
well liked and generaily accepted, would be more attractive from 


Contributed by the Process Industries Division and presented at 
the Semi-Annual Meeting, St. Louis, Mo., June 14-18, 1959, of THE 
AMERICAN SociETY OF MECHANICAL ENGINEERS. 

Norte: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, April 6, 
1959. Paper No. 59—SA-33. 
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Fig. 1 Typical small vibratory feeder with standard flat-pan-type tray 
for controlled feeding of bulk materials 


Fig.2 Typical unit vibrator for use on walls of kins or hoppers fo prevent 
bridging and to assure uniform flow 


the customer’s standpoint if the following features could be im- 
proved, corrected, or added: 

1 -A completely enclosed more powerful motor so that units 
could be safely and confidently installed in wet, dusty, corrosive, 
or sanitary areas. 

2 An enclosed motor that could be installed in hazardous, 
dusty locations where Class II, Group G products would ordi- 
narily be used. 
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3 Elimination of the need for a current rectifier or, in some 
cases, the entire control, and still operate at power-line frequency, 
usually 60 cycles. 

4 Less power consumption to achieve the same amount of 
work. 

5 Asimple method of tuning to avoid loosening or tightening 
bolts or adjusting magnet air gaps to some intangible degree, 
often requiring considerable skill and experience. 

6 A spring material that would not fatigue and eventually 
break. 

7 A spring material that would be noncorrosive and whose 
performance would remain constant over a reasonably wide 
range of temperature variations. 

8 A spring material that would possess inherent damping 
characteristics to provide more stable performance under varying 
load conditions. 

9 A simple design to eliminate existing stacks of steel springs 
and spacers. 

10 A way to prevent packing of foreign materials between the 
springs, which has sometimes been responsible for inconsistent 
feeder performance. 

11 Power-factor correction. 

12 Less front-end flip—excessive tossing or agitation of ma- 
terials at the discharge point of feeder trays. 

13 Eliminate rear-end dead spot—an area at the rear end of 
some feeder trays in which materials merely move about but will 
not move forward. 


Development 


With this imposing list of possible improvements as a target, a 
development program was instituted and carried out. It has 
resulted in the creation of a completely new line of smail vibra- 
tory feeders, through 10 tons per hour! capacity, and a new line of 
unit vibrators. Larger models of both lines are available but 
will not be covered in this paper. Through good fortune, per- 
severance, and ingenuity each separate objective has been at- 
tained. A complete story of the development of these two closely 
related product lines would not be possible within the limits of 
this paper; only the highlights and most important points can 
be covered. 


Fig. 3 Cutaway drawing of vibratory feeder showing enclosed electro- 
permanent magnetic motor and disk-type springs 


A New Type of Enclosed Motor 


The first noteworthy success was an entirely new type of motor 
drive (Fig. 3) which consists of a three-pole electromagnet inter- 
meshed with two poles of an Alnico V permanent magnet. The 
electromagnet is attached to the reaction mass and the permanent 
magnet becomes the armature and an integral part of the tray, 
or moving mass. This system provides a motor with alternate 
driving forces in each direction, the frequency of which coincides 


1 The accepted industry standard for capacity measurement is free 
flowing, dry sand, which weighs 100 lb per cu ft. 
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with the a-c power supply, without current rectification. It 
can be operated with no control. Where adjustment of material 
flow is necessary, a potentiometer or powerstat may be used. By 
driving in both directions and using both halves of each current 
cycle, the resulting motor is far more efficient than the common 
electromagnetic type. Several patent claims have been allowed 
on this new type motor. 

Simplification of the springing system and a simple method of 
enclosing the drive motor was much more difficult. However, 
the problem of enclosing the motor was finally solved by de- 
signing the reaction mass, or housing, in the form of a hollow cube 
which contains the fixed electromagnet and the armature. The 
front and rear ends of the armature were each attached to disk- 
type springs which were in turn bolted to round openings in the 
housing, resulting in a completely enclosed unit. The tray mass 
was then attached securely to the center bolt of the front disk- 
type spring and the rear of the tray was supported by a small 
leaf spring attached to the rear of the housing, the purpose of 
which will be explained later. 

This arrangement proved to be most ideal because it permitted 
the driving forces to be nearly in line with the centers of gravity 
of both the tray and reaction masses, which is prerequisite to a 
smoothly operating feeder. Furthermore, it also eliminated the 
need for large stacks of springs and spacers used on all other 
electromagnetic vibratory feeders, and the harmful packing of 
foreign materials which usually occurs between springs. 


Springs 

The selection of a spring material became increasingly imi- 
portant as it was recognized that the motor in this system is 
merely an exciter and provides a nominal amount of work com- 
pared with the energy supplied by the springs: 

An investigation and search for spring materials with the de- 
sired characteristics were carried on simultaneously with the 
other development activities. The search was intensive and 
broad and was punctuated by disappointments as various ma- 
terials were tested and discarded by development engineers. 
Among the materials tested were brass, beryllium copper, and 
special plywoods imported from Europe. From the very first 
analysis, a laminated fiberglas-epoxy resin product looked en- 
couraging and later exhaustive testing confirmed that it had all 
of the desired attributes. 

The most serious disadvantage in determining that laminated 
fiberglas was the best spring material was the fact that no sig- 
nificant technical information regarding its spring character- 
istics, especially as a disk-type spring, was available. So, 
working in close co-operation with the material supplier, facts 
and figures were formulated so that mathematical calculations 
and predictions could be made. 

Over a period of time, many development prototypes were 
made, observed, analyzed, and tested—sometimes to a point of 
destruction. Some were placed on continuous full-load life test 
and periodic observations were made. Others were deliberately 
dropped on concrete floors, without harm (except to the floors). 
One unit is now running completely submerged in a tank of water. 

The small leaf-type spring previously mentioned, at the rear 
of the feeder tray, not only provides support and stability to the 
tray but is also an effective tuning device. Changing to thicker 
or thinner tuning springs (to compensate for extreme but tolerable 
variations in work loads) is essential to sometimes prevent 
“killing’”’ the feeder and to avoid lowering the natural frequency 
of the system to coincide with that of the power supply. When 
this occurs, the entire unit could become resonant and would be 
self destructive (Fig. 4). Springs and masses are designed so that 
the natural frequency, or resonant point, of each system is from 
six to ten cycles greater than that of the line voltage. This fact, 
together with the inherent damping characteristic of the spring 
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Fig.4 Diagram showing effect of varying the load on a vibratory feeder 


material, results in a product that has stable performance over a 
wide range of working conditions. 

Standard feeders, as tuned at the factory, will perform effi- 
ciently on materials ranging from 50 to 150 lb per cu ft, but 
when maximum output of materials weighing less than 50 lb 
per cu ft is desired, a thinner tuning spring is used. If the product 
weighs more than 150 lb per cu ft or if heavier than standard 
trays are used, thicker or additional tuning springs must be added. 

The laminated fiberglas-epoxy resin spring material, unless 
overstressed, has practically no fatigue point and should last 
indefinitely. Extended life tests tend to bear this out. The 
material is highly resistant to most chemicals and acids and its 
performance does not vary appreciably under reasonably high or 
low-temperature conditions. 

Front-end flip and rear-end dead spot, previously defined, have 
been eliminated in this new design, although by lightening and 
extending the front end of a tray, front-end flip can be introduced 
in cases where agitation at the discharge point is desirable. 


Unit Vibrators 


Up to this point much of the text has been devoted to vibratory 
feeders. However, many of the statements apply as well to unit 
vibrators, whose primary function is to cause maximum safe 
vibratory movement of the walls of bins, hoppers, or chutes to 
assure the flow of bulk materials from the discharge opening at 
the bottom. Properly selected and correctly applied, they can 
reduce the likelihood of packing and arching of material flows. 

The same type of drive motors and spring materials are used in 
the new improved vibrators and all units are completely enclosed. 
Their design (Fig. 5) is such that the armature containing the 
permanent magnet is mounted in the center of a single leaf-type 
spring which is supported at both ends by the base housing. The 
hammer is also attached to the spring and in operation impacts 
either a steel or rubber-faced anvil in the base housing, usually at 
a rate of 3600 cycles per minute. In some of the larger models 
two back-to-back magnetic drives provide ample power to excite 
the spring and heavy hammer. 

Other inherent advantages common to both products are the 
facts that the magnet air gaps are fixed at the factory at their 
most efficient operating point and no further adjustment is nec- 
essary; also that power factor can be easily corrected by in- 
corporating a small capacitor in the control unit. 


Uses and Application of Small Vibratory Feeders 


There is actually no sharp dividing line between small and large 
vibratory feeders. However, feeders whose rated output is ten 
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tons per hour or less are usually classified as small feeders and are 
generally used in applications where space limitations prevail 


and where more precise control of the material flow is necessary. 


It should be emphasized that small vibratory feeders cannot be 
so precisely controlled that they can take the place of automatic 
weighing, packaging, or parts-counting machines. However, the 
degree of success of these automatic machines can usually be 
measured by the performance of the vibratory products which are 
a part of their make-up. 

The versatility of vibratory feeders is accented by the fact that, 
in addition to the most commonly used troughed-pan-type trays, 
they may be fitted with trays of other configurations. Other 
shapes commonly used are tubular, half-round, and V shape, 
although the many possible dimensional variations of these 
shapes make them adaptable to countless applications. 

Trays may have built-in heating elements for heating, drying, 
or maintaining heat in products being conveyed, or may even be 
water-jacketed for feeding materials such as glass batch to a fur- 
nace. They may be made of steel, aluminum, or stainless steel 
with special finishes or coatings such as teflon, glass, or epoxy 
resin when needed. 

The addition of one or more screens to trays is common practice 
for separation, grading, or de-dusting of materials. Still another 
variation is to use corrugated or divided tray bottoms to separate 
products into rows, placing them in alignment for packaging, or 
some other process. 

Another practical and common use is in conjunction with mag- 
netic separators where ferrous materials must be separated from 
other products. A good example is where a vibratory feeder is 
placed directly in front of a magnetic pulley and as the material 
flows from the feeder, the nonmagnetic portion falls straight down 
and the ferrous content is drawn through the air to the pulley 
face. 

Vibratory feeders are probably the best method of conveying 
materials to and from the hoppers of batch-weighing scale ma- 
chines and are also used for feeding processing mills, pulverizers, 
and driers. As an integral part of packaging machines, they are 
used to accurately feed and control flows of macaroni, noodles, 
pills, candy, cookies, potato chips, hardware, and a multitude of 
other products. 

They can be used advantageously for conveying where space 
prohibits the use of belt-type conveyers and for applications where 
several materials must be fed simultaneously in different pro- 
portions for formulating, blending, and mixing. 

Other common uses are for overhead dusting of countless other 


(7) Leaf spring (Fiberglass-epoxy resin) 


Fig.5 Cross-sectional drawing of bin vibrator showing enclosed electro- 
permanent magnetic motor and leaf-type springs 
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products. In this sort of an application it is often possible and 
practical to place the feeder at right angles to the conveyer carry- 
ing the product to be dusted, and by means of long diagonal dis- 
charge-type trays, spread the dusting medium over wide areas. 


Application Engineering, Vibratory Feeders 

Engineering of vibratory feeders to individual applications is 
not difficult, especially if the manufacturer’s standard products 
are used and his suggestions are followed. Particle size, weight, 
and shape, in addition to cohesive and adhesive characteristics, 
and the moisture content of materials, combine to influence and 
determine the rate at which a product can be fed. 

To make an intelligent appraisal and recommendation for any 
application, the following information is required: 


1 Description of material to be fed 
(a) Weight per cubic foot 
(b) Particle size 
(c) Adhesive or cohesive characteristics 
(d) Moisture content 
2 Volume of material to be fed 
3 Manner in which material is to be introduced to tray 
4 Type of tray desired 
(a) Flat pan (open or covered) 
(b) V or half-round 
(c) Tubular 
(d) If for screening, complete specifications including 
mesh size 
5 Description of any connections or attachments 
6 Voltage and cycles of power supply (units will not operate 
on direct current) 
7 Any other pertinent information 
(a) Extreme temperatures or unusual atmospheres 
(b) Method of discharge (alignment, spreading, dusting, 
and so forth) 
(c) Any unusual electrical control 


Accurate “Feeder Capacity Charts” are nonexistent; however, 
rule-of-thumb guideposts have been established for a wide range 
of materials. Laboratory testing of samples of special materials 
is recommended. 

However, even with laboratory tests, it is not always possible 
to predict performance accurately. Some materials assume 
varying characteristics under field conditions as a result of changes 
in humidity or temperature. Variations in column-loading con- 
ditions need also be considered. 

Another feeding aspect that should be recognized is that the 
output from a feeder can be no more consistent than the input to 
the feeder. If materials flow in surges from a hopper into a 
feeder, they will probably be discharged in the same pattern— 
especially with tangled materials or materials with a steep angle 
of repose. Free-flowing products have a tendency to level out. 

Another factor which can influence feed rates is the angle of the 
feeder tray. If the feeder slopes downward, the output will gen- 
erally increase about 5 per cent for each degree of slope; whereas 
if it is inclined upward the output will usually be decreased by 
approximately 5 per cent for each degree (maximum 15 de- 
grees). Tray materials and surface finishes are important for 
maximum efficiency in feeding certain products. However, the 
feeder manufacturer can usually recommend the most satis- 
factory finish. 

Column loads—the weight of products in bins or hoppers that 
actually bears upon the tray—are more serious and detrimental 
to the performance of some vibratory feeders than with this new 
type, which is stable over a wide range of loading. With well 
designed bins or hoppers the material discharge will be baffled or 
deflected in such a way that column loading will be of no conse- 
quence. But where the full column load bears directly upon the 
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feeder tray, its weight should be calculated and added to the 
normal product load being fed, to make certain that the feeder 
is not overworked. 

It is to be expected that, when tubular, half-round, V-shaped, or 
special tray configurations are used, reductions in maximum out- 
put will result. When there is doubt, the best way to determine 
feed characteristics of a material is to submit a sample for labo- 
ratory testing. This is especially true where problems of dusting, 
de-dusting, separating, or alignment are involved. Occasionally, 
solutions to problems can best be found by actually testing the 
feeder under actual working conditions in the user’s plant. 


Some of the most difficult. performances to predict with vibra- 
tory feeders are those encountered in problems of screening. 
Because customer objectives and product behaviors are so varied, 
it is recommended that screening problems be referred to the 
feeder manufacturer for recommendations. Laboratory testing 
of samples may be desirable. 


Special Feeders 


It would be very nice if we could always deal with standard 
products, but this is not possible or practical. Where variations 
from standard are necessary, let’s look at some of the factors that 
should be considered. 

The most frequent variation from standard products occurs 
either in tray size or configuration. Any major departure from 
established standard dimensions or weights could result in poor 
feeder performance unless cartain basic rules of tray design are 
observed. Certain departures can be tolerated. A common 
rule of thumb is that usually a tray may be a maximum of 50 
per cent wider or 50 per cent longer than the largest standard 
size; but not both. The heaviest standard tray or tray load may 
also be exceeded by 50 per cent but it is to be expected that in the 
case of such extremes the total volume output would be some- 
thing less than maximum. Tests are often necessary to check 
performance. 

In the design of vibratory feeders, tray size, weight, and weight 
distribution are important. For best performance, the center cf 
gravity of the tray, or moving mass, should be as nearly as pos- 
sible in line with the line of drive and the center of gravity of the 
reaction mass, or feeder body. It has also been learned that 
the weight of the reaction mass must be from two to four times the 
weight of the moving mass. If a tray is off balance, it can result 
in poor performance because of improper location of the center of 
gravity. If tray and load are overweight and the tuning spring 
stiffness has been increased to regain maximum displacement, the 
weight ratio can be so disturbed that there will be excessive move- 
ment of the reaction mass. The resulting increased relative 
movement between the tray and body can cause excessive de- 
flection and destruction of the springs. 

Feeder trays are usually made of light gage materials to mini- 
mize weight. It is important that flat surfaces be rigid to pre- 
vent tray breakage and spurious vibrations which can cause non- 
uniform movement of the particles on the tray. Therefore 
reinforcing must be added to adequately stiffen the tray without 
adding excessive weight. 

Either flexible or rigid connections or additions to feeders or 
trays, which could restrict or otherwise influence the carefully 
engineered mass-spring relationship, should not be made without 
expert engineering study and advice. Unauthorized modifica- 
tions are a problem with all vibratory-feeder manufacturers. 

Another problem of major concern to manufacturers is the po- 
tential user who, usually for economic reasons, but sometimes 
because of space limitations, desires that a certain size vibratory 
feeder be doctored up in some way to do the work of the next 
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larger size model. Although these products have a built-in 
safety factor, insistence on overworking them, even with the 
hopeful co-operation of an application engineer, often results in 
disappointing and expensive failures. 

A good engineer would certainly not expect a 1-hp motor to 
do a 2-hp task, nor would he expect to continually run most ma- 
chinery at wide-open speeds without experiencing serious dif- 
ficulties. It is far better to select a unit that will safely ac- 
complish a given task rather than to hope that by some stroke of 
good fortune a smaller unit might just happen to do the job. 
The manufacturer’s suggestions should be followed because his 
recommendations will be based on extensive experience with 
many products and applications. 


Application Engineering, Unit Vibrators 


Because of variations in bin designs and material character- 
istics, it is only possible to furnish approximate rules for the se- 
lection of unit vibrators. Nevertheless, one can intelligently 
engineer individual applications by recognition of the most im- 
portant factors that influence vibrator performance. 

Let’s take the design of the bin itself. The slope of the sides 
should always be greater than the angle of repose of the materials 
to be handled; also the discharge opening should be as large as is 
practically possible. In the case of large lumps, the opening 
should be at least 21/2 to 3 times the diameter of the largest 
piece. With fine materials, the opening should be designed, if 
possible, to minimize the possibility of bridging. A gate inter- 
lock switch should be provided so that vibrators cannot be 
turned on while the gate is closed. 

Bin wall thickness is of importance because unit vibrators are 
scientifically designed in a range of sizes so that each size is cal- 
culated to cause maximum safe deflection of a specific bin wall 
thickness. Except in unusual cases, unit vibrators can be con- 
fidently selected to operate on the bin wall thickness for which 
each is designed. Printed literature always recommends sizes 
for suggested wall thickness. 

For best results on large bins with light-gage walls, there are 
two alternatives: 


(a) Selection of a larger vibrator than the one specifically 
designed for that thin material. By welding a reinforcing plate 
to the bin wall to distribute the forces over a greater area, mini- 
mizing the possibility of fatigue cracks, a good effective instal- 
lation is possible. 

(b) Install several smaller units around the bin to spread 
vibration over all of the bin surfaces. (This alternative is usually 
more effective). 


In all cases involving light-gage bin walls, reinforcing pads 
should be welded or bolted in place to prevent bin damage and to 
provide a good mounting surface for the vibrator. 

No factor could be more important than correct location and 
mounting of vibrators on the bins. There are no hard and fast 
tules, but ordinarily a vibrator should be placed one third to one 
quarter of the distance up from the discharge opening on the 
vertical center line of a panel that is free to vibrate. It should 
not be close to a reinforcing member. On conical hoppers or 
bins, channel-type adapters should be used. In all cases the 
bolts should be drawn very tight and should be checked occa- 
sionally. To prevent any possibility of fatigue and breakage, 
heat-treated or alloy-steel bolts may be used. 

Moving a vibrator a short distance upward or downward some- 
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times makes the difference between failure or success of an in- 
stallation. Frequently the addition of another vibrator to the 
opposite side of the bin at a different height is advisable. 

When more than one vibrator is attached to the same bin, care 
should be taken to connect the wiring identically. Connect the 
same color coded wire to the same side of the power supply so 
multiple units will be properly phased and will work as a team, 

It is recommended that in many applications an oversized vi- 
brator be used .. . then by using a control, adjust its amplitude 
to variable operating conditions. This not only prolongs the 
life of the unit, but it also reduces the noise level. 

In some applications where fine materials are being handled, 
static electricity has been a major problem. In some of these 
applications, grounding of the bin has provided a solution. It 
is the nature of static electricity that it must be grounded as 
near its source and in as many places as is practically possible. 
Thin strands of copper wire soldered at 3-ft intervals and to the 
nearest column or beam can solve this kind of problem. 


A complete line of controls has been developed for controlling 
the voltage supply to all models of vibratory feeders and unit 
vibrators. The controls include either a potentiometer or a 
variable autotransformer, a switch, and for some models capaci- 
tors for power-factor correction. By means of these controls, 
volumes can be regulated from a mere trickle to a rushing tor- 
rent. 

One must remember that feeder output is almost directly pro- 
portional to the voltage supply, and that a voltage fluctuation of 
plus 10 per cent is not uncommon. 

Installation of a voltage regulator can reduce this variable to 
plus or minus 1 per cent. : 

Frequently in packaging and automatic scale operations, it is 
desirable to have a feeder operate alternately at high and low 
speeds. Dual control boxes are available for this purpose. Each 
control contains two potentiometers or autotransformers and a 
provision is made for connection of an external relay, by the cus- 
tomer. This relay is used to switch from control by one po- 
tentiometer to the other. One potentiometer is set to a high 
value for “volume” feed and the other at a low value for a 
“trickle” feed. Each of these potentiometers is set by the cus- 
tomer to the relative values desired. 


All of the products mentioned can be made to operate effec- 
tively on 50-cycle current. 

All of the products, including controls, have the acceptance of 
the Mill Mutual Fire Prevention Bureau for installation in 
hazardous, dusty areas where Class II, Group G equipment 
would ordinarily be used. They have also been approved by the 
Canadian Standards Association for General Purpose. Under- 
writers Laboratory explosionproof models have not been de- 
veloped. 

With the great emphasis being placed today on cost reduction 
by automated control in materials handling, it is desirable that 
engineers be well acquainted with one of the most acceptable, 
proven methods of accurately controlling flows of bulk ma- 
terials. It is hoped that the highlights mentioned in this article 
may provide a better understanding of the problems and many 
valuable uses of vibratory products. 


Transactions of the ASME 


a 
t 
t 
Controls 
t 
i 
I 
I 
f 
] 
{ 


SME 


J. R. BAUMGARTEN? 
R. COHEN® 
A. HALL, JR.* 


A Mechanical Time-Delay Device 
Sensitive to Centrifugal Fields 


A proposed design of a mechanism suttable for use as a time-delay device is given. 
The equations of motion of the device are derived and the effects of variation of certain 
design parameters are studied. 


Introduction 


IME-DELAY DEVICES have wide usage in industrial 
and military applications. The most common device used in 
timing events in a mechanical or electrical system is an electronic 
time-delay circuit. There are applications in the control field 
where electrical power is either not available or must be used for 
other purposes. Mechanical power is usually available in indus- 
trial applications in the form of shaft power and is often available 
in missile and other ballistic applications in the form of a rotating 
casing. In addition to the unavailability of electric power, a 
mechanical device may be chosen over the electrical device be- 
cause of the possibility that a particular design will be less ex- 
pensive, easier to fabricate, and more rugged and reliable than the 
corresponding electrical one. 

It is the purpose of this paper to present a device which does not 
require electrical power. This device derives its power from a 
centrifugal field. The idea underlying the mechanism described 
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Nomenclature 


in the following was originally suggested by Mr. W. R. Benson, 
Feltman Research and Engineering Laboratories, Picatinny 
Arsenal, Dover, N. J. 


Description of the Mechanism 

One possible design of a mechanical time-delay device (Fig. 1) 
consists of two thin cylindrical disks, J: and A ball bearing is 
placed between the disks and the disks are free to rotate relative 
to each other. A string (or some other tension member) is 
fastened to the periphery of the outer disk at point A. The string 
is wound about the periphery of the inner disk and a small mass 
m is fixed on the string at an intermediate point between the 
disks. A small helical spring K produces a torque proportional 
to the relative rotation of the two disks. The spring is of small 
mean radius and it is assumed that centrifugal forces of the spring 
have little effect on the system. A pin P locks the two disks to- 
gether with the system in its initial configuration. 

Initially the string is taut and an initial torque due to the 
preset in the coil spring resists relative rotation of the disks. The 
pin is in place and the disks are locked together. The assembly is 
set rotating at a constant angular velocity n and the pin is re- 
moved, subjecting the device to a centrifugal field. The mass 
moves outward as a result of the centrifugal field, causing the 
string to unwind and a relative rotation results between the two 
disks. 


The Equations of Motion 


Using 62 and ¥ as generalized co-ordinates (Fig. 2), the angular 
position of the small disk, J;, can be expressed in terms of the 
generalized co-ordinates as follows: 


a = radius of the arm, a 
b = initial string length from mass to 


outer disks, J; and J2, from 
fixed reference 


F= a = dimensionless inner radius 


int of tangency, a + b = ‘ c 
pe sin~ a = constant G = = dimensionless moment of 
string : p = angle subtended by string and tite 
c = radius of small disk J, 1 perpendicular to the arm, a md? 
d = radius of large disk I; ¢ ‘ 3 H = — = dimensionless small mass 
m = small mass B= im dimensionless spring 
n = initial rotational velocity ore I, Iz = moment of inertia, inner and 
t = time C,,C: = damping constants on disks J; cuter: 
+ = horizontal co-ordinate of the and I; J = — = dimensionless radius of the 
mass, m Cc d 
y = vertical co-ordinate of the mass, D = — = dimensionless damping arm, @ 
m In K = spring constant 
B = initial preset angle Psorngs R = instantaneous radius of the 
Y = instantaneous angular position ae . mass, m 
of the mass, m * In a T, = bearing frictional torque 
6:, 02 = angular position of inner and constant T. = KB = initial preset torque 


Journal of Engineering for Industry 


NOVEMBER 1959 / 398 


n in- 

the 

care 

ly so 

m, 

itude 

3 the 

dled, 

these 

ag 

sible. 

o the 

slling 

unit 

or 

paci- 

trols, 

tor- 

pro- 

on of 

le to 

it is 

1 low 

Each 

ind a 

> cus- 

@ po- 

high 

effec- 

ice of 

ment 

y the 

nder- 

n de- 

iction 

that 

table, 

ma- 

rticle 

many 

_ 


Fig. 2 Notation used in derivation 


=~ VR? -c? 


+ (a+ b — acos y) 
R2 


b 


where R is the instantaneous radial distance of the mass from the 
center of rotation and 


R? = a? + d*? — 2ad cos (y + pw) (2) 


Assuming that the mass of the string is small in comparison to 
the mass, m, the kinetic energy of the system is 


1 
T=— m(6.2d? + a%(62 + 7)? — + 7) cos (y + u)] 
1624 1 (3) 
2 2 
Assuming that gravitational potential is small in comparison 
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with spring potential, and neglecting any elasticity in the string, 
the potential energy of the system is then 


K(6, — 6)? + — 63) (4) 


Energy is dissipated due to shear of the air surrounding the 
disks. Assuming this shear to be viscous, there exists a dissipa- 
tion function’ from which this shear torque can be derived. This 
is 


(C62 + (6) 


In addition to viscous friction, there is present energy dissipa- 
tion due to the friction in the ball bearing. This frictional 
torque opposes relative rotation and is assumed to be constant in 
magnitude with a direction such as to oppose relative rotation. 
The work done against bearing friction is 


W = _7, (0; — (6) 
Y 
Forming Lagrange’s equations for the two generalized co-ordi- 
nates, the equations of motion for the system are 
E cos? p + mae | + cos p + mac — med cos (y + 
2ad sin (y + pu) cos p 
a(R? — c?) — ad? sin? (y + 


—T,cosp— (Ci cor + = eos) 


= —h cos p [sin p- 


cos p — — 62) cos p+ sin (y +m) (7) 


and 
—¥ [mad tan p sin (y + + + md?(tan p cos (y + 4) 
+ sin (y + w)) sin (y + mw) — mad tan p sin (y + p)] 
= md6,2 sin (y + — ¢ sec p) — + 


—mad(b. + sin (y + + — 2) + Te (8) 


where 


sin p = {od sin (y +m) — [a —d cos (y + WR? 


and 


cos p = fela—deos (y + + dsin (y + 
(10) 


It is now advantageous to redefine the independent variable 
time, in the form of dimensionless time. Recalling that the initial 
velocity of the assembly is n, a new variable 7 is defined as 

T=nt (11) 


In making this change of variables, let the notation indicating 
derivatives with respect to T be 


5 Herbert Goldstein, ‘‘Classical Mechanics,” Addison-Wesley 
Publishing Company, Inc., Cambridge, Mass., 1953, p. 21. 
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4) (7) 


dy _ 1 dy » 1 


Using this change of variables, Equations (7) and (8) are re- 
formed in terms of dependent variables and dimensionless 
parameters. Hence, 


D 
(62’)? sin (y + — C08? p 


T, 


2J sin (y + 4) cos p 
R?/d? 
J(R? — + J sin? (y + 


R2/d?2 (R? — c?)/d? 


B B 
py (1 — 62) cos p — cos p) 


G 


+ 


FJ 
E 
[ wo sin (y + pw) (1 seep) 62! 


T, 
JI? 


BB i 
7 (0, — 02) + A (18) 


— H(6.' + sin (y + w) + 


and 
F 


E T 


B B H 
(0: — 62) + | Wy) - (4 (62)? sin (y + 


T, 


B 


2J sin (y + 4) cos p 
R*/d? 


J(R? — c*)/d* + J sin? (y + 
— 


where 


Woy) cost p +H (15) 


H 
= (16) 


Y(y) = —H tan pein (y + 4) (17) 


Zy) = z+ — sin (y + »)[tan p cos (y + 


J 
+ sin (y + w) — J tanp] (18) 


and 
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A = W(y)Z(y) — (19) 


Equations (13) and (14) are the equations of motion of the system 
where Equation (1) gives the relative rotation as a function of . 
Equations (13) and (14) were solved simultaneously using the 
Runge-Kutta method of numerical quadrature, generalized for 
treatment of simultaneous equations. It should be noted here 
that zero values of the determinant A were not encountered in 
the solution. 


The Solution 


The equations of motion, Equations (13) and (14) were pro- 
grammed for digital computation. A basic design was selected for 
the time-delay mechanism and the solution obtained. Dimen- 
sions of the basic mechanism are given in the Appendix. 

The effect of variation of the design parameters on the response 
of the system was studied. Holding all other dimensionless ratios 
constant, the response of the system was noted as a given ratio 
was varied. The effects of variation of the dimensionless spring ~ 
constant B, the dimensionless inertia G, the dimensionless radius 
F, the dimensionless mass H, and the initial preset 8 were studied. 

In general, it was found that the response of the system is os- 
cillatory in nature. A typical solution for such a case (basic di- 
mensions except for 0.4Ho) is shown in Fig. 3 as the solid line. It 
is noted that for this design a peak relative rotation of 153.4 deg 
occurs. Half this relative rotation, namely, 76.7 deg, is reached 
at 0.0265 second. Let this point be referred to as the half- 
amplitude point. This point is used as a criterion for compari- 
son of different cases studied. 


200 - 
—— NUMERICAL SOLN. 
-- APPROX. SOLN. 
150 
w 4 
w 4 
vA AMPLITUDE 
W100 
| HALF-AMPLITUDE 
; 
0.02 004 0.06 008 0.10 014 


TIME - SECONDS 
wa Typical response curve for a periodic case (basic system except 


In those systems where large relative rotation occurred, it was 
noted that the constraint forces in the string became compressive 
at a given point. This is incompatible with the design and the 
solution is valid only to this point. Let this case be referred to as 
the nonperiodic case. A typical solution for such a case (basic 
dimensions except for 0.9 H) is shown in Fig. 4. For the non- 
periodic case, the point of maximum slope (point of maximum 
relative velocity) of the rotation-time curve was calculated. This 
point of maximum relative velocity was taken as the half-ampli- 
tude point of the nonperiodic cases for the purpose of compari- 
son with periodic cases. This point occurs at 382 deg and 0.0400 
sec for the curve of Fig. 4. 

The half-amplitude and the time required to reach the half- 
amplitude are the basis for comparison of the cases studied. These 
half-amplitude points are shown plotted versus time in Fig. 5. 
Those points illustrating periodic cases are shown as circles; 
those points illustrating nonperiodic cases are shown as squares. 

6 ¥. B. Hildebrand, “Introduction to Numerical Analysis,” Mc- 
Graw-Hill Book Company, Inc., New York, N. Y., 1956, p. 237. 
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Fig. 5 Holf-amplitude versus time for 
parameter variation. Curves labeled g in- 
dicate G variation; those labeled b indi- 

% cate B variation; those labeled ¢ indicate 
variation in F, H, and n. 


8 


HALF-AMPLITUDE DEGREES 


8 


Curves have been drawn through points which represent similar 
cases, that is, cases for which only one particular parameter has 
been varied. For example, the curve labeled g; illustrates the 
effect of variation of the dimensionless inertia G Particular 


values of the varied parameters are shown next to the appropriate sii | B05 

point on the curve. The point representing the basic system is | 

noted as “NORM” in Fig. 5. | 
Fig. 6 represents the change in response for systems of varying ade 


initial preset 8. The curve represents half-amplitude points for 
systems with 1.5 times the initial spring rate, a 100 G inertia 
ratio, and with the remaining basic parameters. 


An Approximate Solution 


It is seen from Fig. 3 that the solution to the equations of mo- 
ti on approximates harmonic motion for the periodic case. Thus 


HALF-AMPLITUDE DEGREES 
6 
fo] 


for this case, it can be shown that the system can be approxi- 50 AS 
mated by a second-order linear equation having the solution, 
6: — 62 = A(1 — cos wi) (20) 6 
° 0.05 010 0.15 0.20 028 0.300 
Equation (7) is rewritten with 6; set to zero and with friction and _TIME- SECONDS 
damping neglected. The mac ¥ term must be negligible and two 


. 6 Half-amplitude versus time for variation of preset angle § 
other relations are implied for the approximation. These are ys 


(basic system except 1.5 Bo, 100 G,) 
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where A is the half-amplitude of the approximate solution and 6, 
js the assumed constant velocity of the outer disk. Equation (22) 
can be solved graphically for the approximate half-amplitude. 
The relation for A as a function of y is given by Equation (1). 
The left and right-hand sides of Equation (22) are plotted versus 
y for a given case. The intersection of the two curves yields the 
approximate half-amplitude. Fig. 3 shows a comparison between 
the actual solution and the approximate solution (dashed curve). 


Conclusions and Discussion 


Referring to Fig. 5, it is seen from curves gi, ge, and gz that 
changes in the dimensionless inertia markedly affects the time 
required for the mechanism to reach half-amplitude. This 
parameter has little effect on the magnitude of the relative rota- 
tion of the mechanism. Curve ge indicates the response for sys- 
tems designed with the basic parameters except for 1.5 Bo. 
Curve g; indicates the response for systems with one half the 
basic mass and varied inertia ratios. 

Curves b; and be indicate the response of the systems having 
the basic parameters and varied spring rates. It is seen that the 


value of the dimensionless spring constant markedly affects the _ 


magnitude of responses. For low values of inertia, the spring 
rate has little effect on the time of response. 
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Curves c; and cz of Fig. 5 indicate the effect on response for 
changes in n, the initial velocity, F, the dimensionless radius, and 
H, the dimensionless mass. For low values of inertia, variation 
of these parameters has little effect on time of response. At high 
values of inertia, these parameters markedly affect both magni- 
tude and time of response. 

Fig. 6 indicates that small changes in the initial preset 6 affect 

the magnitude of response only. Large changes in preset affect 
both magnitude and time of response. It is suggested that per- 
haps this is the easiest parameter to vary in the assembled mech- 
anism. 
The designer is given two tools. The first are the equations of 
motion, Equations (7) and (8), and the composite response curve 
of Fig. 5. The second tool is the approximation outlined in 
Equations (20), (21), and (22). The approximate solution is not 
intended as an accurate means of selection of parameters, but 
only an aid in design. Calculations indicate that the approximate 
solution yields results to within 10 per cent error only for systems 
with half-amplitude of 125 deg or less. Because of the assump- 
tions in the derivations a designer is cautioned to use these re- 
sults only as a guide. This study could be supplemented with 
experimental effort in the light of these assumptions. 


APPENDIX 
The dimensions of the basic mechanism are: 
@ = 3.375 in. m = 1.661399 xX 10-* lb sec?/in. 
c = 0.75 in. I, = 7.897587 X 10-5 in lb sec? 
d = 40 in. I, = 2.3639849 X 107? in lb sec? 
K = 0.3 in. lb/rad mo = 2000 rpm 
B = 0.5 rad 
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aie Semi-Conductor Power-Conversion Equip- 
ment, by J. J. Rheinhold 
59—MD-5 The Design and Application Belts, 


and | Gears, 
id W. C. Raridan 


se—MD-6 Clutch, Fluid Coupling, Hydraulic Torque 
considerations and Per- 
by R. W. Bachmann 


59—MD-7 Latest to Resist 
Chemical Corrosion, by S. W.S 


59—MD-8 Recent Materials and 
Finishes to Resist Atmospheric Corrosion, by E. B. 
Friedl, L. J. Nowacki, and W. H. Sefranek 


Systems in Control Applications, by 


Rt ghd Logic Circuits for Machine Control, by 
. V. Weir 


59—?4D-11 Design Curves for Journal Bearings, by 
D. Hays 


oe. Which Bearing and Why? by A. O. 


13 Graphic Method for 
Organization Analysis, by E. M. Ramberg and R. P. 
inic 
59—MD-14 The Man and the Organization, by M. 
M, Bruce 


of the PRODUCTION ENGINEERING CON- 
FERENCE 
Available until March, 1960 


59—Prod-1 The Influence of Surface Residual Stress 
on Fatigue Limit of Titanium, by E. C. Reed and J. A. 
iens 


59—Prod-2 Shear Spinning, by B. N. Colding 


59—Prod-3 Shear-Zone Size, Compressive Stress, and 
Shear Strain in Metal Cutting and Their Effects on 
Mean noggin Stress, by Dimitri Kececioglu 


59—Prod-4 1958 Review of Metal Processing Litere- 
ture Grinding) by George Reichenbach 


59—Prod-5 1958 Review Metal 
ture (Plastic Working), by Boulger 


59—Prod-6 1958 Review of Metal 
(Metal Analysis), Gelimberti 


S. Hahn, H. J. Siekmann, and 6M jomsen 


habe Management Implications of Production 
Controls. . Planning a New Plant, by V. C. Peterson 


cet Product Reliability Measurement, by D. F. 


Electrical Discharge Applied 
i ing for Pr lolm 


idea 0 Process Analysis Via Simulation, by R. W. 
Metzger 


59—Prod-11 The Organization of Central Stores and 
Diep Systems for Productivity. . .A Solution to 
Some of the Problems of Job Shop Manufacture, by 
W. P. Dugan 


59—Prod-12 Chipless Forming of Toothed Parts by 
Rolling, by H. Pelphrey 


; 80¢ per copy to nonmembers 
40¢ per copy to ASME members 


Published by ASME 


by E. S. CL. . 


of the MAINTENANCE AND PLANT EN- 


29 W. 39th Street, 
New York 18, N. Y. 


MATERIALS HANDLING HANDBOOK 


In this Handbook has been condensed and organized virtually the entire literature of modern mates 
rials handling, stressing the economies and short-cuts of the greatest importance throughout industry, 
In consequence, it is the first complete, definitive guide to principles and procedures which are adaptable 
to all types of industrial situations. 

The Handbook gives actual operating experiences of hundreds of companies in all segments of manus 
facturing and industry. It offers advice on choosing the equipment to do the job, timing and coordinating 
job movements, factory planning, production control, industrial packaging, yard handling, warehousing 
truck-rail-marine- and air-terminal handling, etc. It explains the cost methods and includes forms fc 
determining the cheapest shipping methods for varying loads and destinations. 

Supplementing the text discussion are photographs, diagrams, and time-saving picturizations whi 
show typical installations of specialized and general purpose apparatus, innovations in machinery and 
design, standard and unusual equipment applications, and model layouts and floor plans. ’ 

In short, the Materials Handling Handbook provides the means for making a scientific efficiency! 
check on every phase of operation—from hauling and storing raw materials to movoment through 
plant, packaging, stock piling, and safe —y to the buyer. 


1750 Pages 961 Photographs, Line Drawings, and Tables’ 
68-Page Index Flexible Binding $20.00* 


SAFETY CODE FOR CONVEYORS, 
CABLEWAYS, AND RELATED EQUIPMENT 


is your guide for the safe construction, elements of design, their instal 
lation, operation, and maintenance. 


Covering sixty-six different types of conveyors, this standard gives the terms and definitions of their | 
parté’and components. . . .the safety requirements that apply to the design features of conveyors, such 
as belts and chains, pulleys, sprockets, blocks, bearings, backstops and brakes, and overload-protec-3 
tion. . . the recommendations for installation, inspection and maintenance, for safety in connection with 
related structures, for safe operating practices, and for prime movers and controls. Additionally the’ 
code provides specific safety regulations for aerial cableways, slat conveyors, belt conveyors, pneu- 
matic conveyor systems, lsucket conveyors, chain conveyors, roller and wheel conveyors, live roll conveyors) 
vertical conveyors, screw conveyors, and aerial tramways and the traffic control system for them. 
B20.1-1957 $2.00% 


in the industry is . . . 


SAFETY CODE FOR CRANES, DERRICKS, AND HOISTS 


This code applies to the construction, installation, inspection, maintenance, and of cranes 
and derricks driven by steam engines, electric motors, or internal combustion engines to their runweysy 
to simple drum hoists of whatever motive power; to overhead electric hoists and their runways; and 
overhead air hoists; also to hand powered derricks. For every type of equipment included generel 
rules are given covering welded construction, foundation anchorage, cages and machine houses, 
lubricating devices, mechanical guarding, electric equipment, brakes, and runways. On inspection, 
testing, and maintenance, rules are provided for parts subject to wear, locking devices, brakes, chain 
and chain drives, cleaning, and repair. On operation, the recommendations apply to the condu¢ of 
and qualifications for operators, the safe working load of the equipment, signals, brake testing, firing 
boilers, etc. Specific rules included are for overhead and gantry cranes; locomotive, crawler, and motor 
truck cranes; derricks; simple drum hoists; overhead electric hoists and base-mounted electric hoists; 
and air hoists. 


B30.2-1952 


$2.50° 


_* 20% discount to ASME members 


5S9—MPE-5 En Preventive Main- 
Available until. March, 1960 
SP—-MPE-1 The Technique of Efficient MATERIALS. 
59—MH-1 Coal-Haendl Facilities at Power 


59—MPE-3  Practi Protective-Coatin 59—MH-2 Mechanized d Freight Handling et Railroad § 


59—MPE-4_ Maintenance of Parking Lots and Grounds, 
by A. K. Carter 
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